Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,591,974	G→A	A106A (GCG→GCA)	tsr →	methyl‑accepting chemotaxis protein Tsr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,591,974	0	G	A	100.0%	113.5 / NA	41	A106A (GCG→GCA)	tsr	methyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/14);  total (27/14)

TTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGG‑GATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCC  >  NC_000913/4591744‑4592217                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     tttAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCAACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGGGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGGAAAAAACTGGGCAGCTTACGAAGCGTTGCcgcg                                                                                                                                                                                                                                  >  2:272646/1‑251 (MQ=255)     atgtcttcaaGAATGACAAAGATAATTTCACTGTTTTAACAACCATTCGCCAGCAGCAATCCACGCTTAATGGCAGCTGGGTCGCGTCGTTGCAGACGCGTAACACCCGCCACCGCGCGGTTATTCGCGACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCAAGTATTTCGCTGAAACAGGCGGAAAAAAATTGGGCAGATTACGAAGCGTTGCCGCGTGAc                                                                                                                                                                                                                              <  1:224240/243‑1 (MQ=255)        ccTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCcg                                                                                                                                                                                                                           <  1:68567/251‑1 (MQ=255)                      aaaGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCAccgc                                                                                                                                                                                                             >  1:389030/1‑251 (MQ=255)                       aaGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATAGCAGCTGGGTCACGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATAAGAATAATATTGGTAGCGGATCAA‑CCGTTGCTGA‑GCTGATGCCGAGTGCCAGTATTTCGCTGAAAAAGGCGGAACAAAACTGGGCAGATTACGAAGCGTTGCCGAGAGACCCGCGACAGTGAAccgcc                                                                                                                                                                                                            >  2:119096/1‑251 (MQ=255)                            atcctggtgctctcccgatctaTTCCCCAGCAGCAATCCCCGCTGAATGCCAGCGGGGTCGCTTTGTTGCAGACGCGTAACACCCCCACCCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGTAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGc                                                                                                                                                                                                       <  1:316536/230‑1 (MQ=255)                                    tGTTTTACAAACCATTCGCCAGCAGCAATCCAAGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGAGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATCTTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAGTCGTGGCCGCGTGACCCGCGTCAGAGCACCGCCGCGGCGGCagagg                                                                                                                                                                                               >  1:121114/1‑250 (MQ=255)                                       tttACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCTTTGTTGCAGACGCGTAACCCCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACGGGCGGAAAAAAACGGGGTAGATTCCGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATc                                                                                                                                                                                             <  1:40235/250‑1 (MQ=255)                                       tttACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATc                                                                                                                                                                                             >  1:70547/1‑250 (MQ=255)                                        ttACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCCACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCa                                                                                                                                                                                            >  1:424487/1‑250 (MQ=255)                                            aaaCCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAAAGCACCGCCGCAGCGGCAGAGATCaaccg                                                                                                                                                                                        >  1:267046/1‑247 (MQ=255)                                                cATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAAtt                                                                                                                                                                                   >  1:135233/1‑251 (MQ=255)                                                cATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAAtt                                                                                                                                                                                   >  2:224760/1‑251 (MQ=255)                                                 aTTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCCACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGAGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGc                                                                                                                                                                                                       >  2:316536/1‑230 (MQ=255)                                                      ccAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGAAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTCCATGATGGATCAGAATCATATTGGAAGCGGTTCAC‑CCGTTGCTGA‑GCTGATGCAGAGAGCCAGTATTTCGCTGAAACAGGAGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTATGAGCACCGCCGCAGCGGAAAAGATCAAACAGAATTACGat                                                                                                                                                                              >  2:271322/1‑250 (MQ=255)                                                              aaTCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATCTTTATcac                                                                                                                                                                     >  1:127279/1‑251 (MQ=255)                                                                  cACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCTAAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCATTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATg                                                                                                                                                                 <  1:205198/251‑1 (MQ=255)                                                                                     gTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATcc                                                                                                                                               >  1:87268/1‑250 (MQ=255)                                                                                             gttgCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACAGGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTa                                                                                                                                       <  1:19553/250‑1 (MQ=255)                                                                                                 cAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGc                                                                                                                                  >  1:212250/1‑251 (MQ=255)                                                                                                             acCCTCAACCGCGGGGGTATCCGCGACATGATGGATAAGAATAATATTGGTGGTGGTTCAC‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGTTTACGAAGCGTTGCCGCGTGACCCTCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCaa                                                                                                                      <  1:229606/251‑1 (MQ=255)                                                                                                                          agggaaTCCGCTACATGATGGATCAGAGTAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATAAAACGTAATGACGATATTTATCACAATGCGCTGGCGGAGATGATCCACCTATTAGGTGCAGGCAAAATCAACGAGTACTTTGAt                                                                                                         >  2:104732/6‑251 (MQ=255)                                                                                                                                  cGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGAc                                                                                                 >  1:288509/1‑251 (MQ=255)                                                                                                                                   gCTAAATGTTGGATCAGAAAAATATGGGTAGCG‑TTCAAACCGTTGCTGA‑GCGGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATGACGAAGGGTGGCCGCGTGCCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGAcc                                                                                                <  2:135233/251‑1 (MQ=255)                                                                                                                                     tACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGAccc                                                                                               >  1:69395/1‑250 (MQ=255)                                                                                                                                          gttgGATCAAAATAAATTTGGAAGCGGTTCAA‑CCGT‑GGTGAAGCGGATGCAGAGTGCCAGTATTTCGTTGAACAAGGAGAAAAAAAACTGGGAAGATTCGGAGGCGTGGCCGCGTAACCCGCGTCAAAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACGATTAGGTGCATGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGG‑Ga                                                                                         <  2:305767/248‑1 (MQ=255)                                                                                                                                           agggATCAGAATAATATTGGAAGCGGTTAAA‑CAGTTGCTGG‑GCGGAGGAAGAGTGCAAGTATTTAGCTGAAAAGGGAGGAAAAAAATGGGGAAGATTAGGAAGCGTTGCCGATTGACCCGCGTCAGAGAACCGTCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTGCTTTGATCAGCCGACCCAGG‑Gata                                                                                       <  1:336104/249‑1 (MQ=255)                                                                                                                                              gataaGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGCAACAGGCGGACAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCCACTATTAGGTGCAGGCAAAATCAACGGGTTCTTTGCTCAGCCCACCAAGG‑GATATcg                                                                                    >  2:222675/5‑250 (MQ=255)                                                                                                                                              gATCAGAATAATATTGGTAGCGTTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGG‑GATATCa                                                                                    <  1:224760/251‑1 (MQ=255)                                                                                                                                                  aGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCCGG‑GATATCAGGAc                                                                                >  1:350704/1‑251 (MQ=255)                                                                                                                                                   gAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGG‑GATATCAGGACg                                                                               >  1:58504/1‑251 (MQ=255)                                                                                                                                                     aaaataGTGGTAGCGGTCCAA‑CCGTGGCGGA‑GCTGAGGCAGGGGGCCAGTATTTCGCTGAAACAGGCGAAAAAAAACGGGGCAGATGACGAAGCGTGCCCGCGTGCCCCGCGCCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGG‑GATATCAGGACGGt                                                                             <  2:87268/249‑1 (MQ=255)                                                                                                                                                              gTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGG‑GATATCAGGACGGTTTCgaga                                                                      >  1:451408/1‑249 (MQ=255)                                                                                                                                                                            cGTTGCTGA‑GCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGCAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAACCGTAATTACGGTATTTATCACACTGCGCTGGCGGAGCTGATCCCACTATTAGGTGCAGGCACCATCCACGGGTTCTTTGATCGGCAGAACCAGGCGA‑CTCAGGTGGGTTTCGAGAAGCAGTATGTGGCtt                                                       >  2:376678/1‑251 (MQ=255)                                                                                                                                                                                gctga‑gctgaTGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGG‑GATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACAt                                                   >  1:391541/1‑251 (MQ=255)                                                                                                                                                                                                agTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGG‑GATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAGAACGAt                                      >  2:18625/1‑249 (MQ=255)                                                                                                                                                                                                  tGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGG‑GATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGc                                  >  1:391589/1‑251 (MQ=255)                                                                                                                                                                                                   gCCAGTATTTCGCTAAAACAGGCGGAAAAAAACGGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGG‑GATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCt                                 <  1:424314/251‑1 (MQ=255)                                                                                                                                                                                                       gaagttctcggAAACAGGGGGAAAAAAACTGGGCATCTGTCGAAGCGTTGCCGGTTGACCCGTGTCAGATCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAAGGCGCTGGTGGAGCTGATCCAACTATTAGGTGCAGGCAAGAGCGACGAGTTCGTTGGTCAGCCGACCCAGG‑GATATCTGTACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTctat                             <  1:26425/241‑4 (MQ=255)                                                                                                                                                                                                                                 aaCTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGTTCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGG‑GATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGc   >  1:140855/1‑251 (MQ=255)                                                                                                                                                                                                                                  aCTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGACGCAGAGGCAGAGATCAAACGTAATTACGATATTTATCACAATGAGCCGGCGGAGCTGATCCAACTATAAGGTGCAGACAAAATCAACGAGTTCTTTGATCAGCCGACCCAGG‑GATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGcc  >  1:465826/1‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     TTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAA‑CCGTTGCTGA‑GCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGG‑GATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCC  >  NC_000913/4591744‑4592217

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