Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,462,8891.C56.0% 42.5 / 35.6 25intergenic (+243/‑123)fadL/yfdFlong‑chain fatty acid outer membrane channel/bacteriophage T2 receptor/protein YfdF
Reads supporting (aligned to +/- strand):  ref base . (7/4);  new base C (7/7);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCC‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAA  >  NC_000913/2462646‑2463115
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
aTAACGCGTGCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCCGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGGGTTTCCCTTTCTGCATTTCCAGGAT‑ACTcc‑ccccccccctggct                                                                                                                                                                                                                              <  1:39784/251‑8 (MQ=255)
       gTTCGCCTGGATAAAGTTACCTGCATAGCAGGTGACTTTAATTCCACCACTTCACCGAAAGTAGTGTACTCATTGCTTACCTCCTGAGTTTTGAACACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTACCCACAGACAAACCTGAAGGGTAAGGCATATTTTCAGGCGTTCTGCTCGTCGTTCTCAAGGAGTTTACTTTTCTGAATTTCCAGGATCAACCC‑CCCCCCC‑TGGCTATtgtggg                                                                                                                                                                                                                       >  2:283610/1‑249 (MQ=255)
           gCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCCCTGGCTATTGTGCGCTc                                                                                                                                                                                                                    >  1:50233/1‑250 (MQ=255)
             cTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACATCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCC‑TGGCTATTGTGCGCTCATac                                                                                                                                                                                                                >  1:280946/1‑251 (MQ=255)
                   aaagtcacctcaatggctgttGACTTTAACTCCCCCACTGCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTGTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGGCCTTCTCAAAGATTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCC‑TGGCTATTGTGCGCTCATACACTCaa                                                                                                                                                                                                          <  2:372315/231‑1 (MQ=255)
                         aCCTGCATAGCAGGTGAATTTCACTCCCACACTTCACCGAACGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGATTACTTTTCTGCATTTCCAGGAT‑ACTCCACCCCCCC‑TGGCTATTGTGAGCTCCTTCACCCAAAttc                                                                                                                                                                                                      >  2:223977/1‑249 (MQ=255)
                               aTAGCAGGTGACTTTAACTCCCCCACTTCACCGTAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGCCAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTTTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCCCTGGATCTTGTGTGCGCCTACACTCAACTTTAAGactg                                                                                                                                                                                               >  2:260165/1‑248 (MQ=255)
                                           tttAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACATGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCtaaataa                                                                                                                                                                                    >  1:279824/1‑250 (MQ=255)
                                                       ccgtCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCACATCAATAAGCGGTATCAACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCC‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGAGTTTTTTg                                                                                                                                                                      >  2:93998/4‑251 (MQ=255)
                                                           cgccgaaataagtGCCCTCATTGCTTACCGCCTAAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCGTACAAATCAATAAGCGGTAGCCACAGACAAACCGGAAGGGAAAGGCATATTTTCAGGCTTTCTGCTCGTCCGTCTCAAAGAGTTTATTTTTCTGCATTTCCAAGAT‑ACTCC‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAAGTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGAt                                                                                                                                                                    <  1:260165/240‑1 (MQ=255)
                                                                                                   gCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAACGGTATCCACAAACAAACCTGAAGGGAAAGGCATATTGTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCCCTGGCTATTGTTCGCTCATACACTCAAATTAAAGATAGGTTTTAAATAAAGGAGCGTTTTTTGATAGTCTATTTGATTAGGTAATATATATTTGTAACAaatcaa                                                                                                                            <  1:53169/250‑1 (MQ=255)
                                                                                                          gCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACACACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATg                                                                                                                    >  1:307160/1‑251 (MQ=255)
                                                                                                          gCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATg                                                                                                                    >  1:146095/1‑251 (MQ=255)
                                                                                                                            cTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑CCTCC‑CCCCCCC‑TGGCTATTGTGCGCTCATACACTCAACTTACAGTTAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCCTTAGGTAATATATATTTGTAACAAATCAATCCAAATGGCCTAAAATCCTGCTAcc                                                                                                  >  1:20185/1‑250 (MQ=255)
                                                                                                                               tCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCAt                                                                                                <  1:353954/250‑1 (MQ=255)
                                                                                                                                   aTCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCC‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTAttt                                                                                          <  1:340867/251‑1 (MQ=255)
                                                                                                                                             ggTAGCCACAGCAAAACCTGAAAGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑ACCCCCC‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTGTTGATAGTCTATTTCATAAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATcaacaa                                                                                <  2:181654/251‑1 (MQ=255)
                                                                                                                                                            aCCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCtt                                                                  >  1:234426/1‑251 (MQ=255)
                                                                                                                                                                 caGGGAAAGGCATGTATACCGGCGTTCTGCTTCTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGTT‑ACTCC‑CCCCCCC‑TGGCTATTGTGCCATCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTTTTTCATTAGGTATTATATATTTGTACCAAATCAATCAAAATGGAATAAACTCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTAcca                                                              >  2:254659/2‑248 (MQ=255)
                                                                                                                                                                             tatTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTAC‑CCCCCCC‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAAAAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATg                                                >  1:56675/1‑251 (MQ=255)
                                                                                                                                                                                      gCGTTCTGCTCGTCCTTCTCAAAGAGTTTTCTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAAc                                        <  2:146095/251‑1 (MQ=255)
                                                                                                                                                                                       cGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACCCC‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACa                                       <  1:363902/251‑1 (MQ=255)
                                                                                                                                                                                                                        tCTGCATTTCCAGGAT‑ACTCC‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTc                                                                                                                                                          >  1:4568/1‑103 (MQ=255)
                                                                                                                                                                                                                        tCTGCATTTCCAGGAT‑ACTCC‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTc                                                                                                                                                          <  2:4568/103‑1 (MQ=255)
                                                                                                                                                                                                                        tCTGCATTTCCAGGAT‑ACTCC‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCTTCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTc       <  1:223977/250‑1 (MQ=255)
                                                                                                                                                                                                                           gcattgacAGGAT‑ACGCC‑CCCCCCC‑AGGCTATTGTGCGCTCATACACTCAAAGTAAAGATAGGGTCTAAATAAAGGAGCGTTTTTTGATAGTCTATGTCATTAGGTAATATATATTTGTAACAAATCAAGCAAAATGGAAGAAAATCAGGCTAACATCTATTTCAATCAACAAAACCAGCGCAGCTGACCCAGAATCCAGCAATGAAAATAACAAGGATGAAGTTAATGGATGAGTACAAGAGTTCaaaaa  <  2:20185/244‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
ATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGAT‑ACTCC‑CCCCCCC‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCATCTATTTCAATCAACAATACCAGCGCAGCTTACCCAGAATCCATCAATGAAAATAACAATGATGAAGTTAATGGATTAGTACAAGAGTTCAAAAA  >  NC_000913/2462646‑2463115

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