Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,806,2610TG37.9% 60.4 / 18.5 58V585G (GTG→GGG) priAPrimosomal protein N'
Reads supporting (aligned to +/- strand):  ref base T (23/13);  new base G (0/22);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 2.73e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCGGCACCCAGATGCTCGCCAAGGGCCACCACTTCCCGCGGGTGACCCTGGTGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTACTTCGCCTTCGCCGAACAGGCCCTGGACGAACGCCGA  >  KT2440/5806129‑5806397
                                                                                                                                    |                                                                                                                                        
tcgGCACCCAGATGCTCGCCAAGGGCCACCACTTCCCGCGGGTGACCCTGGTGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTggcc                                                                                                                                      >  4:220549/1‑137 (MQ=255)
tcgGCACCCAGATGCTCGCCAAGGGCCACCACTTCCCGCGGGTGACCCTGGTGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTggcc                                                                                                                                      >  7:62674/1‑137 (MQ=255)
        ccgATGCTCGCCACGGGCCACCACTTCCCGCGGGTGACCCTGGGGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGCGCGCATGGCACAACTTATCGTACAGGTggcggggcgggc                                                                                                                              <  6:197356/135‑1 (MQ=255)
        cAGATGCTCGCCAAGGGCCACCACTTCCCGCGGGTGACCCTGGGGGCCATCCTTGATGCCGATGGGGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTggccgggcgggc                                                                                                                              <  4:430288/137‑1 (MQ=255)
            tGCTCGCCAAGGGCCACCACTTCCCGCGGGTGACCCTGGTGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTggccgggcgggccggg                                                                                                                          <  5:419343/137‑1 (MQ=255)
            tGCTCGCCAAGGGCCACCACTTCCCGCGGGTGACCCTGGGGGCCATCCTTGATGCCGAGGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTggccgggcgggccggg                                                                                                                          <  1:52021/137‑1 (MQ=255)
                cGCCACGGGCCACCACTTCCCGCGGGGGACCCTGGTGGCCAACCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTGCCGCGCCAGCGAGCGCATGGCACAACTTATCGAACAGGGGGccgggcgggccgggcggg                                                                                                                      <  5:311065/137‑1 (MQ=255)
                    aaGGGCCACCACTTCCCGCGGGTGACCCTGGTGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGa                                                                                                                  >  5:118094/1‑137 (MQ=255)
                     aGGGCCACCACTTCCCGCGGGTGACCCTGGTGGCCATCCTTGATGCCTATGGCGGTCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGaa                                                                                                                 <  8:32600/137‑1 (MQ=255)
                             aacTTCCCGCGGTGGACCCTGGTGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGGGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCgg                                                                                                         <  2:89637/136‑1 (MQ=255)
                              acTTCCCGCGGGTGACCCTGGTGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGCGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGCGCACGGc                                                                                                        >  8:96454/1‑137 (MQ=255)
                                  cccGGGGGTCACCGTTGTGCCAATCCTTGATGGCTATGGCGGGCTGTTTACCGGGGACTGCCGCCCCCGCGCCCGCATGGCACAACTTATGGTCCCGGGGGCGGGGCGGGCGGGGCGGGCCGAAGAGCCCGGCAAgg                                                                                                    <  2:157610/137‑1 (MQ=255)
                                  cccGCGGGTGACCCTGGTGGCCATCCTTTATGCCGTGGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCAGGGCACAACTTATGGTACCGGGGGCGGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAgg                                                                                                    <  2:89411/137‑1 (MQ=255)
                                  cccGCGGGTGACCCTGGTGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGGGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAgg                                                                                                    <  6:421878/137‑1 (MQ=255)
                                  cccGCGGGTGACCCTGGTGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGGGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAgg                                                                                                    <  8:303680/137‑1 (MQ=255)
                                    cgcgGGTGACCCTGGTGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGtc                                                                                                  >  4:398367/1‑137 (MQ=255)
                                     gcgGGTGACCCTGGTGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGtca                                                                                                 >  8:215498/1‑137 (MQ=255)
                                     gcgGGTGACCCGGGTGGCCATCCTTGATGCCGATGGCGGGATGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGCGCGGGCCGAAGAGCCCGGCAAGGgca                                                                                                 >  1:155189/1‑137 (MQ=255)
                                          tGACCGTGGTGGCCATCGTTTATGCCTTTGGGGGGGTGTTTTCGGTTTACTGCCGCGCCAGCGAGCGCTTGGCAAAAATTTTCTTACTGGGGGCGGGGGGGGCGGGGGGGGCCGAATAGCCCGGCAAGGTCatcatc                                                                                            <  6:184299/137‑1 (MQ=255)
                                          tGACCCGGGTGGCAATCCTTGAGGCCGAGGGGGGGTTGTTTTCGGTTGACTTCCGCGCCAGCGAGCGAATGGAACAAATTATCGTACCGGGGGCGGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCatcatc                                                                                            <  8:236133/137‑1 (MQ=255)
                                               cTGGTGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGAc                                                                                       >  4:429941/1‑137 (MQ=255)
                                                    ggCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCAcc                                                                                  >  7:341204/1‑137 (MQ=255)
                                                    ggCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCAcc                                                                                  >  3:155944/1‑137 (MQ=255)
                                                     gCCTTCCTTGATGCGGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCTGCGAGCGCAGGGCACTACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGTGCCCGGCAATGTCATCATCCAGACCCACCt                                                                                 <  6:172162/137‑1 (MQ=255)
                                                     gCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCt                                                                                 >  2:220200/1‑137 (MQ=255)
                                                      ccATCCTTGATGAATATGTAAAGCTGTTATCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGAccctgtct                                                                                >  5:208891/1‑132 (MQ=255)
                                                        aTCCTTGATGCCGAGGGCGGGCTGTTTTCCGGTGACTGCCGCGCCAGCGCGCGCATGGCACAACTTATCGTCCAGGGGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGc                                                                              <  1:276065/137‑1 (MQ=255)
                                                           cTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCgc                                                                           >  7:464067/1‑136 (MQ=255)
                                                           cTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGa                                                                           >  6:304602/1‑137 (MQ=255)
                                                             tGATGCCGATGGGGGGCTGTTTTCGGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCTTACCGGGGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGacc                                                                         <  7:369696/137‑1 (MQ=255)
                                                             tGAGGCCGATGGGGGGCTGTTGTGCGGGCACTCCCGCGCCAGCGCGCGCATGGCACAAATTCTCGTCCCGGGGGCGGGGCGGGCGGGGGGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGacc                                                                         <  8:437556/137‑1 (MQ=255)
                                                              caggCCGATGGCGGGGTGTTTTCCGGGTACTGCCGCGCCAGGGAGCGATTGCAACAATTTATCGTACAGGGGGCGGGGCGGGCGGGGGGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGacca                                                                        <  7:161947/134‑1 (MQ=255)
                                                                  ccGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGAccaccca                                                                    >  3:373560/1‑137 (MQ=255)
                                                                     aTGGCGGGCTGTTTTCCGGGGACTGCCCCGCCAGCGAGCGCATGGCACAACTTATCGTACCGGGGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTg                                                                 <  4:119880/137‑1 (MQ=255)
                                                                        gCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCCGACCCACCTGGCCGACCACCCACTGTTg                                                              >  2:111341/1‑137 (MQ=255)
                                                                        gCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTg                                                              >  7:430782/1‑137 (MQ=255)
                                                                          gtgtttttttCCGGTGCGTTCCGCGCCGGCTCGCGCTGGTCAAAATTTATCGTACGGGGGGCGGGGGGGGGCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGt                                                            <  8:430782/131‑1 (MQ=255)
                                                                             cTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCa                                                         >  7:276382/1‑137 (MQ=255)
                                                                                tttgcgggTGCCTGCCGCGCCCGCGATCGCATGGCCAAAATTAGCGTCCCGGGGGCGGGGGGGGCGGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCt                                                      <  1:393269/131‑1 (MQ=255)
                                                                                    ccGGGGCCTCCCGCGCCAGCGCGCGCATGCCACAACTTATCGACCTGGGGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAAcc                                                  <  2:143915/137‑1 (MQ=255)
                                                                                     cGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCg                                                 <  1:277452/137‑1 (MQ=255)
                                                                                         cacagccGCCCCAGCGAGCGCATGGAACAAATTATCGTACCGTGGGCGGGGCGGGCGGGGGGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCa                                             <  2:336048/132‑1 (MQ=255)
                                                                                             tCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGc                                         >  4:227557/1‑137 (MQ=255)
                                                                                             tCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGc                                         >  8:358604/1‑137 (MQ=255)
                                                                                               cgcgCCAGCGAGCGCATGGCACAAATTATCTTCCCGGGGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTa                                       <  2:326851/137‑1 (MQ=255)
                                                                                                  gcCAGCGCGCGCATGCAACAACTTATCGTACAGGGGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTActt                                    <  6:83451/137‑1 (MQ=255)
                                                                                                  gcCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTActt                                    <  3:220549/137‑1 (MQ=255)
                                                                                                       cgcgcgcTTGTCACAACTTATCGTCCCGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGGCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTActtcgcct                               <  6:118094/134‑1 (MQ=255)
                                                                                                       cgagcgCATGGCACAACTTATCGTCCAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTActtcgcct                               <  5:450216/137‑1 (MQ=255)
                                                                                                            gcATGTCACAACTTCTCGTCCTGGGGGCGGGGCGGGCGGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTACttcgccttcgcc                          <  7:452823/137‑1 (MQ=255)
                                                                                                            gcATGGCACACCTTTTCGTTCAGGTGGCGGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGGCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTACttcgccttcgcc                          <  3:349176/137‑1 (MQ=255)
                                                                                                            gcATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTACttcgccttcgcc                          <  6:166743/137‑1 (MQ=255)
                                                                                                             cATGGCACAACTTATCGTCCCGGGGCCCGGGCGGGCGGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTACTTCGCCTTCGCCg                         <  3:356834/137‑1 (MQ=255)
                                                                                                              agggCAAAACTAATCGTCCCGGGGGCGGGGGGGGCGGGGGGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTACTTCGCCTTCGCCGa                        <  1:452277/135‑1 (MQ=255)
                                                                                                              aTGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTACTTCGCCTTCGCCGa                        <  1:15513/137‑1 (MQ=255)
                                                                                                                             gTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTACTTCGCCTTCGCCGAACAGGCCCTGGACGa         >  5:326975/1‑137 (MQ=255)
                                                                                                                               aCAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTACTTCGCCTTCGCCGAACAGGCCCTGGACGAAc       >  8:186/1‑137 (MQ=255)
                                                                                                                               aCAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTACTTCGCCTTCGCCGAACAGGCCCTGGACGAAc       >  7:161092/1‑137 (MQ=255)
                                                                                                                                    tGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTACTTCGCCTTCGCCGAACAGGCCCTGGACGAACGCCgc  >  6:237268/1‑136 (MQ=255)
                                                                                                                                    tGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACAACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTACTTCGCCTTCGCCGAACAGGCCCTGGACGAACGCCGa  >  3:395266/1‑137 (MQ=255)
                                                                                                                                    |                                                                                                                                        
TCGGCACCCAGATGCTCGCCAAGGGCCACCACTTCCCGCGGGTGACCCTGGTGGCCATCCTTGATGCCGATGGCGGGCTGTTTTCCGGTGACTTCCGCGCCAGCGAGCGCATGGCACAACTTATCGTACAGGTGGCCGGGCGGGCCGGGCGGGCCGAAGAGCCCGGCAAGGTCATCATCCAGACCCACCTGGCCGACCACCCACTGTTGGTCCAGCTAACCGAGCAAGGCTACTTCGCCTTCGCCGAACAGGCCCTGGACGAACGCCGA  >  KT2440/5806129‑5806397

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: