Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,905,761	G→A	G25D (GGT→GAT)	mntP →	Mn(2(+)) exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,905,761	0	G	A	100.0%	159.9 / NA	46	G25D (GGT→GAT)	mntP	Mn(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/24);  total (22/24)

TGCAAAACGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATGAAGATGAAGAGCCGCGCCGTCGACACGGTTTCTGGC  >  NC_000913/1905520‑1906006                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      tGCAAAACGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGCATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCt                                                                                                                                                                                                                                              >  2:499433/1‑251 (MQ=255)         gTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAAc                                                                                                                                                                                                                                      <  2:130198/251‑1 (MQ=255)                gTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGGCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGAGTATGTGTTTTAGGGATACCCCGGTCAGGACATTGTCATGAATATCGCGGCTACTGTTCTTCTTGCGTTGGGTATGTCGATGGATGCATTTGCTGCATCAATCGTTAAATATGCCATCCTCCATGAACCTAAAtt                                                                                                                                                                                                                               <  2:283486/251‑1 (MQ=255)                   cGGACTGCATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCCTGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGGTGCATTTTCTGCATAAATCGGTAAAGATGCCACCCTCCATAAACCGCAAttttt                                                                                                                                                                                                                            >  2:374040/1‑250 (MQ=255)                     gACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTg                                                                                                                                                                                                                          >  2:410817/1‑251 (MQ=255)                           aTACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAt                                                                                                                                                                                                                    >  1:160471/1‑251 (MQ=255)                            tACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                                                                                                                                                                                                                   >  2:266880/1‑251 (MQ=255)                                   agtcagGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAcc                                                                                                                                                                                                            <  1:73207/251‑1 (MQ=255)                                        ggCGAGACCATATGCACATCAATCGCTATGCCTGCATGCACGCGAACACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATAAACCGGTCAGGCCATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGCTGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAAACGGCCt                                                                                                                                                                                                       >  1:550121/1‑251 (MQ=255)                                             gaCCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTTCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTTAAGCATTGCGAACCGGCCTTAttt                                                                                                                                                                                                  <  1:110795/251‑1 (MQ=255)                                              aCCATATGCACATCAATCGCTATGCCTGCATGAACGCGCAAACCATCACGGCTTGTTGTTCATATCTGCATGTTTGCTGGGGACGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGACTATAACTGCTACTGTTCTTCTTGCGTTTGGTATGAAGATGGAAGCATTTGCTGCATCACTCAGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAtttt                                                                                                                                                                                                 >  1:256177/1‑251 (MQ=255)                                              aCCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGGTCATGAATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGACATGAATATCACGGCTACTGTTCTTCTTGCGTTTGGTATGGCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACACTCCATAAACCGAAATTTTATGAAGCATTGCGAACCGGCCTTAtttt                                                                                                                                                                                                 >  2:475392/1‑251 (MQ=255)                                                   atGCACATCAATCGCTATTCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGTGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGt                                                                                                                                                                                             <  2:342809/250‑1 (MQ=255)                                                            aaTCTCTATGCCTGCATGACCGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGTCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTTCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGaaa                                                                                                                                                                                   <  2:424674/251‑1 (MQ=255)                                                             aTCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTGCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGGATGTCGATGGATGCATTTGCTGCATCAAGCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGGGCCGGCGAAAc                                                                                                                                                                                  <  2:256492/251‑1 (MQ=255)                                                                           aTGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTg                                                                                                                                                                     >  1:143492/1‑250 (MQ=255)                                                                                   gAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTgg                                                                                                                                                            >  1:516442/1‑251 (MQ=255)                                                                                      aaccagaaCGGCTTGTTGTTCATGAATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTGTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCGGACGCCGCGGATCGGCTGGGGa                                                                                                                                                         <  2:529807/244‑1 (MQ=255)                                                                                                   tgtttgtCATGTATGCAGTTTTGCTGGTGGCGATGATGTGTGTTATGGAGACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCTTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATTCCACCCGCCATAAACCGAAATTTTCTGAAGCATTGCGAGCCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTTGGGAATGGGCATGtt                                                                                                                                              <  2:350130/243‑1 (MQ=255)                                                                                                                     gTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCtt                                                                                                                           <  2:458351/250‑1 (MQ=255)                                                                                                                         gCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATg                                                                                                                      <  2:454744/251‑1 (MQ=255)                                                                                                                                 ggatgatGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACt                                                                                                              <  2:160471/250‑1 (MQ=255)                                                                                                                                          gtTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCTAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGt                                                                                                     <  1:288037/251‑1 (MQ=255)                                                                                                                                                tGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGCTGGATGAATTAGCTGCACCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATATTCAGAAGCATTGCGAACCGGCCTTCTATTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGACTGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGtgc                                                                                               >  2:512078/1‑251 (MQ=255)                                                                                                                                                 ggATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGtgct                                                                                              >  1:69590/1‑251 (MQ=255)                                                                                                                                                   aTACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGCGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGtgctgc                                                                                            <  1:410817/251‑1 (MQ=255)                                                                                                                                                     aCCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGGGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGctgctg                                                                                          <  2:516442/251‑1 (MQ=255)                                                                                                                                                     aCCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCGTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCGGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTgctgct                                                                                           <  2:389535/250‑1 (MQ=255)                                                                                                                                                         cGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGCTTGCGTTTGTGCTGCTGatat                                                                                      >  2:559329/1‑251 (MQ=255)                                                                                                                                                          ggTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATAtt                                                                                     >  1:549158/1‑251 (MQ=255)                                                                                                                                                           gTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTTCTGCATCAATCGTTAAAGATGCCACCCGCCATAAACCGAAATTTTCTGAAGCATTGCGAGCCGGCCTTATTTTTGGGGCCGTCGAAACCCTTACGCCGCTGATCTGCGGGGGAAGTGGCAGGTTAGCCAGCGGGGTGGTCCTTGAATGGAACCACTGGATTGCGGTTGTGCTGCTGATATTc                                                                                    <  1:244200/251‑1 (MQ=255)                                                                                                                                                               ggACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTcg                                                                                <  1:10515/251‑1 (MQ=255)                                                                                                                                                                   aTTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTcggcgg                                                                            >  1:359371/1‑251 (MQ=255)                                                                                                                                                                               atCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTACAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGGATGattat                                                                >  1:229492/1‑251 (MQ=255)                                                                                                                                                                                 cACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCGTCAATCGGTAACGATGCCCCCCTCCATAAAACGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCggcggggggcagaattttt                                                              >  2:405394/1‑238 (MQ=255)                                                                                                                                                                                             tcttctGGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGt                                                   <  2:339290/250‑1 (MQ=255)                                                                                                                                                                                                     cGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCTGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAgat                                          <  1:374040/251‑1 (MQ=255)                                                                                                                                                                                                      gTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAgatg                                         <  1:499433/251‑1 (MQ=255)                                                                                                                                                                                                          ggTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGAtgatga                                     <  2:549158/251‑1 (MQ=255)                                                                                                                                                                                                                  gatggatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGGGAATGATTATTGAGGGGTTTTGTGGCGCAGATgatgaagatga                               >  1:495742/1‑249 (MQ=255)                                                                                                                                                                                                                      gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGGTTTCGTGGCGCGGATGATGAAGATGAAGAGcc                         >  1:507232/1‑251 (MQ=255)                                                                                                                                                                                                                          cATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGattatt                                                               <  1:337564/209‑1 (MQ=255)                                                                                                                                                                                                                                  gcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGCACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCggcggg                                                                           >  2:294844/1‑189 (MQ=255)                                                                                                                                                                                                                                  gcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCggcggg                                                                           <  1:294844/189‑1 (MQ=255)                                                                                                                                                                                                                                   cATCAATCGGTAAAGATGCCACCCTCCATCAACCGAACTTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGGGCATGGGTTTTGAGGGTTTTCGTGGCGCAGATTATGCAGGTGACGAGCCGCGCCgggcctcc            >  2:586539/1‑244 (MQ=255)                                                                                                                                                                                                                                             tAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGGTTATTGGGGGGTTTCGTGGCGCAGATGATGCAGATGAAGAGCCGCGCCGTCGACACGGTTGCTGGc  >  1:11852/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      TGCAAAACGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATGAAGATGAAGAGCCGCGCCGTCGACACGGTTTCTGGC  >  NC_000913/1905520‑1906006

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