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
                                                                                                                                    |                                                                                                                                    
cGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACg                                                                                                                               <  2:1568183/139‑1 (MQ=255)
  aTGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGcc                                                                                                                             >  1:1172368/1‑139 (MQ=255)
        ccTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTaaa                                                                                                                       >  1:243390/1‑139 (MQ=255)
                                     cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAGCGa                                                                                                            >  2:981546/1‑121 (MQ=255)
                                     cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGa                                                                                                            <  1:981546/121‑1 (MQ=255)
                                     cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGt                                                                                          >  1:105504/1‑139 (MQ=255)
                                     cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGt                                                                                          >  1:287730/1‑139 (MQ=255)
                                       ggTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTa                                                                                        <  2:1172368/139‑1 (MQ=255)
                                        gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTaaaa                                                                                                                      >  2:1065827/1‑108 (MQ=255)
                                        gTCAATCTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTaaaa                                                                                                                      <  1:1065827/108‑1 (MQ=255)
                                          cAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATg                                                                                     >  2:755390/1‑139 (MQ=255)
                                                  tGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGc                                                                             >  2:407213/1‑139 (MQ=255)
                                                  tGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGc                                                                             >  1:140387/1‑139 (MQ=255)
                                                    gTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTGATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCtt                                                                           <  1:407213/139‑1 (MQ=255)
                                                                  cGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGa                                                                  >  2:62168/1‑134 (MQ=255)
                                                                  cGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGa                                                                  <  1:62168/134‑1 (MQ=255)
                                                                               cTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTTGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACgcgc                                                <  2:1431360/139‑1 (MQ=255)
                                                                                   gCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAAc                                            >  2:1139279/1‑139 (MQ=255)
                                                                                        ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAg                                       >  2:581392/1‑139 (MQ=255)
                                                                                          ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGa                                                                                 <  1:214353/95‑1 (MQ=255)
                                                                                          ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGa                                                                                 >  2:214353/1‑95 (MQ=255)
                                                                                          ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTg                                     >  2:640990/1‑139 (MQ=255)
                                                                                            cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGAc                                                                                                           <  1:579013/67‑1 (MQ=255)
                                                                                            cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGAc                                                                                                           >  2:579013/1‑67 (MQ=255)
                                                                                            cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCa                                   >  1:827076/1‑139 (MQ=255)
                                                                                            cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCa                                   >  2:1576915/1‑139 (MQ=255)
                                                                                            cacaACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCa                                   <  1:640990/139‑1 (MQ=255)
                                                                                              caACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGt                                 <  1:581392/139‑1 (MQ=255)
                                                                                                aCCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCt                               >  2:127389/1‑139 (MQ=255)
                                                                                                 cccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCg                                                                                               >  1:317195/1‑74 (MQ=255)
                                                                                                 cccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCg                                                                                               <  2:317195/74‑1 (MQ=255)
                                                                                                    tGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCa                                                                                                                            >  1:815905/1‑42 (MQ=255)
                                                                                                    tGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCa                                                                                                                            <  2:815905/42‑1 (MQ=255)
                                                                                                       gttgttATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTac                        <  2:105504/139‑1 (MQ=255)
                                                                                                         tgttATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTacac                      <  2:140387/139‑1 (MQ=255)
                                                                                                         tgttATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTacac                      <  2:243390/139‑1 (MQ=255)
                                                                                                           ttATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAg                    >  1:1248691/1‑139 (MQ=255)
                                                                                                              ttGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTa                 <  2:1643305/139‑1 (MQ=255)
                                                                                                                 cGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAAcgc                                                  >  1:11950/1‑103 (MQ=255)
                                                                                                                 cGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAAcgc                                                  <  2:11950/103‑1 (MQ=255)
                                                                                                                   ggTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                                                                                   <  2:1584454/68‑1 (MQ=255)
                                                                                                                   ggTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                                                                                   >  1:1584454/1‑68 (MQ=255)
                                                                                                                        cgcgATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTaa                                                     <  2:333945/93‑1 (MQ=255)
                                                                                                                        cgcgATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTaa                                                     >  1:333945/1‑93 (MQ=255)
                                                                                                                        cgcgATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTa                 >  1:989640/1‑129 (MQ=255)
                                                                                                                        cgcgATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGaa       >  1:262046/1‑139 (MQ=255)
                                                                                                                        cgcgATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGaa       <  1:127389/139‑1 (MQ=255)
                                                                                                                        cgcgATGCAGTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTa                 <  2:989640/129‑1 (MQ=255)
                                                                                                                           gATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATgg    <  2:827076/139‑1 (MQ=255)
                                                                                                                            aTGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGt   <  2:1248691/139‑1 (MQ=255)
                                                                                                                             tGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCCAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATGGTg  <  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: