Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,587,531	G→A	L195F (CTC→TTC)	ugpQ ←	glycerophosphodiester phosphodiesterase UgpQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,587,531	0	G	A	100.0%	134.4 / NA	43	L195F (CTC→TTC)	ugpQ	glycerophosphodiester phosphodiesterase UgpQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/20);  total (23/20)

ATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCG  >  NC_000913/3587286‑3587775                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     aTCAGCAGCAGATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCAGCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACt                                                                                                                                                                                                                                                 <  1:42862/251‑1 (MQ=255)         caTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTg                                                                                                                                                                                                                                         >  2:399411/1‑251 (MQ=255)                             aaaCGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGcagcc                                                                                                                                                                                                                     >  2:26939/1‑251 (MQ=255)                             aaaCGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGcagcc                                                                                                                                                                                                                     >  1:136594/1‑251 (MQ=255)                                                     cAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGCGACGCAGCCCAGCCGCGCGGTCAGTTCGCGcc                                                                                                                                                                                              >  2:126116/1‑250 (MQ=255)                                                       gCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAgt                                                                                                                                                                                           <  1:110814/251‑1 (MQ=255)                                                                           ccGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGc                                                                                                                                                                       >  1:380353/1‑251 (MQ=255)                                                                             gTTGAAACCTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCACCTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCACCCCAGCCGCGCGGTCAGTTCGCGCCAGACGTCGCGCCACTCATCCAGca                                                                                                                                                                      >  2:134499/1‑250 (MQ=255)                                                                                   aaCTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACCCCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAcc                                                                                                                                                               >  2:378752/1‑251 (MQ=255)                                                                                            gtcGTATAGTTCGGACCAATCACGCCAATCGCATCGGTACAGAGGCAACCCACCCCCCAGCGCAGCAACTCTGCTGCGCTCTGGGGTTTATTGACGGTATAATCCAGAATCCGCAGTCCGGCGTCGTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGAGGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCTAGTCGTCGCGCCTTTCATCCAGCAACAAACCGCGCGGCat                                                                                                                                                      <  2:217224/249‑2 (MQ=255)                                                                                                  aaGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTAGGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCgg                                                                                                                                                <  1:355729/251‑1 (MQ=255)                                                                                                        ggACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATACACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCCTCCAGCAACAAACCGCGCGGCCGTTCCGGTgccgc                                                                                                                                          >  1:241879/1‑251 (MQ=255)                                                                                                            cAATCACGTCAATCGCATCGGTACAGATGCAATCCACAACCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCtgt                                                                                                                                      >  2:177570/1‑251 (MQ=255)                                                                                                                     cAATCGCCTCGGTACAGATGCCATCAACACCCCAGCGCAACAACTCTGCTGCGCGCTGGGTTTTATTGACGGTATAAACCAGACTCCGCCGTCCGGAGTCTTTCACCTGCATCCCTCGCGCTTTATTGAATAACTTATGATTGAGCTGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCCGTTCCGGTGCCGCCtgttagatagg                                                                                                                              >  2:396028/1‑243 (MQ=255)                                                                                                                     cAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACGCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCtgtt                                                                                                                                     <  1:396028/243‑1 (MQ=255)                                                                                                                         agCATCGGTACAGATGCAATCCACAACCCAGCGCAGGAACTCAGCGGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTACGGCGTATTTCAACTGCATCACTCGCGCTTTATTGAATAACTTAGGATTGAGATGAATAGAGCCGCAGCCCAGCCGCGCGGTCAGTTCGCGCAAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTGGTGCAGCTTCTa                                                                                                                         <  1:144832/250‑1 (MQ=255)                                                                                                                                    aGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAAtc                                                                                                              <  2:136594/251‑1 (MQ=255)                                                                                                                                         cAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCaaa                                                                                                         <  1:265350/251‑1 (MQ=255)                                                                                                                                                  ccccAGCGCAGCAACTCTTCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCTGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTAGTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATTTCAAACGATGAcag                                                                                                <  1:115206/251‑1 (MQ=255)                                                                                                                                                                     ctgcGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCg                                                                             >  1:158032/1‑251 (MQ=255)                                                                                                                                                                            cTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCcac                                                                       >  1:47030/1‑250 (MQ=255)                                                                                                                                                                               gggTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACGCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCt                                                                   <  2:77152/251‑1 (MQ=255)                                                                                                                                                                                  tttATTGACGGTATAAACCAAAATCCGCAGTACGGCGTCTTTCCACTGCATCACTAGCGATTTATTGAATACCTTATGATTGAGATGAATCGAGACGCAGCACAGCCGCGCGGTCAATTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCACTCTCAAACGATGACAGCAGCGGAGGCGTCATACCGGCCCACAGCTcgc                                                                >  2:28811/1‑251 (MQ=255)                                                                                                                                                                                        gACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTgcc                                                           <  1:173669/250‑1 (MQ=255)                                                                                                                                                                                         actgTATAAACAAGAATCCGCAGTCAGGAGCCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGTGATGAAAAGAGACGCAGCCAAGCCGAGAGGTCAGTTCGCGCCAGTCGTCGCGCCACTCCTCCAGCAACGACACGCGCGGCAGTTACGGTGCCGCCTGTTGGGCAGCTTCTAAAGAATCCATATCAAACGATGACAGCAGCGGCGGAGTCGTAACGGCACACAGcacgagtgacg                                                          >  2:354888/4‑240 (MQ=255)                                                                                                                                                                                               taAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCCTCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCCGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGCCAGCAGCGGCGGCGTCATACCGGCCCCAAGCTCGCGTGCCGCCAgcgc                                                   >  2:374972/1‑251 (MQ=255)                                                                                                                                                                                                  aCCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGccac                                                <  1:173400/251‑1 (MQ=255)                                                                                                                                                                                                     aGAATCCGCAGTCCGGCGTCTTTCATCTGCATCACTCGCGCTTTATTGAATAACTTATGCTTGAGCTAAATAACGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCTACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGACTGTTGTGCATCTTCTAAAGCATCAATCTCACACGATGACAGCAGCGGCGGCGTCATACCGGCACACAGCACGCGTGCCGCCAGCGccaccct                                             >  2:401350/1‑249 (MQ=255)                                                                                                                                                                                                          ccGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCCTTTTGc                                        >  2:242950/1‑251 (MQ=255)                                                                                                                                                                                                               gTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCCACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGtt                                   <  2:172921/251‑1 (MQ=255)                                                                                                                                                                                                                  cGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAAt                                >  1:4508/1‑251 (MQ=255)                                                                                                                                                                                                                      gTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGcc                            >  2:45027/1‑251 (MQ=255)                                                                                                                                                                                                                         tttCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGt                                                              >  1:92880/1‑214 (MQ=255)                                                                                                                                                                                                                         tttCAACTGCATCACTCGCGCTTTAGTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGt                                                              <  2:92880/214‑1 (MQ=255)                                                                                                                                                                                                                          ttCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGccc                           <  1:223325/248‑1 (MQ=255)                                                                                                                                                                                                                             aaCTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccggtgcc                     <  1:326226/251‑1 (MQ=255)                                                                                                                                                                                                                                tGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCAtttt                                          <  2:309621/227‑1 (MQ=255)                                                                                                                                                                                                                                tGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCAtttt                                          >  1:309621/1‑227 (MQ=255)                                                                                                                                                                                                                                     cACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGtggtgg             <  1:242950/251‑1 (MQ=255)                                                                                                                                                                                                                                      aCTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCCTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGCGCCGGTGGTggg            >  2:316154/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aTTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCg  >  1:89278/1‑250 (MQ=255)                                                                                                                                                                                                                                                 aTTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCg  <  1:364451/250‑1 (MQ=255)                                                                                                                                                                                                                                                 aTTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCAGCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTgcc                                                           <  1:329644/193‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     ATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCG  >  NC_000913/3587286‑3587775

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