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.96.4% 215.4 / 1.1 56coding (462/543 nt)zapCFtsZ stabilizer
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (25/29);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATG  >  NC_000913/1006279‑1006542
                                                                                                                                      |                                                                                                                                 
ggCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTggcga                                                                                                                               <  1:1298108/139‑4 (MQ=255)
 gCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTggcgac                                                                                                                              <  1:285335/139‑5 (MQ=255)
  cGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACg                                                                                                                            >  2:1624491/1‑139 (MQ=255)
                tttGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATg                                                                                                              <  1:732977/139‑1 (MQ=255)
                     ttGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGa                                                                                                         <  1:946365/139‑1 (MQ=255)
                     ttGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGa                                                                                                         <  2:1430792/139‑1 (MQ=255)
                                   cGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGc                                                                                                                           >  1:1463806/1‑107 (MQ=255)
                                   cGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGc                                                                                                                           <  2:1463806/107‑1 (MQ=255)
                                       cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCt                                                                                                   >  1:1114885/1‑127 (MQ=255)
                                       cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCt                                                                                                   <  2:1114885/127‑1 (MQ=255)
                                          gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCAt                                                                                                               >  1:135387/1‑112 (MQ=255)
                                          gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCAt                                                                                                               <  2:135387/112‑1 (MQ=255)
                                          gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTg                                                                                                  <  2:419346/125‑1 (MQ=255)
                                          gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTg                                                                                                  >  1:419346/1‑125 (MQ=255)
                                                ttGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGa                                                                              <  1:846466/139‑1 (MQ=255)
                                                     ggTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCt                                                                         <  2:1214326/139‑1 (MQ=255)
                                                         gTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAg                                                                     >  2:1570697/1‑139 (MQ=255)
                                                          ttGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGc                                                                    <  2:934705/139‑1 (MQ=255)
                                                                             cGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGa                                                                                                         >  2:247191/1‑83 (MQ=255)
                                                                             cGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGa                                                                                                         <  1:247191/83‑1 (MQ=255)
                                                                             cGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAgg                                                                                        <  2:249301/100‑1 (MQ=255)
                                                                             cGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAgg                                                                                        >  1:249301/1‑100 (MQ=255)
                                                                                  tATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGcc                                            <  2:778191/139‑1 (MQ=255)
                                                                                     gCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAACGCGCCaa                                          >  2:1396214/1‑139 (MQ=255)
                                                                                          ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCa                                                                                                                         >  1:961404/1‑54 (MQ=255)
                                                                                          ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCa                                                                                                                         <  2:961404/54‑1 (MQ=255)
                                                                                          ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGAc                                                                             <  2:117246/98‑1 (MQ=255)
                                                                                          ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGAc                                                                             >  1:117246/1‑98 (MQ=255)
                                                                                            ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTggc                                                                                                                                 <  1:1383732/45‑2 (MQ=255)
                                                                                            ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTggc                                                                                                                                 >  2:1383732/1‑44 (MQ=255)
                                                                                            ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATg                                                                                                              <  2:420735/63‑1 (MQ=255)
                                                                                            ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATg                                                                                                              >  1:420735/1‑63 (MQ=255)
                                                                                                caACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGt                              <  1:1570697/139‑1 (MQ=255)
                                                                                                   cccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAAc                                                             >  1:1123613/1‑105 (MQ=255)
                                                                                                   cccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAAc                                                             <  2:1123613/105‑1 (MQ=255)
                                                                                                   cccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCtt                           <  1:1456228/139‑1 (MQ=255)
                                                                                                    ccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAACGCGCCAACTTAAGTGCAGTCtt                           <  1:1396214/139‑1 (MQ=255)
                                                                                                       gCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGttta                                                   <  1:827462/111‑4 (MQ=255)
                                                                                                       gCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTa                                                   >  2:827462/1‑111 (MQ=255)
                                                                                                       gCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGt                       >  2:1070623/1‑139 (MQ=255)
                                                                                                       gCGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGt                       >  1:534889/1‑139 (MQ=255)
                                                                                                        cGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGa                                                               <  2:885022/98‑1 (MQ=255)
                                                                                                        cGTTGTTATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGa                                                               >  1:885022/1‑98 (MQ=255)
                                                                                                             ttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGa                                                                              >  2:515800/1‑78 (MQ=255)
                                                                                                             ttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGa                                                                              <  1:515800/78‑1 (MQ=255)
                                                                                                             ttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGAc                                                                             >  2:118408/1‑79 (MQ=255)
                                                                                                             ttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGAc                                                                             <  1:118408/79‑1 (MQ=255)
                                                                                                             ttATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAg                 >  2:661857/1‑139 (MQ=255)
                                                                                                              tATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGc                >  2:1455331/1‑139 (MQ=255)
                                                                                                              tATTGCGGGTCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGc                <  1:1624491/139‑1 (MQ=255)
                                                                                                                       tCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCaaa       >  2:1728685/1‑139 (MQ=255)
                                                                                                                       tCGCGCGATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCaaa       >  1:592164/1‑139 (MQ=255)
                                                                                                                        cgcgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGAc                                                                                                        <  2:317124/41‑1 (MQ=255)
                                                                                                                        cgcgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGAc                                                                                                        >  1:317124/1‑41 (MQ=255)
                                                                                                                          cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGcc                                            >  2:107345/1‑99 (MQ=255)
                                                                                                                          cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGcc                                            <  1:107345/99‑1 (MQ=255)
                                                                                                                          cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCa                                >  1:1428471/1‑111 (MQ=255)
                                                                                                                          cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCa                                <  2:1428471/111‑1 (MQ=255)
                                                                                                                          cgcgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGaa    <  2:1266493/139‑1 (MQ=255)
                                                                                                                            cgATGCAGTT‑GGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATg  <  1:1455331/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                 
GGCGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATG  >  NC_000913/1006279‑1006542

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