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→GGGparCDNA 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.

AAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGTTCTTTCTTCAGCAG  >  NC_000913/3164113‑3164596
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          
aaGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTGCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTa                                                                                                                                                                                                                                            <  2:51181/251‑1 (MQ=255)
    cTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGATGCTCTTACCtt                                                                                                                                                                                                                                        >  1:273707/1‑251 (MQ=255)
      tACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACCGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGGATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTc                                                                                                                                                                                                                                      >  2:368333/1‑251 (MQ=255)
                     gcgcCACCAGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTCCGCTTCCCTGTCACCGCCGCTTTGAAg                                                                                                                                                                                                                       <  1:363794/251‑1 (MQ=255)
                        ccaccaGTCCGTTAAAGGGGCAGACGAAACCGTAACCCGCACCGGAAGCCATCAGCAGTTTCTGACCGTCGCTTTCCATCAGCATATGGTCAACGGCCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCCCGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACCCTACCGGTGGGTTGCTCTTACCTTCCCCCGCCGCTTTGAAGCTa                                                                                                                                                                                                                    <  2:196065/251‑1 (MQ=255)
                         caccaGATCGTTAAAGGTGCAGACGAAACCGTACCCCGCATCGGAGCCCATCAGCAGTTTCTGATCGGCCCTTTCCACCAGCATAGGGACAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGGGAGCGGCTCGCCCTGACCCCGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAACCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTAt                                                                                                                                                                                                                   <  2:382072/251‑1 (MQ=255)
                          accaGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATc                                                                                                                                                                                                                  >  1:25239/1‑251 (MQ=255)
                           ccaGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCCCGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGCATCAACAAACACTACCGGTTGGGTGCTCTTACCTTTCACCGCCCCtttcacgctttcc                                                                                                                                                                                                                 >  2:93228/1‑241 (MQ=255)
                           ccaGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATcc                                                                                                                                                                                                                 >  1:157429/1‑250 (MQ=255)
                            caGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGGTGGTTGCTCTTACCTTTCCCCGCCGCTTTGAAGCTATCAc                                                                                                                                                                                                                >  2:57693/1‑251 (MQ=255)
                               aTCGTTAAAGGTGCAGACGAACCCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTGCACCGCCGCTTTGAAGCTATCACCCg                                                                                                                                                                                                             <  1:393961/251‑1 (MQ=255)
                                tCGTTAAAGGTGCAGACGATACCGTAACCCGCATCGGAACCCATCAGCAGCTTCTGATCGTCGCTTTCCATAAGCATATGGACAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGTTGAGCGCAACGCCCTGGCCACGCGCCGACGGAAGCGTAATCGGGTAAATGGCATAGCTCCCGCCAGTGAAATCAACCAACCCTACCGGTGGGTTGCTCTTATCTTTCATCGCCACTTTgcatcgcccacccgc                                                                                                                                                                                                            >  2:271935/1‑237 (MQ=255)
                                           gaaGACGAAACCGTACCCCGCATCGGAAGCCATCAGCAGTTGCTGATCGCCGCTTCCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCGGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCa                                                                                                                                                                                                 <  2:268115/249‑1 (MQ=255)
                                            cAGACGAACCCGTAACCCGCATCGGAGGCCATCAGCAGTTTTTGATCGTCGCTTCCCATCAGCATAGGGTCAACGGTCGCCCGAGGCGGCAACGTAAATTTGCGGGTGAGCGGCTCGCCCTGACCCCGCGCCGACGGCAGCGTAATCGGGTCAATGGCAAAGCTACGGCCGGTGGAATCAACAAACCCTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAg                                                                                                                                                                                                <  1:106873/251‑1 (MQ=255)
                                             aGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAgg                                                                                                                                                                                               >  1:302816/1‑251 (MQ=255)
                                                 gAAACCGTAACCCGCACCGGCAGCCATCAGCCGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCCCGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTAACGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCCGGCCCg                                                                                                                                                                                           >  2:285231/1‑251 (MQ=255)
                                                      cGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGGGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGcgcg                                                                                                                                                                                      >  1:262906/1‑251 (MQ=255)
                                                           cccGCATCGGAGGCCATCAGCAGTTTCTGATCGTCGCTGTCCATCAGCATAGGGTCAAGGGTCGCCCAAGGCGGAACCGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGAAGCGTAATCGGTTAAATGGCATCGCTACGACCGGTGGAATCAACAAACACTACCGTTGGGTTGCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGat                                                                                                                                                                                 <  1:376490/251‑1 (MQ=255)
                                                              gCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGAAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACACGCTTTATAATTCAGGCCCGGCGCGTCGATATc                                                                                                                                                                              >  1:435536/1‑251 (MQ=255)
                                                                   ggAAGCCATCAGCAGTTCCTGACCGTCGCTTTCCATCGGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTGGCCGGGGAGCGGCTCGCCCTGACCACGGGCCGACGGCAGCGTAAGCGGTTAAATGGCATAGCGACGGCCGGGGGAATCAACAAACATTGCCGGTTGGTTGCTCTTACGTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGc                                                                                                                                                                         <  1:221389/251‑1 (MQ=255)
                                                                      aGCCATCAGCAGTTTCGGATCGTCGCTTGCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGGGAGCGGCTCGCCCTGACCCCGCGCCTACGGCAGCGTAATCGGGTCAATGGCATAGCTACGCCCGGGGGAATCAACAAACACTACCGGTTGGTTGCTCTGCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCtt                                                                                                                                                                      <  1:307956/251‑1 (MQ=255)
                                                                                  tttCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGt                                                                                                                                                           >  1:383911/1‑250 (MQ=255)
                                                                                                 ttccacaaGCATATGACCACCGGCCGCCCCAGGCGACAACGTTAATTTGCCGGGGGGCGGCTCCCCCTGCCCCCGCCCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGAAATCAACAAACACTCCCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTg                                                                                                                                             <  2:273707/242‑1 (MQ=255)
                                                                                                             ggggccaccGGTCGCCCCGGGCGGCACGGTGAATTTGCGGGGGGGCGGCTGCCCCTGACCCCGCGCCGCCGGCGGCGTAATCGGGTCATTGGCATAGCTACGACCGGTGGAACAAACAACCCCTCCCGTTGGGTTGCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATgg                                                                                                                               <  2:342913/243‑1 (MQ=255)
                                                                                                             atGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCCCAATgg                                                                                                                               >  2:467222/1‑251 (MQ=255)
                                                                                                                   aacgggcgaccaaGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGACGCCGGCCGCGTTATCCGGTCACTGGCATAGCTACGACCGGGGGACTCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACTGCCGCTTTGAAGCTATCACCCGCTTTATACTTCAGGCCAGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGACCATCTGCGACAGCACAATGGTGAc                                                                                                                          >  2:48676/13‑249 (MQ=255)
                                                                                                                      ggTCGCCCCAGGCGGCTACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCCCAATGGTGACAGGtt                                                                                                                     >  2:405512/1‑251 (MQ=255)
                                                                                                                               aGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGGTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGgc                                                                                                             >  1:332806/1‑250 (MQ=255)
                                                                                                                               aGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCCGACCAAGCGCCCACCGCAGCGTACTCGGGTCAATGGCATAGCTACGCCCGGTGGACTCAACAAACACTACCGGTTGGTTGCCCTTACCTTTCACCGCCGCTTTGACGCTATCACCCGCTTTATAATTCAGGCCCCGCGCCTCGCTCTCATAGCCTTTACCGCTGCGTACCCCGCCCATCTGCGCCAGCCCCATGGGGCCAGGTTCaaccggc                                                                                                             >  2:288411/1‑244 (MQ=255)
                                                                                                                               aGGCGGCAACGTCAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACCGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCCGCGACAGCACCATGGTGAAAGGGTCCGACGgcc                                                                                                            >  2:250625/1‑250 (MQ=255)
                                                                                                                                  cggcAACGTTAATTTGCCGGTGAGCGGCGCGCCCTGACCACGCGCCCACCACCGCGTAATCGGGTCAATGGCCTAGCAACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGGTTTATAATTAAGGCCCGGCGCGTCGATTTCCTGGCCTTTCGCGCTGCGTACCCAGCCCCTCTGCCACAGCACAACGGTGCCAGGTTCCGGCGgcagc                                                                                                          >  2:233531/1‑250 (MQ=255)
                                                                                                                                        cGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATCGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCCCCCGCTTTATAATTCACGCCCGGCGCGCCGATATCATGGCCTTTATCGCTGCGGACCCAACCCATCTGCGACAGCACAAGGGTGAACAGTTCAGACGGCAGCATGTCGt                                                                                                   >  2:478848/1‑251 (MQ=255)
                                                                                                                                                         gCGGCTCGCCCTGACAACGCGCCGACGCCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACCAACACTACCGGTTGGTTGCTCTTCCCTTTCACCGCCGCTTTGAAGATATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCCTCTGCGACAGCACAATGGTGAAAGGGTCCGACGGCAGCATGTCGTGCTCGCTCAtagctttc                                                                                  >  2:39438/1‑251 (MQ=255)
                                                                                                                                                          cGGCTCGCCCGGACCACGCGCCGACGGCGGCGTAATCGGGTCAATGGCATAGCTACGACGGGTGAAATCAACAAACACTCCCGTTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCAtcgctttcg                                                                                 <  1:120485/251‑1 (MQ=255)
                                                                                                                                                           ggCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCAtcgctttcgc                                                                                >  2:240258/1‑251 (MQ=255)
                                                                                                                                                              ccccccctaCCACGCACCGACGGAAGCGTAATCGGGTTAATGGCATAGCTACGCCCGGGGGAACCAACAACCCCTCCCGGTTGTTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGGACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGcttc                                                                             <  2:157429/243‑1 (MQ=255)
                                                                                                                                                                        cgcccgccgcggcaagggTAATCGGGTCAAGGGCATACCTCCGCCCGGGGGAACCAACAACCACTCCCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTcc                                                                   <  2:99356/234‑1 (MQ=255)
                                                                                                                                                                         ccccgccgcCGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACCCTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCt                                                                  <  1:287602/248‑1 (MQ=255)
                                                                                                                                                                            cgcCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGc                                                                >  1:360217/1‑250 (MQ=255)
                                                                                                                                                                                 aCGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGc                                                          >  2:86549/1‑251 (MQ=255)
                                                                                                                                                                                         gTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTCGTTGCTCTTACCTTTCACCGCCGATTTGAAGCTATCACCAGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCACACCATCTGCGACAGCACACTGGTGACAGGTTCCGCCGGCAGCATGTCGTGCTCGCTCCTCACTTTCGCTTCTTCGCGTTCCTGCAACGGCGAccgacg                                                  >  2:69304/1‑251 (MQ=255)
                                                                                                                                                                                          tAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAacgacg                                                  >  1:59129/1‑250 (MQ=255)
                                                                                                                                                                                           aaTCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGAt                                                >  2:213074/1‑251 (MQ=255)
                                                                                                                                                                                           aaTCGGGTCAATCGCATAGCTACGCACGGCGGACTCAACAAAAAATACAGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCCGGCCCAGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGAGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGAGAACGGCGAt                                                >  2:218998/1‑251 (MQ=255)
                                                                                                                                                                                            aacGGGTCAAGGGCATAGCTACGACCGGTGGAATCAACAACCACTCCCGTTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATc                                               <  1:218924/249‑1 (MQ=255)
                                                                                                                                                                                             gcgggtcaagagcagacacacccccGGTGGAAACAACACACCCCGCCGGTTGGTTGCTCTTACCTTTCCCCGCCGCTTTGAAGCGATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATcc                                              <  2:262906/228‑2 (MQ=255)
                                                                                                                                                                                               gggTCAATGGCATAGCTCCGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTc                                            <  1:368333/251‑1 (MQ=255)
                                                                                                                                                                                                 gTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACCAACGGGTGGTTGCTCTTACCTTTCACCGCAGCTTTGAAGCTATAACCCGCTTTATAATTCAGGCACGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCCTCTGCGACAGAACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCAc                                          >  1:290514/1‑251 (MQ=255)
                                                                                                                                                                                                                  gACCGGTGGAATCAACAAACCCTCCCGGTTGGTTGCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTct                         <  1:205698/251‑1 (MQ=255)
                                                                                                                                                                                                                  gACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTct                         >  1:418190/1‑251 (MQ=255)
                                                                                                                                                                                                                  gACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGACGCTTTGAAGCTATCACCAGCTTTATAATTAAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTct                         >  2:429142/1‑251 (MQ=255)
                                                                                                                                                                                                                  caCCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCGTTATAATTCAGTCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTGCTTCGCGTCCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTct                         <  1:483930/250‑1 (MQ=255)
                                                                                                                                                                                                                     cGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTctgcc                      >  1:129481/1‑251 (MQ=255)
                                                                                                                                                                                                                       ggggaaaaaacaacccctccgggttggttGCTCTACCCTTCCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTctgcctg                    <  2:301723/230‑1 (MQ=255)
                                                                                                                                                                                                                        gggattCAACAAACCCTACCGTTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTctgcctgc                   <  1:240258/246‑1 (MQ=255)
                                                                                                                                                                                                                        ggaaataaacaaccacccccGGGTGGTTGCTCTTCCCTTCCCCCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTctgcctgc                   <  2:59129/233‑1 (MQ=255)
                                                                                                                                                                                                                          aaaaaaacaaacaCTCCCGTTGGTTTCCTCTTCCCTTCCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAg                 <  2:418190/246‑1 (MQ=255)
                                                                                                                                                                                                                           aaaaaaaaaaccCTACCGTTGGGTTCCTCTTACCTTTCACCGGCGCTTTGAAGCTATCACCCGCTTGAAAAGTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTCCCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGGGAACGACGATCGTCACCGTAGGCTTGCGCGCCTGCCTGCAGt                <  1:408619/244‑1 (MQ=255)
                                                                                                                                                                                                                            aacaacaaCCACGACCTGTTGGGTGCTCTGCCCTTTCTCCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATTTCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCTTGCTCGCTCATCGCTTTCGCTTCGTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGtt               <  1:216785/249‑1 (MQ=255)
                                                                                                                                                                                                                             ccaacaaACACTACCGGTTGGTTGCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGGCGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttc              <  1:429142/250‑1 (MQ=255)
                                                                                                                                                                                                                             acccaaaccacccccGGTTGGTTGCTTTTCCCTTTCCCCGCCGCTTTGAAGCTATCACCCGCTGTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttc              <  2:135667/238‑1 (MQ=255)
                                                                                                                                                                                                                               aacaaccaCTCCCGTTTGGTTGCTCTTCCCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGTCCTTTAGCGCTGCGTACCCAGCCCATCTGCGGCAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttctt            <  1:254162/251‑1 (MQ=255)
                                                                                                                                                                                                                               aacaaacaCTACCGTTTGGTTGCTCTTACCTTCCCCCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCAGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttctt            <  1:47316/251‑1 (MQ=255)
                                                                                                                                                                                                                                acaaacaCTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttcttt           >  2:190087/1‑251 (MQ=255)
                                                                                                                                                                                                                                acaaacaCTACCGGTTGGTTGATCTTACCTTTAACCGCCGCTTTGAAGATACCACCCGCTTTATAATTCCGGCCCGGAGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCCGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAAAGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttcttt           >  2:481074/1‑251 (MQ=255)
                                                                                                                                                                                                                                 caaacaCTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttctttc          >  2:45357/1‑251 (MQ=255)
                                                                                                                                                                                                                                 caaacaCTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttc              >  1:135667/1‑247 (MQ=255)
                                                                                                                                                                                                                                  aaacaCTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGttctttct         >  1:467323/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aaccaccccggtggttttctctTCCCTTCCCCCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGTtctttctt        <  2:129481/233‑1 (MQ=255)
                                                                                                                                                                                                                                         cccGGTTGTTTGCTCTCACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGTTCTTTCTTcagcag  <  1:289093/250‑1 (MQ=255)
                                                                                                                                                                                                                                          ccGGTTGGTTGCTCTTACTTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATc                                               <  1:46351/205‑1 (MQ=255)
                                                                                                                                                                                                                                          ccGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATc                                               >  2:46351/1‑205 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                          
AAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAGACGAAACCGTAACCCGCATCGGAAGCCATCAGCAGTTTCTGATCGTCGCTTTCCATCAGCATATGGTCAACGGTCGCCCCAGGCGGCAACGTTAATTTGCCGGTGAGCGGCTCGCCCTGACCACGCGCCGACGGCAGCGTAATCGGGTCAATGGCATAGCTACGACCGGTGGAATCAACAAACACTACCGGTTGGTTGCTCTTACCTTTCACCGCCGCTTTGAAGCTATCACCCGCTTTATAATTCAGGCCCGGCGCGTCGATATCATGGCCTTTAGCGCTGCGTACCCAGCCCATCTGCGACAGCACAATGGTGACAGGTTCAGACGGCAGCATGTCGTGCTCGCTCATCGCTTTCGCTTCTTCGCGTTCCTGCAACGGCGAACGACGATCGTCACCGTAGGCTTGCGCGTCTGCCTGCAGTTCTTTCTTCAGCAG  >  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: