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%	215.1 / NA	58	G25D (GGT→GAT)	mntP	Mn(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/30);  total (28/30)

AAAACGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATGAAGATGAAGAGCCGCGCCGTCGACACGGTTTC  >  NC_000913/1905523‑1906002                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  aaaaCGTGGCACTTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCAATGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                                                                                                                                                                                                        >  2:139048/1‑250 (MQ=255)        ggCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAAc                                                                                                                                                                                                                                  >  2:316178/1‑249 (MQ=255)        ggCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAAc                                                                                                                                                                                                                                  <  1:316178/249‑1 (MQ=255)        ggCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCg                                                                                                                                                                                                                                >  1:178514/1‑251 (MQ=255)                          aCTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                                                                                                                                                                                                               <  2:178514/250‑1 (MQ=255)                          aCTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCAtt                                                                                                                                                                                                               <  1:571287/250‑1 (MQ=255)                                   cagGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGc                                                                                                                                                                                                     <  2:286818/251‑1 (MQ=255)                                   cagGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGCCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGCAATTTTCTGAACCATTGCGAACCGGc                                                                                                                                                                                                     >  2:233612/1‑251 (MQ=255)                                     ggCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCt                                                                                                                                                                                                   <  1:30772/251‑1 (MQ=255)                                          gaCCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGCTGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCCCCCTCCCTAAACCGAAATTTTCTGACGCATTGCGAACCGGCCTTAttt                                                                                                                                                                                              >  2:399310/1‑251 (MQ=255)                                              atatGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTACTTTTgg                                                                                                                                                                                          <  1:394121/251‑1 (MQ=255)                                                     caTCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTGCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAGTTTTGGTGCCGTc                                                                                                                                                                                   <  2:206754/251‑1 (MQ=255)                                                                aTGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCTGTCGAAACCCTg                                                                                                                                                                          <  1:420350/249‑1 (MQ=255)                                                                    cTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCCTGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTCTGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCTTGGTAACCCTGccgccg                                                                                                                                                                    >  2:80469/1‑251 (MQ=255)                                                                      gCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCt                                                                                                                                                                  >  2:367981/1‑251 (MQ=255)                                                                              gcgAAAGCCATCACGGATTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCt                                                                                                                                                          <  2:341643/251‑1 (MQ=255)                                                                               cgAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCAGCTg                                                                                                                                                         <  1:513198/251‑1 (MQ=255)                                                                                    aCCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGGGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGaa                                                                                                                                                    <  1:233612/251‑1 (MQ=255)                                                                                       aTCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATgg                                                                                                                                                 <  1:133390/251‑1 (MQ=255)                                                                                         cACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGCGCGATGATGTGTTTTATGGATACCCCGGTCAGTACATTGTCATGAATATCACTGCTAATGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCTAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAATCCCTGACGCCGCTGATCGGCTGGGGAATggg                                                                                                                                                <  2:302078/250‑1 (MQ=255)                                                                                                tgttgtTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTg                                                                                                                                                                 >  2:209597/1‑226 (MQ=255)                                                                                                tgttgtTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTg                                                                                                                                                                 <  1:209597/226‑1 (MQ=255)                                                                                                   tgttCATGTTTCCATGTTAGCTGGGGGTGAGGATGTTTTTTATGGATACCCCGGTCATGACATTGTCATGAATATCACTGCGACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGATCCGGCCTTAGTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATTGGCATGTTatct                                                                                                                                      <  2:329236/250‑4 (MQ=255)                                                                                                   tgtacATGTATGCATGTTTGCTGGTGGAGATGATGTGTTTTATGGATACGCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTGGCGGTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCCCCCTCCATAAACCGAAATTTTCTGAAGCAGTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTagcc                                                                                                                                      <  2:301115/246‑1 (MQ=255)                                                                                                    gttCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTagccag                                                                                                                                    <  1:367981/251‑1 (MQ=255)                                                                                                            aTGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCt                                                                                                                                                          <  1:194718/221‑1 (MQ=255)                                                                                                            aTGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCt                                                                                                                                                          >  2:194718/1‑221 (MQ=255)                                                                                                            aTGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCGTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGttt                                                                                                                             <  2:196990/250‑1 (MQ=255)                                                                                                                                   ntgtgTTTTATGGATAAAAAGGTCAGGACATTGTTATGAATATTATTGTTAGTGTTATTTTTGAGTTTGGTATGTAGATGGATGCATTTGTTGCATAAATCGGTAAAGATGCCACAATCCATAAACCGAAATTTTATGAAGCATTGCGAACCGGCCTTATTTTTGGTGCAGTCGAAACTCTGATGCCGATGATAGGCTGGGGAATGGGCATGTTAGCTAGCCGGTTTGTCATTGAATGGAACTATTGGAtt                                                                                                     >  2:267017/2‑251 (MQ=255)                                                                                                                                   atgTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGAtt                                                                                                     >  1:65037/1‑251 (MQ=255)                                                                                                                                       gtTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAAc                                                                                                                                                                                                                                  >  2:97056/1‑122 (MQ=255)                                                                                                                                       gtTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAAc                                                                                                                                                                                                                                  <  1:97056/122‑1 (MQ=255)                                                                                                                                        ttttATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAttt                                                                                                                                                                                              <  2:297632/157‑1 (MQ=255)                                                                                                                                        ttttATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAttt                                                                                                                                                                                              >  1:297632/1‑157 (MQ=255)                                                                                                                                           tATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTtgt                                                                                             >  1:73319/1‑251 (MQ=255)                                                                                                                                            aTGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTtgtg                                                                                            >  1:164846/1‑251 (MQ=255)                                                                                                                                            aTGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTtgtg                                                                                            <  1:407344/251‑1 (MQ=255)                                                                                                                                                  aCCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGctgctg                                                                                      <  1:1179/251‑1 (MQ=255)                                                                                                                                                          cAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCt                                                                              >  2:37516/1‑251 (MQ=255)                                                                                                                                                                     cATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCggcgggcgc                                                                    >  1:92208/1‑249 (MQ=255)                                                                                                                                                                      aTGAATATCCCTGCTCCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCGCCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAAt                                                                  <  2:164846/251‑1 (MQ=255)                                                                                                                                                                         aaTATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGaacggt                                                               >  2:92136/1‑247 (MQ=255)                                                                                                                                                                              cACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGattatt                                                           <  1:534775/250‑1 (MQ=255)                                                                                                                                                                              cACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTg                                                          >  1:498316/1‑251 (MQ=255)                                                                                                                                                                               aCTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCGTCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGattattt                                                          >  1:221859/1‑249 (MQ=255)                                                                                                                                                                                     aCTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGttt                                                   <  1:430314/251‑1 (MQ=255)                                                                                                                                                                                     aCTGTTCTTCTTACGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCAGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGAGAATGATTATTGAGGGttt                                                   >  1:408485/1‑251 (MQ=255)                                                                                                                                                                                      cTGTTCTTCTTGCGTTTGGTATGTCGTTGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAAttgttatttcggggttt                                                  >  2:98319/1‑235 (MQ=255)                                                                                                                                                                                        gTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGGTTTCg                                                >  1:415921/1‑251 (MQ=255)                                                                                                                                                                                                   gTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAAt                                                                                                                                                   >  2:212017/1‑141 (MQ=255)                                                                                                                                                                                                   gTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAAt                                                                                                                                                   <  1:212017/141‑1 (MQ=255)                                                                                                                                                                                                           tGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTTGTGGCGCAGAtgatga                                 >  1:479363/1‑247 (MQ=255)                                                                                                                                                                                                            gTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGGATGATTATTTAGGGGTTTCGGGGGGCAGATgatgaagatg                            >  1:250339/1‑251 (MQ=255)                                                                                                                                                                                                                   gatgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATGTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAGCCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATGAAGATGAAGAGcc                     <  2:43120/251‑1 (MQ=255)                                                                                                                                                                                                                    atgCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGGATGATTATTGAGGGGTTTCGTGGCGCAGATGATGAAGATGAAGAGCcg                    >  1:522654/1‑251 (MQ=255)                                                                                                                                                                                                                           ggctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATGAAGATGAAGAGCCGCGCCGTc             <  2:56383/250‑1 (MQ=255)                                                                                                                                                                                                                                     aTCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTGTTGAGGGGTTTCGGTGCGCAGATGGTGAAGATGAAGAGCCGCGCCGTCGACACgggtt   >  1:68358/1‑248 (MQ=255)                                                                                                                                                                                                                                       cGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATGAAGATGAAGAGCCGCGCCGTCGACACGGTTTc  <  1:77760/250‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  AAAACGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATGAAGATGAAGAGCCGCGCCGTCGACACGGTTTC  >  NC_000913/1905523‑1906002

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