Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,335,7170CT22.6% 117.1 / 10.4 62intergenic (+50/‑114)ymiC/acnAprotein YmiC/aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base C (17/31);  new base T (14/0);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.82e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCTTATTACCCATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGT  >  NC_000913/1335469‑1335961
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       
tCTTAGTACCCATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTGAAGGGTGTAGTCTATGATCAACACAAATATGAAATATGGGTCCTGGATGGGCGCTTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTCCTTTGCCACAAGTTTtct                                                                                                                                                                                                                                                        <  2:448206/250‑1 (MQ=255)
         ccATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTTCTTTGCCACAAGGTTTTTCCTCTTTTAt                                                                                                                                                                                                                                              >  1:312165/1‑251 (MQ=255)
              tttAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAAt                                                                                                                                                                                                                                          <  1:197424/250‑1 (MQ=255)
               ttAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTTGCATTTGCTTTGCCACAAGGGTTCTCCTCTTTTATCCAttt                                                                                                                                                                                                                                        >  1:479879/1‑251 (MQ=255)
                tAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGGCATTTCGTTTTCCTCAAGGTTTCTCCTCTTTTATCAATTTg                                                                                                                                                                                                                                       >  1:504263/1‑251 (MQ=255)
                  attaatttaatggCTTTTAATACACCGCAGCAAAAACAGCTTGAGTTATATCAACACAAAAGAATAAACGTTAAGGTTGTAGCCTATGATCAACACAAATAGGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGCCGATGCTCTTCTGGTACGATCTCCGCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTTTCATTTGCTTGGCCACAATGTTTCTCATCTTTTATCAATTTggg                                                                                                                                                                                                                                     <  2:314189/240‑1 (MQ=255)
                    gAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTTGATCATTACTCACTGATCCTTGGCCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGgtt                                                                                                                                                                                                                                   <  1:504562/251‑1 (MQ=255)
                          aaGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATc                                                                                                                                                                                                                            <  1:67590/251‑1 (MQ=255)
                                  ttAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTTCCACAAGGGTTCTCCTCTTTTATCAATTTGGGGTGTTATCAAATCGtt                                                                                                                                                                                                                    >  1:52939/1‑251 (MQ=255)
                                                 aaTAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCtttt                                                                                                                                                                                                                                                <  2:344943/209‑1 (MQ=255)
                                                 aaTAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAAGGTTCTCCTCtttt                                                                                                                                                                                                                                                >  1:344943/1‑209 (MQ=255)
                                                  aTAACAGCTTGAGTTATCTTAAAACAAAATAATAAACGTTAACGGTGTAGCCTATGATCAACACAAATAGGAAATATTTGTCCTGGATGGGCGCGTTTTCTCTGTCGATTCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTGTCTCCTCTTTTATCAATTTGGGTTGGTATCAAATCGTTACGCTATGTTTGTGtt                                                                                                                                                                                                    <  2:525308/251‑1 (MQ=255)
                                                   tAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTTCCACAAGGTTTTTCCTCTTTTATCAATTTTGGTTGTTTTCAAATCGTTTCGCGATGTTTGTGttt                                                                                                                                                                                                   >  1:330613/1‑250 (MQ=255)
                                                     aCAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATc                                                                                                                                                                                                 <  1:534523/251‑1 (MQ=255)
                                                     aCAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATc                                                                                                                                                                                                 <  2:330613/251‑1 (MQ=255)
                                                     aCAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATc                                                                                                                                                                                                 <  1:401177/251‑1 (MQ=255)
                                                     aCAGCGTGAGTTATCACAACACAAAATAATAACCGTGAAGGGTGTAGCCTATGATCAACACAAATATTAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGTGCCGAACTCCTCTGGATCATTACTCACTGAGCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGTTTGTTATCAAAGCGTTACGCGATGGTTGTGTTATc                                                                                                                                                                                                 <  2:145225/251‑1 (MQ=255)
                                                      cAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGCGTCATTTGTTTTTCCCCAAGGGTTTTCCTCTTTTATCACTTTGGGTTGTTTTCAAATCGTTACGCGATTTTTTTGTTATCt                                                                                                                                                                                                >  2:424114/1‑251 (MQ=255)
                                                      aatcTTGAGTTATATCAAAACAAAATAATAACCGTTTAGGGTGTAGCCTATGATCAACACAAATATAAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTGCTGGGCCGAACTCCTCGGGATCATTACTCACTGATCCTTGTCCCCGCTGCGGCGGTGTTGTCATTTGCTTTGCCACAAGGCTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCt                                                                                                                                                                                                <  2:336320/248‑1 (MQ=255)
                                                                       acacAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTTATTTTCTTTTTCCCAAGGGTTTTCCTCTTTTATCAATTTGTGTTGTTTTCACAACGTTACGCGNTGTTTGTGTTTTTTTTAattttttacctgt                                                                                                                                                                                >  2:257064/1‑239 (MQ=255)
                                                                              ataataACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTTTCCTCTTTTTTCAATTTGGGTTGTTATCAAATCGTTTTGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATc                                                                                                                                                                       >  1:240790/1‑251 (MQ=255)
                                                                                 ataaCCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCTCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTGGAATATTCACCCTGAAGAGAATCAgg                                                                                                                                                                    <  1:167207/251‑1 (MQ=255)
                                                                                 ataaCCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAgg                                                                                                                                                                    <  2:479879/251‑1 (MQ=255)
                                                                                       gTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGGTGTTTTCAAAACGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCg                                                                                                                                                             >  1:344388/1‑251 (MQ=255)
                                                                                               tgtAGCCTATGATCAACACAAATATGAAGTATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGACCTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCCTTTTCTTTTCCACAACGGTTCTTCTCTTTTTTCAATTTTGGCTGTTTTTAAAACGCGACTCAGTGTTTGCTTTTTCTTTAATATTTACCCTGAAGAGAAACCAGGGTTCGCCACCCTg                                                                                                                                                     >  2:312379/1‑251 (MQ=255)
                                                                                                  aGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCCTTTGCTTTGCCCCAAGGTTTCTCCTCTTTTTTCAATTTTGGGTGTTATCAAATCGGTACGCGATGTTTTTGTTATCTTTAATATTTACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTcc                                                                                                                                                  >  1:161163/1‑250 (MQ=255)
                                                                                                       aTGATCAACACAACTATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGGAATTTGTGTTGCCACAAAGTTTTTCCTCTTTTATCAATTTGGGTTGTTTTTCAAACGTTACGCGGTGTTTTTTTTTTCTTTTATATTTACCCTGAAGAGAATCAGGGATTTGCCACCCTGTCAttta                                                                                                                                              >  2:78194/1‑248 (MQ=255)
                                                                                                       aTGATCAACACAAATATGATATATTTGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTTTTCTGGGCCGAACTCCTCTGGATCATTACTCACTTATCCTTGACCCCTCTGCGGCGGGGTTGTCATTTGCTTTGCCACAATGTTTCTCCTCTTTTATCAATTTTGGTTGGTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCAttt                                                                                                                                               <  2:68412/249‑2 (MQ=255)
                                                                                                          aTCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACACAGTTTCTCCTCTTTTATCAATTTGGGTTGTTATTAAAACGTTTCGCGATGTTTGTGTTTTCTTTAATATTCACCCTGAAGAGAATTAGGGCTTCGCAACCCTGTCATTAagg                                                                                                                                            >  1:499739/1‑249 (MQ=255)
                                                                                                           tCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACCAGGTTTCTCCTCCTTTTTCAATTTTGGTTGTTTTCAAATCGTTACGCCGTGTTTGTGTTATTTTTTATATTTACCCTGAAGAGAAATCGGGCTTCGCAACCCTGTCATTattgagg                                                                                                                                         >  2:501234/1‑245 (MQ=255)
                                                                                                            cAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACGCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAaggagg                                                                                                                                         <  1:47246/250‑1 (MQ=255)
                                                                                                               cacaAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCATGGCTTCGCAACCCTGTCATTAAGGAGGAGCt                                                                                                                                     <  1:199021/251‑1 (MQ=255)
                                                                                                               cacaAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTTATTTGCTTTTCCACAAGGTTTTTCCTCTTTTATCAATTTTGGTTGTTTTCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAAGGAATCAGGGCTTTCCCAACCTGTCATTAAGGAGGAGCt                                                                                                                                     >  1:264057/1‑251 (MQ=255)
                                                                                                                   aaTATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTTCTTTTCCACAAGGTTTCTCCTCTTTTTTCAATTTTGGGTGTTATTAAATCTGTACGCGGTGTTTGTGTTATCTTTAATATTTACCCTGTAGAGAATCAGGGCTTCGCAACCCTTTCATTAAGGAGGAGCTATgt                                                                                                                                 >  1:398930/1‑251 (MQ=255)
                                                                                                                      atGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATgtcgt                                                                                                                              <  1:4304/251‑1 (MQ=255)
                                                                                                                      atGAAATATTGGTACTGGATTGGCGCGTTTTCTCTGTCGATGCTCTTCTTGGCCGAACTCCTCTGTATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTTCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATgtcgt                                                                                                                              <  1:270244/251‑1 (MQ=255)
                                                                                                                           atatTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTTCCACAAGGGTTTTTCCCTTTTATCAATTTGGGGTGTTATCCAATCGTTACGCGATGTTTGTGTTTTCTTTTATATTTACCCTGAAGGGAATCAGGGCGTCGCAACCCCGTCCTTACAGAGGCGCAATGTCGTCaaac                                                                                                                         >  1:387728/1‑249 (MQ=255)
                                                                                                                           atatTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCCCCAGGGTTTTCCTCTTTTTTCAATTTGGGGTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTTACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAAcc                                                                                                                         >  1:3898/1‑251 (MQ=255)
                                                                                                                            tatTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCCTTTGCTTTTCCACAAGGTTTTTCCTCTTTTTTCAATTTGGGTTGGTATCAAAACGTTACGCGATGTTTGTGTTTTTTTTAATATTCACCCCGTAGAGAATCAGGGGTTTGCAACCCTTTCATTTAGGAGGAGCTATGTCGTCAAccc                                                                                                                        >  2:526164/1‑251 (MQ=255)
                                                                                                                                ggTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTTCCACAAAGTTTTTCCTCTTTTATCAATTTGGGGTGTTATCAAATCGTTACGCGATGTTTGTGTTTTTTTTAATATTCAACCTGAAAGGAATCAGGGCTTCCCCACCCTGTCATTTAGGGGGGGGTTTGTCGTCaaacccacg                                                                                                                    >  1:308670/1‑244 (MQ=255)
                                                                                                                                      ggATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAg                                                                                                               <  1:526164/250‑1 (MQ=255)
                                                                                                                                          gggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAAt                                                                                                                                                                                                                                          <  2:277217/126‑1 (MQ=255)
                                                                                                                                          gggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCCTTTTCTTTGCCACAAGGGTTCTCCTCTTTTATCAAt                                                                                                                                                                                                                                          >  1:277217/1‑126 (MQ=255)
                                                                                                                                                gTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTTTCATTTGCCTTTCCACCAAGTTTCTCCTCTTTTATCAATTTGGGTGTTTATCAAATCGTTACGCGCTGTTTGTGTTATCTTTAATATTTACCCTGACGAGAATCAGGGCTTCGCACCCCCGGCATTAAGGAGGAGCTATGTCGGCAAccccccgaagagccagcaaggcc                                                                                                    >  2:27988/1‑231 (MQ=255)
                                                                                                                                                gTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTTCCACCAGGTTTCTCCTCTTTTATCAAATTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACAAGAAGCCAGTCAGGa                                                                                                     >  1:315001/1‑250 (MQ=255)
                                                                                                                                                         gTCGATGCTCTTCTGGGCCGAACTCCTCTGGAACATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTGGGGTTGTTATCAAATCGTTACGCGATGGTTGTGTTATCTGTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCTTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAg                                                                                           <  1:193400/251‑1 (MQ=255)
                                                                                                                                                            gATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAACTTGGGGTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCCGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAAGACAACTTGCAGGcc                                                                                        >  1:76434/1‑251 (MQ=255)
                                                                                                                                                               gCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCGTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTAGCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGGGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCaaa                                                                                     <  1:78194/251‑1 (MQ=255)
                                                                                                                                                                 tcttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTa                                                                                                                                                                                                   <  1:485543/141‑1 (MQ=255)
                                                                                                                                                                 tcttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCCTTTTCTTTTCCACCAGGTTTTTCCTCTTTTTTTAATTTTGGGTGTTATTAAATCGTTACGCGGTGTTTTTGttt                                                                                                                                                                                                   >  2:485543/1‑140 (MQ=255)
                                                                                                                                                                  cttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAg                                                                                    <  1:225923/249‑1 (MQ=255)
                                                                                                                                                                    tctGGGCCGAACTCCTCTGGATCATTACCCACTTATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTAACAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATa                                                                                 <  1:363974/250‑1 (MQ=255)
                                                                                                                                                                            gAACTCCTCTGGATCATTACTCACTGATCCTNGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTAcc                                                                        <  2:264057/251‑1 (MQ=255)
                                                                                                                                                                                      ggATCATTACTCACTGGTCCTTGACCCCNCTGCGGCTGGGTTGTCAGTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGGTTGTGTTATCGTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCActacta                                                                 <  2:266829/248‑1 (MQ=255)
                                                                                                                                                                                       gATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCCTTTGCTTTTCCACAAGGTTTTTCCTCTTTTTTAAATTTTGGATGTTTTCAAATCGTTACGCCATGATTTTGTTATTTTTTATATTCTCCCCGAAGGGAATCAGGGCATCGCAACCCTGTCATTAAGGAGGAGCTATTTCGTCAAACCTAAGTGAAGCCAGGAAGGGCACGATGCAGGACAACGATAAAACTTACCACTACTAAAGc                                                             >  1:405173/1‑251 (MQ=255)
                                                                                                                                                                                        aTCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGGTTTTCCTCTTTTATCAATTTGGGGTGTTTTCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTTCCCCTGAAGAGAATCAGGGCTTTGCAACCCTGCCATTTAGGAGGAGCTATGTAGTCAAACCTACGAGAAGCCAGTAAGGACAAGTTGCAGGCCCAAGATAAAACTTAACACCACTACAGcc                                                            >  1:101903/1‑251 (MQ=255)
                                                                                                                                                                                                 cACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGGTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCGTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACGTACCACTACTACAGCCTGCCGCTtg                                                   <  1:173857/251‑1 (MQ=255)
                                                                                                                                                                                                     gATCCTTGCCCCCGCTGCGGCGGGGTTGTCATTTGTTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGATATGTCGTCAACCCTACGAGAAGCCAGTAAGGATACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTtgct                                                 <  1:484085/249‑1 (MQ=255)
                                                                                                                                                                                                     gATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCGGCCGCTtgctgc                                               <  1:110560/251‑1 (MQ=255)
                                                                                                                                                                                                     gATCCTTGACCCCGCTGCGGCGGGGGTTTCCTTTTCTTTGCCCCAAGGTTTCTCCTCTTTTTTCAATTTTGGTTGTTATCAAAATGTTACGCGCTGTTTGTGTTTTCTTTAATATTCTCCCCGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACACTACGAGAAACCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAACCTGCCGCTtgctgc                                               >  1:383213/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCt                    >  1:450917/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tCATTTGCTTTGCCACAAGGTTTCTCCTCTNTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCt                    <  2:240790/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   aGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGCCACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGt  >  2:497131/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       
TCTTATTACCCATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCTACCCAAGTCACTCAAAGT  >  NC_000913/1335469‑1335961

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