Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,164,3630AC22.5% 122.8 / 10.7 71G537G (GGT→GGG) parCDNA topoisomerase IV subunit A
Reads supporting (aligned to +/- strand):  ref base A (37/18);  new base C (1/15);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.96e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGTTCTTTCTTCAGCAG  >  NC_000913/3164113‑3164596
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          
aaGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTGCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTa                                                                                                                                                                                                                                            <  2:51181/251‑1 (MQ=255)
    cTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGAT‑GCTCTTACCtt                                                                                                                                                                                                                                        >  1:273707/1‑251 (MQ=255)
      tACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACCGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGGATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTc                                                                                                                                                                                                                                      >  2:368333/1‑251 (MQ=255)
                     gcgcCACCAGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑CCGCTTCCCTGTCACCGCCGCTTTGAAg                                                                                                                                                                                                                       <  1:363794/251‑1 (MQ=255)
                        ccaccaGTCCGTTAAAGGGGCAGACGAAACCGTAACCCGCACCGGAAGCCATCAGCAGTTTCTGACCGTCGCTTTCCATCAGCATATGGTCAACGGCCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCCCGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACCCTACCGGTGGGTT‑GCTCTTACCTTCCCCCGCCGCTTTGAAGCTa                                                                                                                                                                                                                    <  2:196065/251‑1 (MQ=255)
                         caccaGATCGTTAAAGGTGCAGACGAAACCGTACCCCGCATCGGAGCCCATCAGCAGTTTCTGATCGGCCCTTTCCACCAGCATAGGGACAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGGGAGCGGCTCGCCCTGACCCCGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAACCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTAt                                                                                                                                                                                                                   <  2:382072/251‑1 (MQ=255)
                          accaGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATc                                                                                                                                                                                                                  >  1:25239/1‑251 (MQ=255)
                           ccaGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCCCGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGCATCAACAAACACTACCGGTTGGGT‑GCTCTTACCTTTCACCGCCCCtttcacgctttcc                                                                                                                                                                                                                 >  2:93228/1‑241 (MQ=255)
                           ccaGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATcc                                                                                                                                                                                                                 >  1:157429/1‑250 (MQ=255)
                            caGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGGTGGTT‑GCTCTTACCTTTCCCCGCCGCTTTGAAGCTATCAc                                                                                                                                                                                                                >  2:57693/1‑251 (MQ=255)
                               aTCGTTAAAGGTGCAGACGAACCCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTGCACCGCCGCTTTGAAGCTATCACCCg                                                                                                                                                                                                             <  1:393961/251‑1 (MQ=255)
                                tCGTTAAAGGTGCAGACGATACCGTAACCCGCATCGGAACCCATCAGCAGCTTCTGATCGTCGCTTTCCATAAGCATATGGACAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGTTGAGCGCAACGCCCTGGCCACGCGCCGACGGAAGCGTAATCGGGTAAATGGCATAGCTCCCGCCAGTGAAATCAACCAACCCTACCGGTGGGTT‑GCTCTTATCTTTCATCGCCACTTTgcatcgcccacccgc                                                                                                                                                                                                            >  2:271935/1‑237 (MQ=255)
                                           gaaGACGAAACCGTACCCCGCATCGGAAGCCATCAGCAGTTGCTGATCGCCGCTTCCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCGGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCa                                                                                                                                                                                                 <  2:268115/249‑1 (MQ=255)
                                            cAGACGAACCCGTAACCCGCATCGGAGGCCATCAGCAGTTTTTGATCGTCGCTTCCCATCAGCATAGGGTCAACGGTCGCCCGAGGCGGCAACGTAAATTTGCGGGTGAGCGGCTCGCCCTGACCCCGCGCCGACGGCAGCGTAATCGGGTCAATGGCAAAGCTACGGCCGGTGGAATCAACAAACCCTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAg                                                                                                                                                                                                <  1:106873/251‑1 (MQ=255)
                                             aGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAgg                                                                                                                                                                                               >  1:302816/1‑251 (MQ=255)
                                                 gAAACCGTAACCCGCACCGGCAGCCATCAGCCGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCCCGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTAACGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCCGGCCCg                                                                                                                                                                                           >  2:285231/1‑251 (MQ=255)
                                                      cGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGGGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGcgcg                                                                                                                                                                                      >  1:262906/1‑251 (MQ=255)
                                                           cccGCATCGGAGGCCATCAGCAGTTTCTGATCGTCGCTGTCCATCAGCATAGGGTCAAGGGTCGCCCAAGGCGGAACCGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGAAGCGTAATCGGTTAAATGGCATCGCTACGACCGGTGGAATCAACAAACACTACCGTTGGGTT‑GCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGat                                                                                                                                                                                 <  1:376490/251‑1 (MQ=255)
                                                              gCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGAAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACACGCTTTATAATTCAGGCCCGGCGCGTCGATATc                                                                                                                                                                              >  1:435536/1‑251 (MQ=255)
                                                                   ggAAGCCATCAGCAGTTCCTGACCGTCGCTTTCCATCGGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTGGCCGGGGAGCGGCTCGCCCTGACCACGGGCCGACGGCAGCGTAAGCGGTTAAATGGCATAGCGACGGCCGGGGGAATCAACAAACATTGCCGGTTGGTT‑GCTCTTACGTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGc                                                                                                                                                                         <  1:221389/251‑1 (MQ=255)
                                                                      aGCCATCAGCAGTTTCGGATCGTCGCTTGCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGGGAGCGGCTCGCCCTGACCCCGCGCCTACGGCAGCGTAATCGGGTCAATGGCATAGCTACGCCCGGGGGAATCAACAAACACTACCGGTTGGTT‑GCTCTGCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCtt                                                                                                                                                                      <  1:307956/251‑1 (MQ=255)
                                                                                  tttCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGt                                                                                                                                                           >  1:383911/1‑250 (MQ=255)
                                                                                                 ttccacaaGCATATGACCACCGGCCGCCCCAGGCGACAACGTTAATTTGCCGGGGGGCGGCTCCCCCTGCCCCCGCCCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGAAATCAACAAACACTCCCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTg                                                                                                                                             <  2:273707/242‑1 (MQ=255)
                                                                                                             ggggccaccGGTCGCCCCGGGCGGCACGGTGAATTTGCGGGGGGGCGGCTGCCCCTGACCCCGCGCCGCCGGCGGCGTAATCGGGTCATTGGCATAGCTACGACCGGTGGAACAAACAACCCCTCCCGTTGGGTT‑GCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATgg                                                                                                                               <  2:342913/243‑1 (MQ=255)
                                                                                                             atGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCCCAATgg                                                                                                                               >  2:467222/1‑251 (MQ=255)
                                                                                                                   aacgggcgaccaaGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGACGCCGGCCGCGTTATCCGGTCACTGGCATAGCTACGACCGGGGGACTCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACTGCCGCTTTGAAGCTATCACCCGCTTTATACTTCAGGCCAGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGACCATCTGCGACAGCACAATGGTGA‑c                                                                                                                          >  2:48676/13‑249 (MQ=255)
                                                                                                                      ggTCGCCCCAGGCGGCTACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCCCAATGGTGA‑CAGGtt                                                                                                                     >  2:405512/1‑251 (MQ=255)
                                                                                                                               aGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGGTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGgc                                                                                                             >  1:332806/1‑250 (MQ=255)
                                                                                                                               aGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCCGACCAAGCGCCCACCGCAGCGTACTCGGGTCAATGGCATAGCTACGCCCGGTGGACTCAACAAACACTACCGGTTGGTT‑GCCCTTACCTTTCACCGCCGCTTTGACGCTATCACCCGCTTTATAATTCAGGCCCCGCGCCTCGCTCTCATAGCCTTTACCGCTGCGTACCCCGCCCATCTGCGCCAGCCCCATGGGGC‑CAGGTTCaaccggc                                                                                                             >  2:288411/1‑244 (MQ=255)
                                                                                                                               aGGCGGCAACGTCAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACCGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCCGCGACAGCACCATGGTGA‑AAGGGTCCGACGgcc                                                                                                            >  2:250625/1‑250 (MQ=255)
                                                                                                                                  cggcAACGTTAATTTGCCGGTGAGCGGCGCGCCCTGACCACGCGCCCACCACCGCGTAATCGGGTCAATGGCCTAGCAACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGGTTTATAATTAAGGCCCGGCGCGTCGATTTCCTGGCCTTTCGCGCTGCGTACCCAGCCCCTCTGCCACAGCACAACGGTGC‑CAGGTTCCGGCGgcagc                                                                                                          >  2:233531/1‑250 (MQ=255)
                                                                                                                                        cGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATCGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCCCCCGCTTTATAATTCACGCCCGGCGCGCCGATATCATGGCCTTTATCGCTGCGGACCCAACCCATCTGCGACAGCACAAGGGTGAACAG‑TTCAGACGGCAGCATGTCGt                                                                                                   >  2:478848/1‑251 (MQ=255)
                                                                                                                                                         gCGGCTCGCCCTGACAACGCGCCGACGCCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACCAACACTACCGGTTGGTT‑GCTCTTCCCTTTCACCGCCGCTTTGAAGATATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCCTCTGCGACAGCACAATGGTGA‑AAGGGTCCGACGGCAGCATGTCGTGCTCGCTCAtagctttc                                                                                  >  2:39438/1‑251 (MQ=255)
                                                                                                                                                          cGGCTCGCCCGGACCACGCGCCGACGGCGGCGTAATCGGGTCAATGGCATAGCTACGACGGGTGAAATCAACAAACACTCCCGTTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCAtcgctttcg                                                                                 <  1:120485/251‑1 (MQ=255)
                                                                                                                                                           ggCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCAtcgctttcgc                                                                                >  2:240258/1‑251 (MQ=255)
                                                                                                                                                              ccccccctaCCACGCACCGACGGAAGCGTAATCGGGTTAATGGCATAGCTACGCCCGGGGGAACCAACAACCCCTCCCGGTTGTTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGGACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGcttc                                                                             <  2:157429/243‑1 (MQ=255)
                                                                                                                                                                        cgcccgccgcggcaagggTAATCGGGTCAAGGGCATACCTCCGCCCGGGGGAACCAACAACCACTCCCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTcc                                                                   <  2:99356/234‑1 (MQ=255)
                                                                                                                                                                         ccccgccgcCGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACCCTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCt                                                                  <  1:287602/248‑1 (MQ=255)
                                                                                                                                                                            cgcCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGc                                                                >  1:360217/1‑250 (MQ=255)
                                                                                                                                                                                 aCGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGc                                                          >  2:86549/1‑251 (MQ=255)
                                                                                                                                                                                         gTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTCGTT‑GCTCTTACCTTTCACCGCCGATTTGAAGCTATCACCAGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCACACCATCTGCGACAGCACACTGGTGA‑CAGGTTCCGCCGGCAGCATGTCGTGCTCGCTCCTCACTTTCGCTTCTTCGCGTTCCTGCAACGGCGAccgacg                                                  >  2:69304/1‑251 (MQ=255)
                                                                                                                                                                                          tAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAacgacg                                                  >  1:59129/1‑250 (MQ=255)
                                                                                                                                                                                           aaTCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGAt                                                >  2:213074/1‑251 (MQ=255)
                                                                                                                                                                                           aaTCGGGTCAATCGCATAGCTACGCACGGCGGACTCAACAAAAAATACAGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCCGGCCCAGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGAGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGAGAACGGCGAt                                                >  2:218998/1‑251 (MQ=255)
                                                                                                                                                                                            aacGGGTCAAGGGCATAGCTACGACCGGTGGAATCAACAACCACTCCCGTTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATc                                               <  1:218924/249‑1 (MQ=255)
                                                                                                                                                                                             gcgggtcaagagcagacacacccccGGTGGAAACAACACACCCCGCCGGTTGGTT‑GCTCTTACCTTTCCCCGCCGCTTTGAAGCGATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATcc                                              <  2:262906/228‑2 (MQ=255)
                                                                                                                                                                                               gggTCAATGGCATAGCTCCGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTc                                            <  1:368333/251‑1 (MQ=255)
                                                                                                                                                                                                 gTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACCAACGGGTGGTT‑GCTCTTACCTTTCACCGCAGCTTTGAAGCTATAACCCGCTTTATAATTCAGGCACGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCCTCTGCGACAGAACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCAc                                          >  1:290514/1‑251 (MQ=255)
                                                                                                                                                                                                                  gACCGGTGGAATCAACAAACCCTCCCGGTTGGTT‑GCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTct                         <  1:205698/251‑1 (MQ=255)
                                                                                                                                                                                                                  gACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTct                         >  1:418190/1‑251 (MQ=255)
                                                                                                                                                                                                                  gACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGACGCTTTGAAGCTATCACCAGCTTTATAATTAAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTct                         >  2:429142/1‑251 (MQ=255)
                                                                                                                                                                                                                  caCCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCGTTATAATTCAGTCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTGCTTCGCGTCCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTct                         <  1:483930/250‑1 (MQ=255)
                                                                                                                                                                                                                     cGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTctgcc                      >  1:129481/1‑251 (MQ=255)
                                                                                                                                                                                                                       ggggaaaaaacaacccctccgggttggtt‑GCTCTACCCTTCCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTctgcctg                    <  2:301723/230‑1 (MQ=255)
                                                                                                                                                                                                                        gggattCAACAAACCCTACCGTTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTctgcctgc                   <  1:240258/246‑1 (MQ=255)
                                                                                                                                                                                                                        ggaaataaacaaccacccccGGGTGGTT‑GCTCTTCCCTTCCCCCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTctgcctgc                   <  2:59129/233‑1 (MQ=255)
                                                                                                                                                                                                                          aaaaaaacaaacaCT‑CCCGTTGGTTTCCTCTTCCCTTCCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAg                 <  2:418190/246‑1 (MQ=255)
                                                                                                                                                                                                                           aaaaaaaaaaccCTACCGTTGGGTT‑CCTCTTACCTTTCACCGGCGCTTTGAAGCTATCACCCGCTTGAAAAGTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTCCCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGGGAACGACGATCGTCACCGTAGGCTTGCGCGCCTGCCTGCAGt                <  1:408619/244‑1 (MQ=255)
                                                                                                                                                                                                                            aacaacaaCCACGACCTGTTGGGT‑GCTCTGCCCTTTCTCCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATTTCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCTTGCTCGCTCATCGCTTTCGCTTCGTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGtt               <  1:216785/249‑1 (MQ=255)
                                                                                                                                                                                                                             ccaacaaACACTACCGGTTGGTT‑GCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGGCGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttc              <  1:429142/250‑1 (MQ=255)
                                                                                                                                                                                                                             acccaaaccacccccGGTTGGTT‑GCTTTTCCCTTTCCCCGCCGCTTTGAAGCTATCACCCGCTGTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttc              <  2:135667/238‑1 (MQ=255)
                                                                                                                                                                                                                               aacaaccaCTCCCGTTTGGTT‑GCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGTCCTTTAGCGCTGCGTACCCAGCCCATCTGCGGCAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttctt            <  1:254162/251‑1 (MQ=255)
                                                                                                                                                                                                                               aacaaacaCTACCGTTTGGTT‑GCTCTTACCTTCCCCCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCAGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttctt            <  1:47316/251‑1 (MQ=255)
                                                                                                                                                                                                                                acaaacaCTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttcttt           >  2:190087/1‑251 (MQ=255)
                                                                                                                                                                                                                                acaaacaCTACCGGTTGGTT‑GATCTTACCTTTAACCGCCGCTTTGAAGATACCACCCGCTTTATAATTCCGGCCCGGAGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCCGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAAAGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttcttt           >  2:481074/1‑251 (MQ=255)
                                                                                                                                                                                                                                 caaacaCTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttctttc          >  2:45357/1‑251 (MQ=255)
                                                                                                                                                                                                                                 caaacaCTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttc              >  1:135667/1‑247 (MQ=255)
                                                                                                                                                                                                                                  aaacaCTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttctttct         >  1:467323/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aaccaccccggtggttttctctTCCCTTCCCCCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGTtctttctt        <  2:129481/233‑1 (MQ=255)
                                                                                                                                                                                                                                         cccGGTTGTTT‑GCTCTCACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGTTCTTTCTTcagcag  <  1:289093/250‑1 (MQ=255)
                                                                                                                                                                                                                                          ccGGTTGGTT‑GCTCTTACTTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATc                                               <  1:46351/205‑1 (MQ=255)
                                                                                                                                                                                                                                          ccGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATc                                               >  2:46351/1‑205 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          
AAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTT‑GCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGA‑CAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGTTCTTTCTTCAGCAG  >  NC_000913/3164113‑3164596

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