Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	6,239	G→T	A74E (GCG→GAG)	yaaA ←	DNA binding and peroxide stress response protein YaaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	6,239	0	G	T	89.4%	119.0 / ‑5.3	49	A74E (GCG→GAG)	yaaA	DNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base T (29/13);  minor base C (1/4);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 5.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTT  >  NC_000913/5992‑6481                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   gTTCAGCTTGTTGGTGAGGATCTCTCCCCAGAATTGAAACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCGcc                                                                                                                                                                                                                                                  <  2:18648/250‑1 (MQ=255)    cAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGc                                                                                                                                                                                                                                               <  1:24666/250‑1 (MQ=255)      gCTTGTTGGTGATGATATCCCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGCCGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTg                                                                                                                                                                                                                                            <  2:243869/251‑1 (MQ=255)                  tgatATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCCCCGGGGTGAAATCTGGCTGc                                                                                                                                                                                                                                >  2:409484/1‑251 (MQ=255)                                 aTTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAg                                                                                                                                                                                                                             >  1:323797/1‑239 (MQ=255)                                 aTTGATACAGATCTTCCCCTCGGGCATTCCCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAg                                                                                                                                                                                                                             <  2:323797/239‑1 (MQ=255)                                   tGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAagcgagc                                                                                                                                                                                                               >  1:58918/1‑251 (MQ=255)                                     aTACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAACGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAgcgg                                                                                                                                                                                                             >  1:374535/1‑251 (MQ=255)                                       aCAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACTCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGATTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAAAGCCAGAATCGCCTGGAGGGCAT‑TCTCCGACATGAAATCTGGGTGACAGTCATGaaggagacctgcg                                                                                                                                                                                                           >  2:373642/1‑240 (MQ=255)                                        cAGATCTTTCCATCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGGGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGt                                                                                                                                                                                                          <  1:373642/251‑1 (MQ=255)                                            tCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAACCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGa                                                                                                                                                                                                       <  2:8798/250‑1 (MQ=255)                                               ttCCCTCGGGCATTCTCAAGACGGATCCCAATATCCAGACGAATAGGCTGCATTAAAGCGAGGGGTCGCAGTACGCCAGACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTGTCGGCCTGCAAGCTGGTGTAGACATCACCTTTAAACTCCAGAATCGCCTGGCGGGCAG‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATAc                                                                                                                                                                                                    <  2:482593/250‑1 (MQ=255)                                               ttCCCCCGGGCTTTCTCAAGACGGATCCCCTTTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATa                                                                                                                                                                                                     <  2:140159/249‑1 (MQ=255)                                               ttCCCCCGGGCATTCCCAAGGCGGACCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAACGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATa                                                                                                                                                                                                     <  2:374535/249‑1 (MQ=255)                                                tCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATAccc                                                                                                                                                                                                  >  1:170476/1‑251 (MQ=255)                                                  ccTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGACGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTTGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGc                                                                                                                                                                                                >  1:241193/1‑251 (MQ=255)                                                   cTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGcc                                                                                                                                                                                               >  1:457179/1‑251 (MQ=255)                                                    ccGGGCATCCTCAAGCCGGTTCCCAATTTCCAGGCGAAAAGGCGGCATTAAATCGGGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTGGGGCAAAATCGAAACCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAG‑TCTCCGGCGTGAAATCTGGCGGCCAGTCAGGAAAGCGAGCGGCGTTGATACCCGCCa                                                                                                                                                                                              <  2:479206/250‑1 (MQ=255)                                                         cATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTg                                                                                                                                                                                         >  1:376440/1‑251 (MQ=255)                                                          aTTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTATATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGAAACCCGCCCGTTTGt                                                                                                                                                                                        >  2:326153/1‑251 (MQ=255)                                                            tctcAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAACTCGAAATCGTCTTCGCTGAAGGTTTCGGCATGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCACCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGACCGGCGTTGATACCCGCCAGTTTGTCg                                                                                                                                                                                      >  1:227567/1‑251 (MQ=255)                                                                  gACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCCGTTTGTCGCTGATg                                                                                                                                                                                >  1:167076/1‑251 (MQ=255)                                                                            cATTCCCAGACGATAAGTCTGCATTAAACCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGTTGTTGGGCAAAATCGAAATCTTCTCCGCTGAAGTTTTCGGCCTGCAAGCCGGTGTAGACATCCCCTTTAAACGCCAGAAGCGCCTGGCGGGCAT‑TCTCCGGCGTAAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCg                                                                                                                                                                      <  1:353416/251‑1 (MQ=255)                                                                                   ataCGATAAGGCTCAATTAAAGCGATCGGGCGCATTCCGCCTTACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCAAAATCGTCTTCGCTGAAGTTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAACCTGGCTGCCAGTCATGAAAGCGAGCGGGGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAAt                                                                                                                                                               <  1:119191/249‑1 (MQ=255)                                                                                       gATAAGGCTGAATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTCAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCACTGGTGCGCATCAGAGTGCTAATCTGc                                                                                                                                                           >  2:308388/1‑251 (MQ=255)                                                                                              cTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCg                                                                                                                                                    >  1:383187/1‑251 (MQ=255)                                                                                                 cATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCCTGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAggagtac                                                                                                                                                 >  2:462138/1‑246 (MQ=255)                                                                                                      aaTCGAGCGGGCGCGGTACGCCATACAAGCCGGAAAGCATTAGCAAATGCTGTTGGGCAAAATCGACATAGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGGGTAGACCTCACATTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCAGCGTTGATACCAGCCAGTTTGTCGCTGGTGCGCATCAGCGTGCTAATcggcgggggagtaaggttc                                                                                                                                            >  2:218812/1‑238 (MQ=255)                                                                                                      aaTCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTc                                                                                                                                            >  1:28122/1‑251 (MQ=255)                                                                                                           aGCGGGCGCAGTACGCCATACAAGCCGGCAAGCATTCGCAACTGCTGTTGGGCAAAATCGCAATCGTCTTCGCTGCCGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGACAGCGACCGGCGTTGATACCCGCCAGTTTGTCGCGGATGCGCATCAGCGTGCTCATCTGCGGAGGCGTCAGATTgcgcgc                                                                                                                                       >  2:21704/1‑251 (MQ=255)                                                                                                                 cgcaGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATg                                                                                                                                 >  1:216321/4‑251 (MQ=255)                                                                                                                     gTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGACTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGGTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCCGTTTCCGCGCTTCATGGATc                                                                                                                             >  1:175522/1‑251 (MQ=255)                                                                                                                              aCAACCCGTAAAGCATCCGCAATTGCTGTGTGGCACAATCGAATTCTTCTTCGCTTAAGTTTTCGCCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGACTCGCCGGGCGGGCAT‑TCTCCGGCGTAAAATCGGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgctg                                                                                                                     <  1:321392/250‑1 (MQ=255)                                                                                                                                   ccGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAAtt                                                                                                               >  1:446879/1‑251 (MQ=255)                                                                                                                                     ggAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGt                                                                                                             >  1:323322/1‑251 (MQ=255)                                                                                                                                       aaaGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCt                                                                                                           >  1:119746/1‑251 (MQ=255)                                                                                                                                               cGCAAATGCTGTGGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCAAGAATCGCCGGGGGGGTAT‑TCTCCGTCGTGAAATCTGTCTCCCAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCGACTGCTGGGAATTGTCTAACAGCTc                                                                                                   <  1:462138/251‑1 (MQ=255)                                                                                                                                                 cAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGCTGCGCATCCGCGTGCTAATCTGCGGGGGCGTCGGTTTCCGCGCTTCATGGGTTAACTGATGGGCATTGTCTGACAGCTCCg                                                                                                 >  2:84389/1‑251 (MQ=255)                                                                                                                                                    aTGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTcc                                                                                                                                           >  1:436995/1‑206 (MQ=255)                                                                                                                                                                                  gCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCACGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGGTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCttt                                                                >  2:103536/1‑251 (MQ=255)                                                                                                                                                                                   cTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTACCAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTg                                                               >  2:17798/1‑251 (MQ=255)                                                                                                                                                                                      agggttccgcccTGAAAGCGGGTGTAGACATCACCTTAAACCGCCAGAACCGCCGGGCGGGCAT‑TCCCCGGCGGGAAATCGGGCTGCCAGTCATGAAAGCGAGGGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTa                                                            <  2:93578/241‑1 (MQ=255)                                                                                                                                                                                              ggccggaaagcctgtggaaaaaTCACCTTTAAACGCCAGAATCGCCTGGGGGGCATCTC‑CCGGGGTGAAAGGGGGCTGCCATGAATGAAAGCGAGCGGCGGTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAAGTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCg                                                    <  2:175522/230‑1 (MQ=255)                                                                                                                                                                                                           gtgtaacaatccccTTTAAACGCCAGACCCGCCGGGGGGGTAT‑TTCCCGGGGGGAAAGCTGGCGGCCAGTAATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGa                                       <  2:58918/239‑1 (MQ=255)                                                                                                                                                                                                                     cactgtataaccctaGAATCGCTCGGCGGGCAT‑TCCCCGGCGTGAAATCTGCCTCCAAGTCTTGAAAGCGCGGGGCGTTGACCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATc                             <  2:387777/237‑1 (MQ=255)                                                                                                                                                                                                                        cTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGc                          >  1:325209/1‑251 (MQ=255)                                                                                                                                                                                                                           tAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgctg                                                                                                                     >  2:31334/1‑157 (MQ=255)                                                                                                                                                                                                                           tAAACGCCAGAACCGCCGGGCGGGCAT‑TCCCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgctg                                                                                                                     <  1:31334/157‑1 (MQ=255)                                                                                                                                                                                                                                                gggCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAAttt  >  1:48326/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   GTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTT  >  NC_000913/5992‑6481

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