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)
  cGATGCGGCCTGCGTTTGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACg                                                                                                                            >  2:1624491/1‑139 (MQ=255)
                tttGGTTGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATg                                                                                                              <  1:732977/139‑1 (MQ=255)
                     ttGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGa                                                                                                         <  1:946365/139‑1 (MQ=255)
                     ttGAGCGATACTCACGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGa                                                                                                         <  2:1430792/139‑1 (MQ=255)
                                   cGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGc                                                                                                                           >  1:1463806/1‑107 (MQ=255)
                                   cGAGCAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGc                                                                                                                           <  2:1463806/107‑1 (MQ=255)
                                       cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCt                                                                                                   >  1:1114885/1‑127 (MQ=255)
                                       cAGGTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCt                                                                                                   <  2:1114885/127‑1 (MQ=255)
                                          gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCAt                                                                                                               >  1:135387/1‑112 (MQ=255)
                                          gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCAt                                                                                                               <  2:135387/112‑1 (MQ=255)
                                          gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTg                                                                                                  <  2:419346/125‑1 (MQ=255)
                                          gTCAATTTGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTg                                                                                                  >  1:419346/1‑125 (MQ=255)
                                                ttGCTGGTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGa                                                                              <  1:846466/139‑1 (MQ=255)
                                                     ggTCGTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCt                                                                         <  2:1214326/139‑1 (MQ=255)
                                                         gTTGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAg                                                                     >  2:1570697/1‑139 (MQ=255)
                                                          ttGAATCTGGCGAAAATGCCGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGc                                                                    <  2:934705/139‑1 (MQ=255)
                                                                             cGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGa                                                                                                         >  2:247191/1‑83 (MQ=255)
                                                                             cGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGa                                                                                                         <  1:247191/83‑1 (MQ=255)
                                                                             cGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAgg                                                                                        <  2:249301/100‑1 (MQ=255)
                                                                             cGCACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAgg                                                                                        >  1:249301/1‑100 (MQ=255)
                                                                                  tATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGcc                                            <  2:778191/139‑1 (MQ=255)
                                                                                     gCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAACGCGCCaa                                          >  2:1396214/1‑139 (MQ=255)
                                                                                          ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCa                                                                                                                         >  1:961404/1‑54 (MQ=255)
                                                                                          ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCa                                                                                                                         <  2:961404/54‑1 (MQ=255)
                                                                                          ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGAc                                                                             <  2:117246/98‑1 (MQ=255)
                                                                                          ctgGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGAc                                                                             >  1:117246/1‑98 (MQ=255)
                                                                                            ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTggc                                                                                                                                 <  1:1383732/45‑2 (MQ=255)
                                                                                            ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTggc                                                                                                                                 >  2:1383732/1‑44 (MQ=255)
                                                                                            ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATg                                                                                                              <  2:420735/63‑1 (MQ=255)
                                                                                            ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATg                                                                                                              >  1:420735/1‑63 (MQ=255)
                                                                                                caACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGt                              <  1:1570697/139‑1 (MQ=255)
                                                                                                   cccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAAc                                                             >  1:1123613/1‑105 (MQ=255)
                                                                                                   cccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAAc                                                             <  2:1123613/105‑1 (MQ=255)
                                                                                                   cccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCtt                           <  1:1456228/139‑1 (MQ=255)
                                                                                                    ccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAACGCGCCAACTTAAGTGCAGTCtt                           <  1:1396214/139‑1 (MQ=255)
                                                                                                       gCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGttta                                                   <  1:827462/111‑4 (MQ=255)
                                                                                                       gCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTa                                                   >  2:827462/1‑111 (MQ=255)
                                                                                                       gCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGt                       >  2:1070623/1‑139 (MQ=255)
                                                                                                       gCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGt                       >  1:534889/1‑139 (MQ=255)
                                                                                                        cGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGa                                                               <  2:885022/98‑1 (MQ=255)
                                                                                                        cGTTGTTATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGa                                                               >  1:885022/1‑98 (MQ=255)
                                                                                                             ttATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGa                                                                              >  2:515800/1‑78 (MQ=255)
                                                                                                             ttATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGa                                                                              <  1:515800/78‑1 (MQ=255)
                                                                                                             ttATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGAc                                                                             >  2:118408/1‑79 (MQ=255)
                                                                                                             ttATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGAc                                                                             <  1:118408/79‑1 (MQ=255)
                                                                                                             ttATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAg                 >  2:661857/1‑139 (MQ=255)
                                                                                                              tATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGc                >  2:1455331/1‑139 (MQ=255)
                                                                                                              tATTGCGGGTCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGc                <  1:1624491/139‑1 (MQ=255)
                                                                                                                       tCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCaaa       >  2:1728685/1‑139 (MQ=255)
                                                                                                                       tCGCGCGATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCaaa       >  1:592164/1‑139 (MQ=255)
                                                                                                                        cgcgcgATGCAGTTGGCGACGCCATTAAAATCATGAACGAc                                                                                                        <  2:317124/41‑1 (MQ=255)
                                                                                                                        cgcgcgATGCAGTTGGCGACGCCATTAAAATCATGAACGAc                                                                                                        >  1:317124/1‑41 (MQ=255)
                                                                                                                          cgcgATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGcc                                            >  2:107345/1‑99 (MQ=255)
                                                                                                                          cgcgATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGcc                                            <  1:107345/99‑1 (MQ=255)
                                                                                                                          cgcgATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCa                                >  1:1428471/1‑111 (MQ=255)
                                                                                                                          cgcgATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCa                                <  2:1428471/111‑1 (MQ=255)
                                                                                                                          cgcgATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGaa    <  2:1266493/139‑1 (MQ=255)
                                                                                                                            cgATGCAGTTGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTTTAACGCGCCAACTTAAGTGCAGTCTTCGGTACACAGCTACAGCAAAGAATg  <  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: