Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,006,415 Δ1 bp coding (464/543 nt) zapC → FtsZ stabilizer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,006,4130G.98.0% 205.8 / ‑3.0 51coding (462/543 nt)zapCFtsZ stabilizer
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (26/24);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTG  >  NC_000913/1006281‑1006545
                                                                                                                                    |                                                                                                                                    
cGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACg                                                                                                                               <  2:1568183/139‑1 (MQ=255)
  aTGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGcc                                                                                                                             >  1:1172368/1‑139 (MQ=255)
        ccTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTaaa                                                                                                                       >  1:243390/1‑139 (MQ=255)
                                     cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAGCGa                                                                                                            >  2:981546/1‑121 (MQ=255)
                                     cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGa                                                                                                            <  1:981546/121‑1 (MQ=255)
                                     cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGt                                                                                          >  1:105504/1‑139 (MQ=255)
                                     cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGt                                                                                          >  1:287730/1‑139 (MQ=255)
                                       ggTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTa                                                                                        <  2:1172368/139‑1 (MQ=255)
                                        gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTaaaa                                                                                                                      >  2:1065827/1‑108 (MQ=255)
                                        gTCAATCTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTaaaa                                                                                                                      <  1:1065827/108‑1 (MQ=255)
                                          cAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATg                                                                                     >  2:755390/1‑139 (MQ=255)
                                                  tGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGc                                                                             >  2:407213/1‑139 (MQ=255)
                                                  tGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGc                                                                             >  1:140387/1‑139 (MQ=255)
                                                    gTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTGATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCtt                                                                           <  1:407213/139‑1 (MQ=255)
                                                                  cGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGa                                                                  >  2:62168/1‑134 (MQ=255)
                                                                  cGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGa                                                                  <  1:62168/134‑1 (MQ=255)
                                                                               cTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTTGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACgcgc                                                <  2:1431360/139‑1 (MQ=255)
                                                                                   gCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAAc                                            >  2:1139279/1‑139 (MQ=255)
                                                                                        ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAg                                       >  2:581392/1‑139 (MQ=255)
                                                                                          ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGa                                                                                 <  1:214353/95‑1 (MQ=255)
                                                                                          ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGa                                                                                 >  2:214353/1‑95 (MQ=255)
                                                                                          ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTg                                     >  2:640990/1‑139 (MQ=255)
                                                                                            cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGAc                                                                                                           <  1:579013/67‑1 (MQ=255)
                                                                                            cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGAc                                                                                                           >  2:579013/1‑67 (MQ=255)
                                                                                            cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCa                                   >  1:827076/1‑139 (MQ=255)
                                                                                            cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCa                                   >  2:1576915/1‑139 (MQ=255)
                                                                                            cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCa                                   <  1:640990/139‑1 (MQ=255)
                                                                                              caACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGt                                 <  1:581392/139‑1 (MQ=255)
                                                                                                aCCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCt                               >  2:127389/1‑139 (MQ=255)
                                                                                                 cccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCg                                                                                               >  1:317195/1‑74 (MQ=255)
                                                                                                 cccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCg                                                                                               <  2:317195/74‑1 (MQ=255)
                                                                                                    tGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCa                                                                                                                            >  1:815905/1‑42 (MQ=255)
                                                                                                    tGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCa                                                                                                                            <  2:815905/42‑1 (MQ=255)
                                                                                                       gttgttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTac                        <  2:105504/139‑1 (MQ=255)
                                                                                                         tgttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTacac                      <  2:140387/139‑1 (MQ=255)
                                                                                                         tgttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTacac                      <  2:243390/139‑1 (MQ=255)
                                                                                                           ttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAg                    >  1:1248691/1‑139 (MQ=255)
                                                                                                              ttGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTa                 <  2:1643305/139‑1 (MQ=255)
                                                                                                                 cGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAAcgc                                                  >  1:11950/1‑103 (MQ=255)
                                                                                                                 cGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAAcgc                                                  <  2:11950/103‑1 (MQ=255)
                                                                                                                   ggTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                                                                                   <  2:1584454/68‑1 (MQ=255)
                                                                                                                   ggTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                                                                                   >  1:1584454/1‑68 (MQ=255)
                                                                                                                        cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTaa                                                     <  2:333945/93‑1 (MQ=255)
                                                                                                                        cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTaa                                                     >  1:333945/1‑93 (MQ=255)
                                                                                                                        cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTa                 >  1:989640/1‑129 (MQ=255)
                                                                                                                        cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGaa       >  1:262046/1‑139 (MQ=255)
                                                                                                                        cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGaa       <  1:127389/139‑1 (MQ=255)
                                                                                                                        cgcgATGCAG‑TGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTa                 <  2:989640/129‑1 (MQ=255)
                                                                                                                           gATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATgg    <  2:827076/139‑1 (MQ=255)
                                                                                                                            aTGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGt   <  2:1248691/139‑1 (MQ=255)
                                                                                                                             tGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCCAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTg  <  1:1139279/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                    
CGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTG  >  NC_000913/1006281‑1006545

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑