Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,118,621 (A)7→8 coding (107/246 nt) zapB → cell division factor ZapB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,118,6141.A95.7% 156.0 / ‑0.7 48coding (100/246 nt)zapBcell division factor ZapB
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (27/18);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGG  >  NC_000913/4118485‑4118731
                                                                                                                                  |                                                                                                                       
tCCATCAACGACTTGCAATTCAGGCGAGGTATAACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTgaaagaaa                                                                                                                      <  2:222989/133‑1 (MQ=255)
   aTCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGCTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGaaaaaa                                                                                                                  >  1:411508/1‑133 (MQ=255)
     cAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAGGAGCTGAAAGaaaaaaaa                                                                                                                >  1:178504/1‑132 (MQ=255)
     cAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTgaaagaaa                                                                                                                      <  1:233384/128‑1 (MQ=255)
     cAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTgaaagaaa                                                                                                                      >  2:233384/1‑128 (MQ=255)
            ttGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCa                                                                                                         <  1:37577/133‑1 (MQ=255)
            ttGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCa                                                                                                         <  1:110911/133‑1 (MQ=255)
                          agGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCAAAGAGCTGAAAGAAAATAAACAACTGACTGTTGCAGGAAGt                                                                                           >  2:298415/1‑133 (MQ=255)
                            gTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTc                                                                                         >  1:155816/1‑133 (MQ=255)
                                gACAATGTAATTAGAAGTATTTGAGAAACTTGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCaaaa                                                                                     <  1:228793/133‑1 (MQ=255)
                                 aCAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAAt                                                                                    >  2:386857/1‑133 (MQ=255)
                                 aCAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACCCAATCACTGTCTCAGGCCGTTCaaagt                                                                                    >  2:291580/1‑131 (MQ=255)
                                    aTGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGcc                                                                                 <  1:386857/133‑1 (MQ=255)
                                    aTGTCATTAGAAGTGCTTGAGAAATTGGCAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGcc                                                                                 <  2:178504/133‑1 (MQ=255)
                                     tGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGccc                                                                                >  1:454076/1‑133 (MQ=255)
                                         aTTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCa                                                                            >  1:425047/1‑133 (MQ=255)
                                         aTTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCa                                                                            <  1:335201/133‑1 (MQ=255)
                                            aGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGAAAGATCCAAATGCCCAGCATCa                                                                         >  1:284309/1‑133 (MQ=255)
                                               aGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAgcg                                                                      >  2:130418/1‑133 (MQ=255)
                                                 tgtTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGcgcg                                                                    <  2:397667/133‑1 (MQ=255)
                                                 tgtTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGcgcg                                                                    >  1:604215/1‑133 (MQ=255)
                                                   tttGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCTGGAAGTTCAAAATGCCCAGCATCAGcgcg                                                                    >  1:78172/1‑131 (MQ=255)
                                                   tttGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCTGGAAGTTCAAAATGCCCAGCATCAGcgcg                                                                    <  2:78172/131‑1 (MQ=255)
                                                   tttGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGaa                                                                  <  2:368935/133‑1 (MQ=255)
                                                     tGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCCCTGTCGCAGGAAGTTCAAAATGCCCGGCATCAGCGCGAaga                                                                >  2:227740/1‑133 (MQ=255)
                                                     tGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAaga                                                                >  1:101694/1‑133 (MQ=255)
                                                       agaAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGc                                                              <  2:373013/133‑1 (MQ=255)
                                                             tGGAAGCAGAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGc                                                        >  2:471537/1‑133 (MQ=255)
                                                             tGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGc                                                        <  1:606451/133‑1 (MQ=255)
                                                             tGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGc                                                        <  1:563619/133‑1 (MQ=255)
                                                                   cAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTgaga                                                  >  2:447485/1‑133 (MQ=255)
                                                                      aaGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGaaca                                               >  1:443548/1‑133 (MQ=255)
                                                                        gTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGaacaac                                             <  1:471537/133‑1 (MQ=255)
                                                                        gTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGaacaac                                             <  1:447485/133‑1 (MQ=255)
                                                                               agGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGa                                      <  1:8536/133‑1 (MQ=255)
                                                                                ggCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCATCATCAGCGGGAAAAGCTGGAGCGTGAGAACAACCATCTGaa                                     >  2:210640/1‑133 (MQ=255)
                                                                                ggCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGaa                                     >  2:114128/1‑133 (MQ=255)
                                                                                          aCCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAAcagcag                           >  2:479303/1‑133 (MQ=255)
                                                                                             aTCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAAc                        >  1:367402/1‑133 (MQ=255)
                                                                                             aTCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAAc                        >  1:122851/1‑133 (MQ=255)
                                                                                              tCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACg                       >  2:89834/1‑133 (MQ=255)
                                                                                                   ctGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTgg                  >  2:316320/1‑133 (MQ=255)
                                                                                                    tGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGc                 <  1:89834/133‑1 (MQ=255)
                                                                                                      ttGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAAGCATCTGAAAGAACATCAGAACGGCTGGCAg               >  1:213312/1‑133 (MQ=255)
                                                                                                      ttGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAg               >  2:323296/1‑133 (MQ=255)
                                                                                                      ttGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAg               >  1:371309/1‑133 (MQ=255)
                                                                                                       tGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGa             >  1:11157/1‑133 (MQ=255)
                                                                                                              ggAAATCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCt       >  2:460759/1‑133 (MQ=255)
                                                                                                                aaaTCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAACAACAGCACAACGGCAGGCAGGAACGTCTGc     >  1:257519/1‑133 (MQ=255)
                                                                                                                aaaTCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGATGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGc     >  2:601518/1‑133 (MQ=255)
                                                                                                                   tCGAAGAGCTGAAAGTAAAAAAACCACTAACTATCGCAGGACGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAgg  <  2:213312/133‑1 (MQ=255)
                                                                                                                   tCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAgg  <  2:363688/133‑1 (MQ=255)
                                                                                                                   tCGAAGAGCTGAAAGAAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAgg  <  1:130418/133‑1 (MQ=255)
                                                                                                                                  |                                                                                                                       
TCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGAGCGTGAGAACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGG  >  NC_000913/4118485‑4118731

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: