Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,322,662 (A)7→8 intergenic (‑157/+71) secG ← / ← glmM Sec translocon subunit SecG/phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,322,6551.A96.5% 93.8 / ‑3.2 29intergenic (‑150/+78)secG/glmMSec translocon subunit SecG/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (12/16);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

TTCCTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCT  >  NC_000913/3322516‑3322787
                                                                                                                                            |                                                                                                                                     
ttCCTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCgg                                                                                                                                        >  2:595846/1‑140 (MQ=255)
ttCCTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCgg                                                                                                                                        >  1:563781/1‑140 (MQ=255)
  ccTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGGGAATGACCGTCGAAAGAGCAATTTTGTCGCGTTGTATCAACTGCGGaa                                                                                                                                     >  2:173033/1‑140 (MQ=255)
       cGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGaaaaaaa                                                                                                                               <  2:388008/140‑1 (MQ=255)
            ttGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGaaaaaaaacagc                                                                                                                          <  2:245450/140‑6 (MQ=255)
             tGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCa                                                                                                                          <  1:124035/140‑1 (MQ=255)
                  cgTACCAATGCTTCAACCAATACAGCGGGAGCTTTCAATCTTTCCCACAGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAAcc                                                                                                                     <  2:179538/140‑1 (MQ=255)
                  cgTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAAcc                                                                                                                     <  2:343732/140‑1 (MQ=255)
                     aCCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGAt                                                                                                                  <  2:71253/140‑1 (MQ=255)
                                     aTAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGaaaa                                                                                                  >  1:330196/1‑140 (MQ=255)
                                                    ttAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgcc                                                                                   <  1:268754/140‑1 (MQ=255)
                                                     tAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgccg                                                                                  <  1:518690/140‑1 (MQ=255)
                                                     tAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATcgccg                                                                                  <  2:563781/140‑1 (MQ=255)
                                                             cccACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTa                                                                          <  2:550248/140‑1 (MQ=255)
                                                               cACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGc                                                                        <  2:400452/140‑1 (MQ=255)
                                                                         ggTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTa                                                              <  1:376245/140‑1 (MQ=255)
                                                                                    tAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAAc                                                   >  1:506012/1‑140 (MQ=255)
                                                                                     aaCCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAAc                                                            >  2:7542/1‑130 (MQ=255)
                                                                                     aaCCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAAc                                                            <  1:7542/130‑1 (MQ=255)
                                                                                      aCCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACgg                                                 >  1:591147/1‑140 (MQ=255)
                                                                                       ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGc                                                <  1:69017/140‑1 (MQ=255)
                                                                                       ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGc                                                <  2:330196/140‑1 (MQ=255)
                                                                                       ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGc                                                <  1:155900/140‑1 (MQ=255)
                                                                                       ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGc                                                >  1:118413/1‑140 (MQ=255)
                                                                                            gCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTa                                           >  1:31993/1‑140 (MQ=255)
                                                                                                aTGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGc                                       >  2:47677/1‑140 (MQ=255)
                                                                                                        cGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGa                               >  1:588151/1‑140 (MQ=255)
                                                                                                         gCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTAATCTAACCCTTTAAACGGCTTTTACTGCATCGGCGAt                              >  1:274981/1‑140 (MQ=255)
                                                                                                            aGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCg                           <  1:595846/140‑1 (MQ=255)
                                                                                                             gggCAATTTTATCGCCTTGTATGAATTGCCGAAAAAACACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCgg                          <  1:173033/140‑1 (MQ=255)
                                                                                                                         cGCATTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAATGGCTTTTGCTGTATGGGCGATGCGGTGTGCAAATTCa              >  1:201020/1‑140 (MQ=255)
                                                                                                                            aTTGTATCAACTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACTCTTTAAACAGCTGTTAATGCATCGGGGATGCTGTGTGCAAATTCAGTc           >  2:295339/1‑140 (MQ=255)
                                                                                                                                     aCTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCt  <  2:338242/140‑1 (MQ=255)
                                                                                                                                     aCTGCGGAAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCAGTATGCAAAGTCAGTCACCTGCGCt  >  1:309225/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                     
TTCCTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCT  >  NC_000913/3322516‑3322787

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