Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	225,875	G→A	noncoding (117/2904 nt)	rrlH →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	225,875	0	G	A	100.0%	117.3 / NA	57	noncoding (117/2904 nt)	rrlH	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/30);  total (27/30)

CTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCT  >  NC_000913/225628‑226116                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  cTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAAt                                                                                                                                                                                                                                                >  2:732381/1‑251 (MQ=21)   gCGAAGTTTTGCGCTTTAAACATCTGGATCAAGCTGAAAATTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGACCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGTCGATTTCCAAATgg                                                                                                                                                                                                                                              <  2:499190/251‑1 (MQ=255)           ttGCTCTTTCAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAAccc                                                                                                                                                                                                                                      >  1:216698/1‑251 (MQ=21)                  ttAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAgtg                                                                                                                                                                                                                                  <  1:461073/248‑1 (MQ=25)                   tAAAAATCTGGATCAAGCTGACAATTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGAAGTCAGAGGCGATGAAGGACGTGCTACTCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGGGtt                                                                                                                                                                                                                              >  2:619638/1‑251 (MQ=255)                      aaaTCTGGATCAAGCTGAAAATTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAGCCCAGTGTGTTTCg                                                                                                                                                                                                                           <  1:122774/251‑1 (MQ=21)                       aaTCTGGATCAAGCAGAAAATTGAAACTCTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATGTGTGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGATGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAAGGGGGAAACCCAGTGTGTTTCGa                                                                                                                                                                                                                          <  2:398354/251‑1 (MQ=11)                         tCTGGATCAAGCTGAAAATTGACACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGac                                                                                                                                                                                                                         >  1:196812/1‑250 (MQ=25)                         tCTGGATCAAGCTGAAAATTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGac                                                                                                                                                                                                                         >  2:78009/1‑250 (MQ=25)                                cAAGCTGAAAAGTGAAACACGGAACAATGAAAGTTATTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGCCGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATTACCTGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATc                                                                                                                                                                                                                 <  2:648543/251‑1 (MQ=11)                                 aaGCTGCAAATTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCCAGAAACATCTTCGGGATGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATACCCGGCGATTTCCAAATGGGGAAACCCAGAGTGTTTCGACACACTATc                                                                                                                                                                                                                 >  2:692297/1‑250 (MQ=255)                                     tGAAAATTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTa                                                                                                                                                                                                             >  1:151389/1‑250 (MQ=25)                                                 cacTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAg                                                                                                                                                                                                <  1:435665/251‑1 (MQ=21)                                                  acTGAACAATGAAAGTTGTTCGTGAGTCTCCCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTCAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGCCACACTATCATTACCTGAATCCATAgg                                                                                                                                                                                               >  2:701696/1‑251 (MQ=255)                                                            gAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAgg                                                                                                                                                                                      <  1:109961/250‑1 (MQ=18)                                                               aGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGaa                                                                                                                                                                                  <  2:559743/251‑1 (MQ=18)                                                                gttgttCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGTCACACTATCATTAACTGAATCCATAGGTGAATGAGGCGAAc                                                                                                                                                                                 <  2:90626/251‑1 (MQ=18)                                                                  tgttCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAAcc                                                                                                                                                                                <  1:709327/250‑1 (MQ=18)                                                                   gttCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGATAATGAGGCGAACCgg                                                                                                                                                                              >  1:335046/1‑251 (MQ=18)                                                                       gTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTACGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCG‑TAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGa                                                                                                                                                                          >  1:62770/1‑250 (MQ=16)                                                                                    aTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTa                                                                                                                                                             <  1:557924/251‑1 (MQ=18)                                                                                    aTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTa                                                                                                                                                             <  1:26093/251‑1 (MQ=18)                                                                                           cAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGACCTGGCAGTCAGAGGCGATGAAGGCCGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTAccc                                                                                                                                                      >  1:734896/1‑251 (MQ=18)                                                                                            aaCACGATGCTGGATCGCAAGAAACCTCTTCGGGTTGTGCGGTTAAGCGACTAAGCGTACACGATGGATGCCCTGGCAGTCAGAGGCGCTGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAACCATCTAAGTAccac                                                                                                                                                     >  2:426333/1‑249 (MQ=11)                                                                                            aaCACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTGCCATATGGTGAAACCCAGTGTGTTTCGACACGCTATCATTAACTGAAGCCATAGGTTAATGAGGCGAACCGGGTGATCTGGAACATCGAAGTAcccc                                                                                                                                                     <  2:737796/251‑1 (MQ=11)                                                                                                cGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGACCCGGGGGAACTGAAACATCTAAGTACCCCGAgg                                                                                                                                                 >  2:535248/1‑251 (MQ=18)                                                                                                  atgatgGATCGCAAGAATCATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGaa                                                                                                                                               <  2:439611/251‑1 (MQ=18)                                                                                                    gatgGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGaaa                                                                                                                                              >  2:440653/1‑250 (MQ=18)                                                                                                     atgGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGaaa                                                                                                                                              >  1:370915/1‑249 (MQ=18)                                                                                                                                                  gCGTACACTGTTGATGCCCTTGCAGTCAGATGTGATGAAGGACGTGCTAATCTGCGATAAGTGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGATGAAAAGAAATCAACCTAGATTCCCCCAGTAGCGGCGAGCTAACGGagagga                                                                                               <  2:116067/251‑1 (MQ=11)                                                                                                                                                     tACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATAAGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGCACTGAAAAATCTAAGTACCCCGAGGAACAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGcc                                                                                            >  1:355880/1‑251 (MQ=9)                                                                                                                                                      acacGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTACGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAAGGGGgaggaggcc                                                                                           >  2:576943/1‑248 (MQ=9)                                                                                                                                                          ggTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCagag                                                                                       <  2:196812/251‑1 (MQ=14)                                                                                                                                                              gATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTg                                                                                   >  2:198762/1‑251 (MQ=14)                                                                                                                                                                  cccTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATc                                                                               <  1:41553/251‑1 (MQ=14)                                                                                                                                                                   ccTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCa                                                                              <  1:732381/251‑1 (MQ=14)                                                                                                                                                                   ccTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCAGAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAATTACCCCGAGTAAAAGAAATCAACCGAGATTCCTCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCa                                                                              <  2:734896/251‑1 (MQ=9)                                                                                                                                                                                 gCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGt                                                                <  2:151389/251‑1 (MQ=14)                                                                                                                                                                                          gACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTc                                                       <  2:97691/251‑1 (MQ=14)                                                                                                                                                                                           aCGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCt                                                      <  1:619638/251‑1 (MQ=14)                                                                                                                                                                                            cGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTg                                                     >  2:245669/1‑251 (MQ=14)                                                                                                                                                                                             gTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGGGTGTGTGTTAGTGGCAGCGTCTgg                                                    >  2:165663/1‑251 (MQ=14)                                                                                                                                                                                                       gCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGgcgc                                            <  1:78009/249‑1 (MQ=14)                                                                                                                                                                                                        cGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGTGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGGTAAGGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGCAATCAACCGAGAGTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGATGCAGTGTGTGTGTTAGTGTAATCGTCTGGAAAGgcgcgtg                                         <  2:668199/251‑3 (MQ=11)                                                                                                                                                                                                                  cGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGAAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTg                               <  2:395421/251‑1 (MQ=14)                                                                                                                                                                                                                    gTAAGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGGGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGAc                             >  2:403210/1‑251 (MQ=14)                                                                                                                                                                                                                      aaGGTGATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAg                           <  2:63671/251‑1 (MQ=14)                                                                                                                                                                                                                          ggATATGAACCGTTATAACCGGCGATTTCCAAATGGGGAAAACCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGGACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGcccc                       <  1:44531/250‑1 (MQ=9)                                                                                                                                                                                                                            atatGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCCTTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGAGTTAGTGGAAGCGTCTGGAACGGCGCGCGATACAGGGTGACAGCCCCgg                     >  1:544376/1‑250 (MQ=9)                                                                                                                                                                                                                              atGAACCGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTac                   >  1:103043/1‑251 (MQ=14)                                                                                                                                                                                                                                    cGTTATAACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACaaaa             <  1:254458/251‑1 (MQ=14)                                                                                                                                                                                                                                    agTTATAACCGGCGATTTCCAAACGGGGAAACCCAGAGTGTTTCGACCCACTATCATTAACTGAATCCATAGGTGAATGAGGCGAACCGGGGGAACTGAAACATCTCAGTACCACGTGGAAAAGAACTCAACCGAGATTCCCCAAGTAGCGGCGAGCGAACGGAGAGGAGCCCAAAGCCTGTATCAGCGTGAGAGTTACTGGAAGCGTCTGGAAAGGCGACCGATACAAGGTGACAGCTCCgcaccaaaa              >  2:614399/2‑242 (MQ=11)                                                                                                                                                                                                                                       tataACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGCACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGCGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGAGACAGCCCCGTACACAACAATg          >  1:658723/1‑251 (MQ=9)                                                                                                                                                                                                                                         taACCGGCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGc         >  1:53985/1‑250 (MQ=14)                                                                                                                                                                                                                                           acctgCGATTTGCAAATGGGGAAACCCAGTGTGTGTCGACACACTATCATAAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGAGTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGGCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGcaca      <  2:62770/247‑1 (MQ=11)                                                                                                                                                                                                                                               gCGATTTCCAAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAAAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCt  >  1:33330/1‑251 (MQ=14)                                                                                                                                                                                                                                               gCGATTTCCAAATGGGGAAACCCAGTGTGTGTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCt  <  1:576082/251‑1 (MQ=14)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  CTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCAACCGAGATTCCCCCAGTAGCGGCGAGCGAACGGGGAGGAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCT  >  NC_000913/225628‑226116

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑