Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913124,7440T.23.3% 112.2 / 31.4 43coding (1728/2664 nt)aceEpyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base T (16/17);  new base . (5/5);  total (21/22)
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.16e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGTCTGTTCCGTCAGATTGGTATTTACAGCCCGAACGGTCAGCAGTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATC  >  NC_000913/124607‑124891
                                                                                                                                         |                                                                                                                                                   
ggTCTGTTCCGTCAGATTGGTATTTACAGCCCGAACGGTCAGCAGTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCctggct                                                                                                                                          <  8:3370/149‑1 (MQ=255)
                 tGGTATTTACAGCCCGAACGGTCAGCAGTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGG‑TGTTCCTGGCTGGCAGCGGCGACCTCTTa                                                                                                                        >  2:41903/1‑149 (MQ=255)
                     aTTTACAGCCCGAACGGTCAGCAGTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAg                                                                                                                     >  1:96511/1‑149 (MQ=255)
                     aTTTACAGCCCGAACGGTCAGCAGTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAg                                                                                                                     >  3:116963/1‑149 (MQ=255)
                           aGCCCGAACGGTCAGCAGTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGG‑TGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACcaac                                                                                                              <  3:75035/149‑1 (MQ=255)
                              ccGAACGGTCAGCAGTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACcaacaa                                                                                                            >  5:34412/1‑149 (MQ=255)
                              ccGAACGGTCAGCAGTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACcaacaa                                                                                                            <  4:116963/149‑1 (MQ=255)
                                gAACGGTCAGCAGTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGAGCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATc                                                                                                          <  4:76792/149‑1 (MQ=255)
                                   cGGTCAGCAGTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGGCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGGTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGc                                                                                                       <  1:139935/149‑1 (MQ=255)
                                     gTCAGCAGTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGG‑TGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCga                                                                                                    <  1:41903/149‑1 (MQ=255)
                                      tCAGCAGTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGG‑TGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCgat                                                                                                   >  4:35394/1‑149 (MQ=255)
                                         gcagTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTTTAAAGAAGACGAGAAAGGTCAGATTCTGCAGTAAGGAATCAACGAGAAGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGCCCTCTTACAGCACCAACAATCTGCCcatg                                                                                                  >  2:94424/1‑148 (MQ=255)
                                                    cGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTAc                                                                                      <  2:96511/149‑1 (MQ=255)
                                                     gCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTAc                                                                                      >  4:49114/1‑148 (MQ=255)
                                                     gCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTAc                                                                                      <  3:49114/148‑1 (MQ=255)
                                                                    ggTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTAc                                                                             <  5:152174/142‑1 (MQ=255)
                                                                    ggTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTAc                                                                             >  6:152174/1‑142 (MQ=255)
                                                                     gTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCTACCTCTTACAGCACCAACAATCTTCCGATTATCCCGTTCTACATCTATTACTCGATGtt                                                                     <  1:94424/149‑1 (MQ=255)
                                                                        gCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCg                                                                   <  6:34412/148‑1 (MQ=255)
                                                                        gCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCGTACAGCACCAACACTCTGCCGAGGCTCCCGTTCTACATCTATTACTCGATGTTCgg                                                                  >  8:7928/1‑149 (MQ=255)
                                                                         cTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACcaa                                                                                                               <  4:135151/103‑1 (MQ=255)
                                                                         cTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACcaa                                                                                                               >  3:135151/1‑103 (MQ=255)
                                                                            aCTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTc                                                              >  6:26594/1‑149 (MQ=255)
                                                                             cTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACa                                                                                     <  4:36319/125‑1 (MQ=255)
                                                                             cTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACa                                                                                     >  3:36319/1‑125 (MQ=255)
                                                                              tataAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTc                                                              >  7:98829/1‑147 (MQ=255)
                                                                               ataAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGG‑TGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCa                                                            >  2:8343/1‑147 (MQ=255)
                                                                               ataAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGG‑TGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCa                                                            <  1:8343/147‑1 (MQ=255)
                                                                                   agaagaCGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTa                                                       >  5:22602/1‑149 (MQ=255)
                                                                                            gaAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGG‑TGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCt                                             >  8:47558/1‑149 (MQ=255)
                                                                                                 gTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGc                                         <  6:22602/149‑1 (MQ=255)
                                                                                                    aGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTgg                                      >  1:132931/1‑149 (MQ=255)
                                                                                                            cAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTgg                              <  5:26594/149‑1 (MQ=255)
                                                                                                             aGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGc                             <  8:98829/149‑1 (MQ=255)
                                                                                                              ggAAGGGAGCAACGAGCTGGGCGCAGGTGGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCg                            <  7:7928/149‑1 (MQ=255)
                                                                                                               gAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGa                           >  6:87834/1‑149 (MQ=255)
                                                                                                               gAAGGGATCAACGAGCTGGGCGCAGG‑TGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGAc                          >  3:141734/1‑149 (MQ=255)
                                                                                                                aaGGGATCAACGAGCTGGGCGCAGG‑TGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGAcc                         <  4:141734/149‑1 (MQ=255)
                                                                                                                 aGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGAc                          >  6:99621/1‑148 (MQ=255)
                                                                                                                  gggATCAACGAGCTGGGCGCAGG‑TGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGTCGCCAGCCGCCCAg                       <  3:35394/149‑1 (MQ=255)
                                                                                                                         aCGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAgcg                 >  7:97460/1‑149 (MQ=255)
                                                                                                                             gCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTg             <  2:132931/149‑1 (MQ=255)
                                                                                                                                       ggTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGAt   <  5:99621/149‑1 (MQ=255)
                                                                                                                                        gttgttCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATc  >  2:37489/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
GGTCTGTTCCGTCAGATTGGTATTTACAGCCCGAACGGTCAGCAGTACACCCCGCAGGACCGCGAGCAGGTTGCTTACTATAAAGAAGACGAGAAAGGTCAGATTCTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATC  >  NC_000913/124607‑124891

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: ‑