Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,350,3600TG37.7% 71.3 / 11.8 53R28R (CGT→CGGCJ019_04090hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (6/27);  new base G (20/0);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.14e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCTACACGGATGTGTGGATCTGCGCCGATCCGCAAGGCCACCTGCAAGCCACC  >  KT2440/4350228‑4350495
                                                                                                                                    |                                                                                                                                        
aGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCGGGCGCGCCCGtttga                                                                                                                                      >  5:339349/1‑135 (MQ=255)
         aTTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGGGCGCCCGGGTGATCTACCTgg                                                                                                                             >  8:281076/1‑137 (MQ=255)
         aTTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGGTTTATCTACCTgg                                                                                                                             >  8:321079/1‑137 (MQ=255)
         aTTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGGTTCATATACCAgg                                                                                                                             >  6:208566/1‑137 (MQ=255)
           tGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGAt                                                                                                                           >  6:224308/1‑137 (MQ=255)
             tCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGc                                                                                                                         <  2:163733/137‑1 (MQ=255)
                   gATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATgg                                                                                                                  <  1:258927/137‑1 (MQ=255)
                   gATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATgg                                                                                                                  <  1:136530/137‑1 (MQ=255)
                     tCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGcc                                                                                                                <  7:119660/137‑1 (MQ=255)
                        gTGACTTGATCGAGGCCCCCCGCCCTTGCTTGACTGCCCGCTGCCGCGCACGCTGGACTACGTTGACCACAGTCAGCGGGGGCTGACCCGTCGGCGCTGGCGCGACCGTTTCTTCTACCTGGATGCCGATGGCCAAc                                                                                                             <  6:45108/137‑1 (MQ=255)
                                   gAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGCCCCTTTCATCTACCTGGATGCCGATGGCCAACGGGTAcgcgcc                                                                                                  >  8:214031/1‑135 (MQ=255)
                                   gAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCGGGCGCGACCGTTTACTCAACCAGGATGCCCATGGCCCACGGGGACGCGAc                                                                                                  >  7:32075/1‑137 (MQ=255)
                                     ggCCCGCCGCCCATGCTGGCCTGCCCCCGGCCGCGCACCCTGCACTACGTAGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACTGGTACGCGACAg                                                                                                <  7:403249/137‑1 (MQ=255)
                                      gcacgccgcCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGCCCGGGTGCTCGACCGGGATGGCGAGGGCCCACAGGGGCGCGGCAGc                                                                                               >  7:172412/4‑137 (MQ=255)
                                        ccgccgCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCGGGCGCGCCCGGTTGCTCTACCTGGATGCCGATGGCCCACGGGTTCGCGCCcgcgc                                                                                             >  7:211162/1‑136 (MQ=255)
                                         cgccgcCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCGGGGGCCCGCGGTTCATATAACAGGATTGCCATGGCCAACCCGTAAGGGACAGCgca                                                                                            >  2:79007/1‑135 (MQ=255)
                                                 aTGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTgc                                                                                    <  1:176191/137‑1 (MQ=255)
                                                      tGAGTGCCCGCTGCCGCGCACCCTGACGTACGTTGACGTGAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACGGTTTCTTCTACCTGGATGCCGTTCGACAACGGGTACGCGACAGCGAAACCCTTgcgcgct                                                                               <  3:372340/137‑2 (MQ=255)
                                                          tGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATcgc                                                                           <  8:31788/137‑1 (MQ=255)
                                                           gCCCCCTGCCGCGCACCCTGCACTAGGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATcgcc                                                                          <  6:48501/137‑1 (MQ=255)
                                                           gCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGCCCGGTTGATCGACCTGGATGACGCAGGGCAACGCGTTCGCGTCAGCGAACACCCTGCGCGCCTcccc                                                                          >  4:339874/1‑137 (MQ=255)
                                                               ccTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGGGGGGGCGCCCGGTTCATGGACCTGGGTGGCGCTGGGCCACGGGGACGCGCCCGCGCAACCCGTGCGCGCATcggcgccg                                                                      >  2:62162/1‑135 (MQ=255)
                                                               ccTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGGGCGCCCGTTTAATCTACCAGGATGCCGATGGCCAACGGGTTCGCGACAGCGAAACCCTTGCGCGCAAGGCCGCaa                                                                      >  6:231259/1‑136 (MQ=255)
                                                                   ccGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGt                                                                  <  5:51872/137‑1 (MQ=255)
                                                                        cACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATcc                                                             <  8:311532/137‑1 (MQ=255)
                                                                         aCCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGCCCGGTTCATCTACCTGGATGCCGAGGGGCCAAGGGTACGCGCCAGCCAAAACCATGCCCGCAACGCCCCATTGGTgtaccc                                                            >  4:52252/1‑132 (MQ=255)
                                                                          cccTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATcccc                                                           <  4:403256/137‑1 (MQ=255)
                                                                          cccTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATcccc                                                           <  7:281076/137‑1 (MQ=255)
                                                                          cccTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATcccc                                                           <  7:59505/137‑1 (MQ=255)
                                                                           ccTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGCCCGGTTAATCTACCTGGATGCCGATTGCCAACGGGTACGCGACAGCGAAACCCATGCGCGCCTCGCCGCATTTGTGGTacccc                                                          >  1:54990/1‑137 (MQ=255)
                                                                            cTGCACTACGTTGACGACAGTCAGCGGGGCCTGACCCGTCGGCGCTGGCGCGACCTTTTCATCTACCTGGATGCCGATGGCCAACGTGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATcccccc                                                         <  1:79007/137‑1 (MQ=255)
                                                                                aCTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCGGGCGCGGCCGGTTAATATACCTGGATGCCGATTGCCAACCGGGACGCGCCAGCGAAACCCATTCGCCCCACGCCCCAATTGTGATCCCCCCGGcc                                                     >  4:296879/1‑137 (MQ=255)
                                                                                 cTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCt                                                    <  7:129710/137‑1 (MQ=255)
                                                                                  tACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGCCCGGTTCATCTTCCGGGATGCCGATGGCCCACGGGTACGCGCCCGCGCAACCCTTGCGCGCCTCGCCGCCTTGGTGATCCCCCCGGCCTa                                                   >  4:281504/1‑137 (MQ=255)
                                                                                  tACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGCCCGGTTAATCTACCTGGATTCCGATGGCCAACCGGTACCCGCCAGCGCAACCCATTCGCCCAAAGCCGCAATGGTGGTACCCCCCGccga                                                   >  4:298289/1‑135 (MQ=255)
                                                                                    ggTTCACGACAGTCAGCGGGGCCTGACCCGGCGGCGGTGGCGCGACCGTTTCATCTACGTGGATCCCGATGCCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCTaca                                                 <  1:397634/136‑1 (MQ=255)
                                                                                    cGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCGTaca                                                 <  7:242455/137‑1 (MQ=255)
                                                                                    cGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCTaca                                                 <  1:243985/137‑1 (MQ=255)
                                                                                    cGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCTaca                                                 <  7:440431/137‑1 (MQ=255)
                                                                                        gacgacAGTCAGCCGGGCCTGACCCGTCGGCGCTGGGGCGGCCGGTTGATCTACCTGGATTCCGATGGCCCACCGGTAAGCGACAGCGAAACCCATGCGCGCATCGCCGCCTTGGTGATCCCCCCCGCCTACACGGa                                             >  2:1130/1‑137 (MQ=255)
                                                                                        gacgacAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGCCCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGCCAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCTACACGGa                                             >  2:81654/1‑137 (MQ=255)
                                                                                        gacgacAGTCAGCCGGGCCTGACCCGTCGGCGCGGGCCGGCCCCGTTGATCTAACTGGAAGCCCATGGCCAAACGGTACCGCCCCCCGAAACCCATGCCCCCCTCGCCCCCATGTTGGACCCCCCCGCCTACACCGa                                             >  8:440431/1‑137 (MQ=255)
                                                                                         acgacAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCTACACGGAt                                            <  6:339349/137‑1 (MQ=255)
                                                                                         acgacAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCTACACGGAt                                            <  5:224308/137‑1 (MQ=255)
                                                                                          cgacAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGGGACCGGTTGATCTACCTGGCTTCCGCTGGGCCACCGGGTCGCCACCGCGAAACCCATGCGCGCCACCGCGCATTGGGGATACCCCCCGCCAAAACGGAtg                                           >  8:230278/1‑137 (MQ=255)
                                                                                          cgacAGTCAGCCGGGCCTGACCCGTCGGCGCGGGCGCGCCCCGTTAATCTACCATGCTTCCGCTTGGCAACGGGGTAGCGACCGCGAAACCCAGGCGCGCATCGCCCCATTTGTGATGCCCCCCGCCAACACGGAtg                                           >  8:264683/1‑137 (MQ=255)
                                                                                            acAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGCCCGGTTGATCGAACGGGCTTGCCATGGCCCACCCGTACGCGCCACCGCAAACCATTCGCGCAACGCCGCATTGGTGATGCCCCCCGCCTACACgcaagtt                                         >  2:32330/1‑131 (MQ=255)
                                                                                                   gCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCTACACGGATGTGTGGATCt                                  <  4:197659/137‑1 (MQ=255)
                                                                                                   gCCGGGCCTGACCCGTCGGCGCGGGCGCGCCCCGGGCATCGAACCGGCTGCCCCTGGCCCACGGGGACGCCCCCGCGCATCCCGTGCGCGCAAAGCCGGCGTGGTGGTGCCCCCGGCCTACCCGGGTGTGTGGGTcg                                  >  6:418563/1‑136 (MQ=255)
                                                                                                            gACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCTACACGGATGTGTGGATCTGCGCCGATc                         <  6:135435/137‑1 (MQ=255)
                                                                                                               ccGTCGGCGCGGTCGCGACCGGTTAATCTTACATGATTCCGATTGCCCAACGGTAAGCGACAGCGCAACCCGATCGCGCATAGACGCATTGGTGCTAGCCCCGGCCTACACGGATGTGTGTGTCTTCGCCGCTACGc                      >  8:241217/1‑137 (MQ=255)
                                                                                                                 gTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCTACACGGATGTGTGGATCTGCGCCGATCCGCaa                    <  6:155092/137‑1 (MQ=255)
                                                                                                                  tCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCTACACGGATGTGTGGATCTGCGCCGATCCGCAAg                   <  8:32075/137‑1 (MQ=255)
                                                                                                                                   gTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCTACACGGATGTGTGGATCTGCGCCGATCCGCAAGGCCACCTGCAAGCCAcc  >  2:345523/1‑137 (MQ=255)
                                                                                                                                    |                                                                                                                                        
AGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGGTACGCGACAGCGAAACCCTTGCGCGCATCGCCGCATTGGTGATCCCCCCGGCCTACACGGATGTGTGGATCTGCGCCGATCCGCAAGGCCACCTGCAAGCCACC  >  KT2440/4350228‑4350495

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: