Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi3,546,3910CT53.8% 6.9 / 45.7 39intergenic (+192/‑168)nfuA/gntTFe/S biogenesis protein, putative scaffold/chaperone protein/gluconate transporter, high‑affinity GNT I system
Reads supporting (aligned to +/- strand):  ref base C (9/9);  new base T (12/9);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAGTCGTGGTGCCGCAATATCTCTGTTCCCGATTGGGATAATTAGAGTTTGTCGTCAGAAAATTGACGTTACCCATAA  >  NC_000913_3_bme_pgi/3546244‑3546536
                                                                                                                                                   |                                                                                                                                                 
cATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCCGt                                                                                                                                                 <  1:339438/150‑1 (MQ=255)
  tttCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCTGTAt                                                                                                                                               >  4:100714/1‑150 (MQ=255)
                  aaCATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCTGTATCATTCACATGGGTAcc                                                                                                                               >  3:222199/1‑150 (MQ=255)
                         gCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATa                                                                                                                        >  3:125248/1‑150 (MQ=255)
                            aGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTc                                                                                                                     <  1:287015/150‑1 (MQ=255)
                                       cagcaTTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGa                                                                                                          <  4:222199/150‑1 (MQ=255)
                                           aTTGCTTTGTCAAATCCTCGGTTCATTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTAc                                                                                                      <  3:100714/150‑1 (MQ=255)
                                                tttGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACaataa                                                                                                 >  2:290263/1‑150 (MQ=255)
                                                    tCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAg                                                                                             <  4:125248/150‑1 (MQ=255)
                                                      aaaTCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGtt                                                                                           >  1:63936/1‑150 (MQ=255)
                                                         tCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTa                                                                                        >  3:177173/1‑150 (MQ=255)
                                                            tCGGTTCGTTTTCCCGCGATAATATGACCAAATACTCATAATTTAAATTTAACCCGCTCTGGTGATTCTCAAACGCCAGATGTTACACGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTg                                                                                     >  2:476475/1‑150 (MQ=255)
                                                              ggTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTggg                                                                                   >  4:172584/1‑150 (MQ=255)
                                                                    tttCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAg                                                                             >  2:426380/1‑150 (MQ=255)
                                                                                       ccgaCCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCCGTATCATTca                                                                                                                                         >  2:132748/4‑71 (MQ=255)
                                                                                       ccgaCCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCCGTATCATTca                                                                                                                                         <  1:132748/68‑1 (MQ=255)
                                                                                               tcATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAGTCg                                                                          <  2:217920/126‑1 (MQ=255)
                                                                                               tcATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAGTCg                                                                          >  1:217920/1‑126 (MQ=255)
                                                                                               tcATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAGTCGTGGTGCCGCAATATCTCTGTTccc                                                  >  4:208569/1‑150 (MQ=255)
                                                                                                   aatttaaatttaCCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAGTCGTGGTGCCGCAATATCTCTGTTCCCGAtt                                              <  4:177173/150‑1 (MQ=255)
                                                                                                      ttaaatttaCCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCTGTATCATTCACATGGGTACCAAAc                                                                                                                           >  2:41171/1‑70 (MQ=255)
                                                                                                      ttaaatttaCCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCTGTATCATTCACATGGGTACCAAAc                                                                                                                           <  1:41171/70‑1 (MQ=255)
                                                                                                                       tGGTGATTCTCAAACGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAGTCGTGGTGCCGCAATATCTCTGTTCCCGATTGGGATAATTAGAGTTTgtcg                          >  3:278484/1‑150 (MQ=255)
                                                                                                                         gTGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAg                                                                                             >  2:67388/1‑81 (MQ=255)
                                                                                                                         gTGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAg                                                                                             <  1:67388/81‑1 (MQ=255)
                                                                                                                          tGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCt                                                                                                                   <  3:203372/58‑1 (MQ=255)
                                                                                                                          tGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCt                                                                                                                   >  4:203372/1‑58 (MQ=255)
                                                                                                                          tGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAGTCGTGGTGCCGCAATATCTCTGTTCCCGATTGGGATAATTAGAGTTTGTCGTCa                       <  3:172584/150‑1 (MQ=255)
                                                                                                                           gATTCTCAAACGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAGTCGTGGTGCCGCAATAtct                                                          <  3:43784/114‑1 (MQ=255)
                                                                                                                           gATTCTCAAACGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTACCAGTCGTGGTGCCGCAATAtct                                                          >  4:43784/1‑114 (MQ=255)
                                                                                                                             ttCTCAAACGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAGTCGTGGTGCCGCAATATCTCTGTTCCCGATTGGGATAATTAGAGTTTGTCGTCAGaa                    <  1:426380/150‑1 (MQ=255)
                                                                                                                                 cAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAGTCGTGGTGCCGCAATATCTCTGTTCCCGATTGGGATAATTAGAGTTTGTCGTCAGAAAAtt                <  3:208569/150‑1 (MQ=255)
                                                                                                                                  aaaCGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGa                                                                                                                 >  3:252490/1‑52 (MQ=255)
                                                                                                                                  aaaCGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGa                                                                                                                 <  4:252490/52‑1 (MQ=255)
                                                                                                                                  aaaCGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAGTCGTGGTGCCGCAATATCTCTGTTCCCGATTGGGATAATTAGAGTTTGTCGTCAGAAAATTg               >  2:25349/1‑150 (MQ=255)
                                                                                                                                  aaaCGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAGTCGTGGTGCCGCAATATCTCTGTTCCCGATTGGGATAATTAGAGTTTGTCGTCAGAAAATTg               >  2:25213/1‑150 (MQ=255)
                                                                                                                                   aaCGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGa                                                                                                                 <  2:41604/51‑1 (MQ=255)
                                                                                                                                   aaCGCCAGATGTTACCTGTATCATTCACATGGGTACCAAACATACTCCTGa                                                                                                                 >  1:41604/1‑51 (MQ=255)
                                                                                                                                               tACCTGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAGTTATTAGTTTTAGTGGGTATCAGTCGTGGTGACACAATAGCTCTGTTCCCGATTGGGATAATTAGAGTTTGTCGTCAGAAAATTGACGTTACCCATaa  <  2:63936/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CATTTCGGGCCTACGTTGAACATGTGCCAGACGTTGGCGCAGCATTGCTTTGTCAAATCCTCGGTTCGTTTCCCCGCGATAATATGACCAACCTCTCATAATTTAAATTTACCCCGCTCTGGTGATTCTCAAACGCCAGATGTTACCCGTATCATTCACATGGGTACCAAACATACTCCTGACATCTGACTACAATAATTAGTTTTAGTGGGTATCAGTCGTGGTGCCGCAATATCTCTGTTCCCGATTGGGATAATTAGAGTTTGTCGTCAGAAAATTGACGTTACCCATAA  >  NC_000913_3_bme_pgi/3546244‑3546536

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑