Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,462,6770AT21.1% 79.5 / 15.2 38intergenic (+31/‑335)fadL/yfdFlong‑chain fatty acid outer membrane transporter/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (15/15);  new base T (4/4);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATA  >  NC_000913/2462530‑2462820
                                                                                                                                                   |                                                                                                                                               
gACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCa                                                                                                                                                 >  3:35528/1‑148 (MQ=255)
    ttGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCa                                                                                                                                            <  4:52732/149‑1 (MQ=255)
        tgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGt                                                                                                                                         >  1:177163/1‑148 (MQ=255)
          tttCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGa                                                                                                                                       <  2:42691/148‑1 (MQ=255)
             cTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACtt                                                                                                                                    <  8:75988/148‑1 (MQ=255)
                 tatGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGACTTTatc                                                                                                                                <  2:177163/148‑3 (MQ=255)
                                      aaTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAagtag                                                                                                          >  5:110280/1‑149 (MQ=255)
                                       aTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAagtagt                                                                                                         <  4:35528/149‑1 (MQ=255)
                                                     aTACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTa                                                                                           >  5:9265/1‑149 (MQ=255)
                                                        ccAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTAcct                                                                                        >  2:168300/1‑149 (MQ=255)
                                                           gTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTAcctcct                                                                                     >  3:104403/1‑149 (MQ=255)
                                                                  tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGtttt                                                                              <  6:110280/149‑1 (MQ=255)
                                                                    tGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGACTTTat                                                                                                                                 <  2:7431/96‑2 (MQ=38)
                                                                    tGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGACTTTat                                                                                                                                 >  1:7431/1‑95 (MQ=38)
                                                                     gAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGc                                                                            >  3:59868/1‑148 (MQ=255)
                                                                           aaaGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGc                                                                     >  1:180471/1‑149 (MQ=255)
                                                                              gCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTg                                                                  >  6:122919/1‑149 (MQ=255)
                                                                                     tGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAg                                                           <  2:180471/149‑1 (MQ=255)
                                                                                       ttCGGTACTAACTTTAACCACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGcc                                                                    >  4:59300/1‑138 (MQ=255)
                                                                                       ttCGGTACTAACTTTAACCACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGcc                                                                    <  3:59300/138‑1 (MQ=255)
                                                                                              cTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCtt                                                  >  8:4235/1‑149 (MQ=255)
                                                                                                 aCTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCttt                                                 <  5:122919/147‑1 (MQ=255)
                                                                                                 aCTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCttt                                                 <  6:9265/147‑1 (MQ=255)
                                                                                                 aCTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTc                                                >  7:67329/1‑148 (MQ=255)
                                                                                                         tACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCAACGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAAt                                        >  7:115821/1‑148 (MQ=255)
                                                                                                             cgTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCg                                   <  7:4235/149‑1 (MQ=255)
                                                                                                                  tGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCAtt                                                                                                >  4:133799/1‑83 (MQ=255)
                                                                                                                  tGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCAtt                                                                                                <  3:133799/83‑1 (MQ=255)
                                                                                                                    aTAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATc                              <  1:168300/147‑1 (MQ=255)
                                                                                                                    aTAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCca                            >  3:105560/1‑149 (MQ=255)
                                                                                                                                    gATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGaa              <  4:104403/147‑1 (MQ=255)
                                                                                                                                      tAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAggg           >  5:30972/1‑148 (MQ=255)
                                                                                                                                      tAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCACACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGAATCCACAGACAAACCAAAAggg           >  8:102327/1‑148 (MQ=255)
                                                                                                                                        aaGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaaa        <  4:105560/149‑1 (MQ=255)
                                                                                                                                          gTCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCAc                                                                                                                   <  7:2452/40‑1 (MQ=17)
                                                                                                                                          gTCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCAc                                                                                                                   >  8:2452/1‑40 (MQ=17)
                                                                                                                                          gTCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCcaca                          <  4:129851/129‑1 (MQ=255)
                                                                                                                                          gTCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCcaca                          >  3:129851/1‑129 (MQ=255)
                                                                                                                                              ccTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCata  <  6:30972/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATA  >  NC_000913/2462530‑2462820

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