Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,335,7050GT25.4% 102.4 / 10.6 63intergenic (+38/‑126)ymiC/acnAprotein YmiC/aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base G (15/32);  new base T (16/0);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 8.20e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCATCATTCTTATTACCCATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCT  >  NC_000913/1335462‑1335943
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
tcatccTTCATATTAACCATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAATAGCTTGAGTGATCTCAACACAAAATAATAACCGTTAAGGGTGCCGCCTATGATCCACACAAATATGAGATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGTACGCCTCTGTACCATTACTCACGAATCCTTGCCACCGCGGCGGCGGCGTTGAGATTTGCATTGGcacctg                                                                                                                                                                                                                                             >  2:542802/1‑248 (MQ=255)
       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                                                                                                                                                                               <  2:330613/251‑1 (MQ=255)
                                                            aCAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATc                                                                                                                                                                               <  1:401177/251‑1 (MQ=255)
                                                            aCAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATc                                                                                                                                                                               <  1:534523/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)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
TCATCATTCTTATTACCCATTTTTAATGAATTAAAGGGCTTTTAATACACCGCAGCAATAACAGCTTGAGTTATCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTATGTCGTCAACCCTACGAGAAGCCAGTAAGGACACGTTGCAGGCCAAAGATAAAACTTACCACTACTACAGCCTGCCGCTTGCTGCTAAATCACTGGGCGATATCACCCGTCT  >  NC_000913/1335462‑1335943

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