Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	402,615	G→A	W290* (TGG→TAG)	phoA →	alkaline phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	402,615	0	G	A	100.0%	151.3 / NA	45	W290* (TGG→TAG)	phoA	alkaline phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/24);  total (21/24)

GGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGATCTCGATGA  >  NC_000913/402371‑402852                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              ggAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAgg                                                                                                                                                                                                                                         <  1:476175/251‑1 (MQ=255)    aaaGGATCGATTACCCAACACCGTCTTAAACCTCGCGCCGACGTTACGTTTGTCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCTCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACTGGCGCGTGTTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGGTACGGAAGCGAATCAGCAAATCCCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGAcc                                                                                                                                                                                                                                      <  2:229322/251‑1 (MQ=255)     aaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCg                                                                                                                                                                                                                                     <  2:178711/251‑1 (MQ=255)      aGGATCGATTACCGAACAGCTGCTTAACGCTCGTGACGACGTTACGCTTGGCGGAGGCGCAAAAACATTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAACACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCg                                                                                                                                                                                                                                     >  1:226392/1‑250 (MQ=255)               tACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACg                                                                                                                                                                                                                            <  1:88329/250‑1 (MQ=255)                   gAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTAcc                                                                                                                                                                                                                        >  1:90199/1‑250 (MQ=255)                       aGCTGCTTAACGCTAGTGCCGACGATACCCTTCGCGGCGGCGCAAAAACCTTTGCTGAAACGGACACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACCGGCGCGTGGTTCTCAGTTGGTGAGCGATGCTGCATCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCCTGGc                                                                                                                                                                                                                   >  2:48919/1‑251 (MQ=255)                           gctTAACGCACGTGCCGACGTTACGCTTGGCGGCGGAGAAAAAACTTTTGCTGAACCGGCAACCGCTGGTGAATGGCAGGGAAAAAAGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAATCCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACTGAAAGCAACGTACCATGGCAata                                                                                                                                                                                                               <  1:131421/251‑1 (MQ=255)                                cgtcctatctcGACGTTCGGCTTGGCGGCCGCGCAAAAACCGTGGCGGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAGCAGGCACAGGCGCGTGGTGATCAGTTGTTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATAAGCAAAAACCCCTGCTTGGCCGGGCTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGAt                                                                                                                                                                                                          <  2:26616/241‑1 (MQ=255)                                 cGCTCGTGCCGACGGTACGCTGGGCGGCGGCGCAAAAAACTTTGCTGAAACGGCAACCGCTGGTGACTGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGAGGATAGAAGTTGGTGAGCGTTGCTGCCACACTGAATTCGGTGACGGAAACGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAAAAAGCCAGTGCGCAACCTAGGACCGAAAACAACGTACCATAGAAATATCGAt                                                                                                                                                                                                          >  2:27012/1‑250 (MQ=255)                                          cGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGAt                                                                                                                                                                                                          >  1:26616/1‑241 (MQ=255)                                                 aCGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCt                                                                                                                                                                                          <  1:116938/250‑1 (MQ=255)                                                       ggcggcggcGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACTGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGc                                                                                                                                                                                    <  1:284098/250‑1 (MQ=255)                                                                            gCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAg                                                                                                                                                              <  2:52353/251‑1 (MQ=255)                                                                                            cTGGTGAATTGCAGGGAAAAACGCTGCGTGAACAGGCACAGTCGCGTGGTTATCCGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCCATATTCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAAt                                                                                                                                                                               >  2:119525/1‑218 (MQ=255)                                                                                            cTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAAt                                                                                                                                                                               <  1:119525/218‑1 (MQ=255)                                                                                             tgatGAATGGCAGGGAAAAACGCAGCGTGCACAGGCCCAGGCGCGTGGTTATCAGTTGGCGAGAGATGCTGCATCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAAACCCTGCATGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCACCGTACCATGGCAATATCGATAAGCACGCAGTCACCTGAACGCCAAATACGCAACGTAATGACAGTGTACCAACACTGgcgc                                                                                                                                             >  2:297017/4‑251 (MQ=255)                                                                                              ggTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCa                                                                                                                                            <  1:333172/251‑1 (MQ=255)                                                                                                   aTGGCAGGGAAAAACGCTGCGTGAACAGGCACCGGCGCGGGGTTATCAGTTGGGGAGCGATGCTGCCTAACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGCTAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATg                                                                                                                                        >  2:58306/1‑250 (MQ=255)                                                                                                   aTGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGa                                                                                                                                       <  1:152964/251‑1 (MQ=255)                                                                                                              aaaCGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGcc                                                                                                                            >  1:282284/1‑251 (MQ=255)                                                                                                                      cGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGCGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACCAAGCCATTGAAtt                                                                                                                    >  2:452531/1‑251 (MQ=255)                                                                                                                       gTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCgctagcta                                                                                                                                                                                                                                           >  1:132548/1‑131 (MQ=255)                                                                                                                       gTGAACAGGCACAGGCGCGTGGTGATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCgctagcta                                                                                                                                                                                                                                           <  2:132548/131‑1 (MQ=255)                                                                                                                         gAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAAttgtt                                                                                                                 <  1:477217/251‑1 (MQ=255)                                                                                                                           acaggcacaggcGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGa                                                                                                               <  2:282284/251‑1 (MQ=255)                                                                                                                            caggcacaggcGCGAGGTTATCAGTTGGTGAGAGATGCTGCCTCACTGACTTCGGTGAAGGAAGCGAATCAGCAAAAACCACTGCTTGGCCTGTTTGCTGCCGGCCATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCCAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACAGACAAAGCCATTGAAttgttgcg                                                                                                              >  2:234239/1‑249 (MQ=255)                                                                                                                             aggcacaggcGCGTGGTTATCAGTTGGTGAGCGCTGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCCGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGt                                                                                                             >  1:72182/1‑251 (MQ=255)                                                                                                                                      cgcgTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGAGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATg                                                                                                      >  2:456171/1‑249 (MQ=255)                                                                                                                                           ggTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAgg                                                                                               <  1:208303/251‑1 (MQ=255)                                                                                                                                                  aGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTc                                                                                        >  1:465611/1‑251 (MQ=255)                                                                                                                                                      ggTGAGCGATGCTGCCTCACTGAAGTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTg                                                                                     <  2:465611/250‑1 (MQ=255)                                                                                                                                                         gAGCGATGCTGCCTCACGGAATGCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAg                                                                                 <  2:90199/251‑1 (MQ=255)                                                                                                                                                                    ccTCACTTAATTCGGTGACTGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACTGCAATATGCCAGTGCGCTAGCTAGGTCCGAAAGCAACTTACCATGGCAATATCGAGAAGCCCGCAGTCCCCTGTACGCCAAATCCGCAACGTAATTACAGTGTAGCAACCCTGGCGCAGATGACCGACAATGCCATTGAATTGTTGAGTAAAAATGAGAAAGCCTTTTTCCTGCAAGTTGAAGGTGCg                                                                      <  2:541828/251‑1 (MQ=255)                                                                                                                                                                        aCTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCa                                                                   >  1:113287/1‑250 (MQ=255)                                                                                                                                                                                      cGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGACAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGAt                                                    >  2:362877/1‑251 (MQ=255)                                                                                                                                                                                       ggAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGATAGAGAGCTAGCTCGAAACGAAAGCAACGTACAATGGCACTATAGATAAACCCGCAGTCACCTGTACGCCAAATACGAAACGTAATGACAGTGTACCAACCCTGGAGCAGATGACCGACAAAGACATTGAAATGTTGAGTAAAAATGAGAAAGGATTTTTAATGCAAGTTGAAGGAACGTCAATCGATAAACAGGATc                                                   >  1:502836/1‑251 (MQ=255)                                                                                                                                                                                           gCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCAtgc                                               >  2:344863/1‑251 (MQ=255)                                                                                                                                                                                            cGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCAtgct                                              <  2:389925/251‑1 (MQ=255)                                                                                                                                                                                                        aaCCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTtgt                                  <  1:456171/251‑1 (MQ=255)                                                                                                                                                                                                           cccTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCGAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCATATGACCGACAAATCCATTTAATTGTTGAGTAAAAAGGAGAAAGGCTTTGTCCTGCAGGTTTAATGTGCGTCAATGGATAAACATGCTCATGCTGCGGATCCTTGTggg                               <  2:307949/251‑1 (MQ=255)                                                                                                                                                                                                               gCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAAGTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCaa                            <  1:303960/250‑1 (MQ=255)                                                                                                                                                                                                                cTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAAt                          >  1:403901/1‑251 (MQ=255)                                                                                                                                                                                                                   gtccTGTTTGCTGATTGCCATATGCTAGTGCGCTAGCTAGGACCTAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTGCTCCAAATCCGCACCGGTATGACAGGGTACCAACCCTTACGCAGATGACCGACAAAGCCAGTGAATTGTTGAGTATAAATGAGATAGGTGTTTTCCGGCAAGTTGAAGGTGCGGCAATCGATAAACAGGATCATGCTTCGAATCCTTTTTGGCAAATTg                        <  2:348485/248‑1 (MQ=255)                                                                                                                                                                                                                               gACGGCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTc             <  1:362877/249‑1 (MQ=255)                                                                                                                                                                                                                                   gCAATATGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGAtctc       >  1:531836/1‑251 (MQ=255)                                                                                                                                                                                                                                         tGCCAGTGCGCTAGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGATCTCGATGa  >  1:2631/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              GGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGATCTCGATGA  >  NC_000913/402371‑402852

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: ‑