Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143188,0120TA34.7% 40.6 / 25.1 49G957G (GGT→GGACJ019_00185hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/17);  new base A (0/17);  total (15/34)
Fisher's exact test for biased strand distribution p-value = 6.62e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCT  >  ACB143/187886‑188146
                                                                                                                              |                                                                                                                                       
cAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACACAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACCCCACAGGCGGTAATTTCGAg                                                                                                                               >  6:55535/1‑136 (MQ=255)
 aaCGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                                                                                                >  3:180205/1‑134 (MQ=37)
      ccAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACAGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTg                                                                                                                         >  1:268174/1‑136 (MQ=37)
                       aaGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGAcccc                                                                                                        >  1:291560/1‑136 (MQ=37)
                       aaGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACAGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAACCGAccc                                                                                                         >  3:122864/1‑135 (MQ=21)
                        aGCCGTCAACACCCAGGGCAGCGTCGAATTCCAGACCCCGGCTAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCg                                                                                                       <  1:257642/136‑1 (MQ=255)
                         gCCGGCAAAACCCAGGGCAGCGTCGAGTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACAGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCgg                                                                                                      >  8:265685/1‑136 (MQ=21)
                              cAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGGCTACAACAACGGTTCGGCTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGACCCCGCCCCCGGCTCAg                                                                                                 >  6:177374/1‑136 (MQ=255)
                               aaaaCCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAACCGACCCCGGCTCAGa                                                                                                >  8:41181/1‑136 (MQ=37)
                                           aGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGCGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACg                                                                                    <  7:141992/136‑1 (MQ=37)
                                           aGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCa                                                                                                                             <  3:77661/95‑1 (MQ=33)
                                           aGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAACCGACCCCGGCTCAGa                                                                                                >  3:96593/1‑124 (MQ=25)
                                                cGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                               >  7:64932/1‑136 (MQ=37)
                                                  aCTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCgg                                                                              >  4:180789/1‑135 (MQ=37)
                                                                          tCTACAACAACGTTTCGACGGTCAGCGGCACCATCCACAACGCAACCGGCGTAAATTTCGAGCACATGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGAc                                                     <  1:25633/136‑1 (MQ=11)
                                                                              caacaaCGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAAGCCACGACCGCGACCCTGACGGc                                                  >  2:251716/1‑135 (MQ=21)
                                                                                                  gcgcCACCATCCACACCCCCACTGGGGGAAATTTCCACCAACTCACCCCCAACCCGACCCCGGCCCAGACCACCATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAg                              <  3:190820/135‑1 (MQ=11)
                                                                                                  gCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAgg                             >  6:91927/1‑136 (MQ=32)
                                                                                                  gCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAGGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAg                              >  3:95699/1‑135 (MQ=17)
                                                                                                    gTCACCATCGACACGGCCCCTGGCGTTATTTTCGAGCACCTCACCCCGAACCCGACCCCGGCTCACCCCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                           <  6:250078/136‑1 (MQ=11)
                                                                                                     tCACCATCGCCAACCCCACGGGCGTAAATTACGACCAACTGACCCCCAACCCGACCCCGGCGCAGACCCCCATCACCGACTCGGTCGAAGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                           <  4:239959/135‑1 (MQ=11)
                                                                                                     tCACCATCGACAACGCCACAGGCGGTAATTTCGCGCAACTGACCCCGAACCCGACCCCGCCTCAGACCCCGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                           <  4:316521/135‑1 (MQ=255)
                                                                                                      cacaaccGAGACCGCCACCGGCGTAAATTTCGCCCAACTCCCCCCCAACCCGCCCCCGGCCCACACCACGACCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                         <  6:320706/130‑1 (MQ=11)
                                                                                                      cACCATCGAGACCGCCCCGGGGGGAAATTCCTACCACCGGACCCCGAACCCGCCCCCGGCTCAGCCCCCTATCAACGACTCGGTCGATCCCACAACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                         <  8:306764/136‑1 (MQ=11)
                                                                                                      cACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                         <  2:61302/136‑1 (MQ=17)
                                                                                                      cACCATCGACACCGCCACCGGGGGTAATTTCGACCACCTGCCCCCGAACCCGCCCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                         <  5:89492/136‑1 (MQ=11)
                                                                                                      cACCATCGACACCGCCACCGGCGGAAATTTCGACCACCTCACCCCTAACCCGCCCCCGGCCCACACCACGAACAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                         <  6:1058/136‑1 (MQ=11)
                                                                                                      cACCATCGACACCGCCACAGGCGTTAATTTCGAGCAACTGCCCCCCAACCCCACCCCGGCCCCGCCCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                         <  7:312871/136‑1 (MQ=11)
                                                                                                       caaaacGAGAACGCCACTGGCGGAAATTTCCACCACCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                         <  2:55149/130‑1 (MQ=11)
                                                                                                       accaccGACACCCCCACAGGGGGTAATTACGACCAACTCACCCCCAACCCGACCCCGGCTCCCACCACTATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                         <  3:49559/135‑1 (MQ=11)
                                                                                                       aacATCTAGAACCCCACCGGCGTAAATTTCGACCACATCACCCCGAACCCGCCCCGGGCTCACACCACTATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                         <  4:183869/133‑1 (MQ=11)
                                                                                                       aCCATCGAGAACGCCACCGGCGGAAATTTCCACCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGAACAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                         <  6:180348/135‑1 (MQ=11)
                                                                                                       aCCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                         <  7:269914/135‑1 (MQ=32)
                                                                                                       aCCATCGACAACGCCACCGGCGTAAATTACCACCAACTCACCCCGAACCCGACCCCGGCGCACACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                         <  4:373356/135‑1 (MQ=11)
                                                                                                       aCCATCCAGAACGCCACAGGCGGAAATTTCGAGCAACTCCCCCCGAACCCCCCCCCGGCTCAGACCACGATCCACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                         <  4:58273/135‑1 (MQ=11)
                                                                                                        ccATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                       <  8:292598/136‑1 (MQ=32)
                                                                                                        ccATCGACACCGCCACCGGCGGAAATTGCGACCACCTCCCCCCCAACCCGCCCCCGGCTCAGACCACGAGCAACGACTCGGGCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                       <  2:45129/136‑1 (MQ=11)
                                                                                                             gagaACGCCACCGGCGTAAATTTCCACCACCTCCCCCCGAACCCGACCCCGGCACACACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATc                   <  6:141823/135‑1 (MQ=11)
                                                                                                             gagaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATc                   <  7:27565/135‑1 (MQ=32)
                                                                                                              agaACGCCACCGCGGTAAATTGGGAGCAACTGACCCCGAACCCGACCGCGGCTCAGACCACGATCAACGACTGGGTCGATGCCACCACCGCTACCCTTACGGCGAGCCCGTCGGTCACGGAAGGGGGCGTTATCAc                 <  2:259395/136‑1 (MQ=11)
                                                                                                               gaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                 <  8:64932/135‑1 (MQ=32)
                                                                                                               gaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                 <  3:316071/135‑1 (MQ=32)
                                                                                                                aaCGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAATCCACAACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                 >  3:79472/1‑134 (MQ=12)
                                                                                                                   gCCACTGGCGGAAATTTCGCCCACCTCCCCCCGAACCCGCCCCCGGCTCCGACCACGAACAACGACTCGGGCGATGCCACCCCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTaca            <  1:316967/136‑1 (MQ=11)
                                                                                                                   gCCACTGGCGGAAATTTCGACCACCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTaca            <  8:108504/136‑1 (MQ=255)
                                                                                                                    cccccggCGTTAATTTCGAGCACCTGACCCCGAACCCGACCCCGGCACACACCCGTATCAACGACTCGGGCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTaca            <  4:18286/130‑1 (MQ=11)
                                                                                                                    ccACAGGCGGTAATTTCCACCAACGCCCCCCAACCCCGCCCCCGCCTCACACCACTATCAACGACTCGGGCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTaca            <  4:138307/135‑1 (MQ=11)
                                                                                                                           cGGAAATTTCGACCACCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg       <  4:311218/133‑1 (MQ=12)
                                                                                                                           cGGAAATTTCGACCACCTCACCCCGAACCCGACCCCGGCTCAGACCCCGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg       <  3:128671/133‑1 (MQ=11)
                                                                                                                              tAATTTCCAGCAACTCACCCCTAACCCGCCGCCGCCGCTTAACACCAACAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCt  <  4:111220/135‑1 (MQ=11)
                                                                                                                              |                                                                                                                                       
CAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCT  >  ACB143/187886‑188146

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