Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,040,433	C→A	A319D (GCC→GAC)	msrP →	protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,040,433	0	C	A	97.5%	122.1 / NA	40	A319D (GCC→GAC)	msrP	protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (21/18);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACG  >  NC_000913/2040194‑2040679                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       gATTGTGCAGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAACCTGCCCCGCGAGCGTCCGCCAACCACCTGGCATCTGGCAGCGCCTGCCGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACCGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGGTGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGAGGCAt                                                                                                                                                                                                                                             >  2:122224/1‑251 (MQ=255)  aTTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATc                                                                                                                                                                                                                                            <  1:259864/251‑1 (MQ=255)      tgCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTg                                                                                                                                                                                                                                        >  2:192872/1‑251 (MQ=255)                          aaaGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCg                                                                                                                                                                                                                    <  1:4751/251‑1 (MQ=255)                                aTTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGgaga                                                                                                                                                                                                               >  1:57951/1‑250 (MQ=255)                                aTTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGaa                                                                                                                                                                                                              <  1:147246/251‑1 (MQ=255)                                aTTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGaa                                                                                                                                                                                                              <  2:49797/251‑1 (MQ=255)                                  tAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAAt                                                                                                                                                                                                             <  1:248908/250‑1 (MQ=255)                                     aTCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCt                                                                                                                                                                                                         >  2:266348/1‑251 (MQ=255)                                           aGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATg                                                                                                                                                                                                    <  1:76796/250‑1 (MQ=255)                                           aGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCAACTGGAATCTGGCAGCGCCTGACGACTACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGc                                                                                                                                                                                                   >  1:161982/1‑251 (MQ=255)                                                gTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGATGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTg                                                                                                                                                                                              >  2:2225/1‑251 (MQ=255)                                                               gcgAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTg                                                                                                                                                                               >  1:224406/1‑251 (MQ=255)                                                                              ccaccTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGtt                                                                                                                                                                >  1:72088/1‑251 (MQ=255)                                                                                                    tGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGACCGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGCGGCATCGCTGTATCGTGGCCTGGCTTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGAttgt                                                                                                                                          >  2:41244/1‑251 (MQ=255)                                                                                                            aCGGTTGTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGttt                                                                                                                                  <  2:243466/251‑1 (MQ=255)                                                                                                                tttttACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTctct                                                                                                                              >  1:112235/1‑251 (MQ=255)                                                                                                                      cGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTTCCGTTTCTctggcg                                                                                                                          <  2:249041/249‑2 (MQ=255)                                                                                                                      cGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTctggct                                                                                                                          >  1:249041/1‑249 (MQ=255)                                                                                                                        ccAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCCACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTc                                                                                                                      >  1:131334/1‑251 (MQ=255)                                                                                                                         cAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACATCAAAACAGGTGACATGGCGGAAAGTTTGCCTGCATCTTGCCGGGTTGGTGCCGTTTCTCTGGCTGGTCt                                                                                                                     <  2:118183/251‑1 (MQ=255)                                                                                                                              ttAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCg                                                                                                                >  2:193082/1‑251 (MQ=255)                                                                                                                                     gTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGTTTACGCCGACCAGGTGGCATCGCTGGATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGGCTGGGCTatcaatc                                                                                                         <  2:161982/251‑1 (MQ=255)                                                                                                                                                gcccGCGCTGGGCGCAGGCTAGGGAAAGATTTATTGGTTCAGGCGGCACCCTCGATTTACAAAGCAAACCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTATGGGCGATCAATAACGGAGGACTg                                                                                              >  2:143206/2‑251 (MQ=255)                                                                                                                                                         ggTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGa                                                                                      <  1:2225/250‑1 (MQ=255)                                                                                                                                                            cacaGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCg                                                                                  <  2:207474/251‑1 (MQ=255)                                                                                                                                                            cacaGGCTACCGAACGATTTATTGGTGCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGTGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGGTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCGGGTCTGGGCGATCAATCACGGTGGCCTGGGTGCCGATCCg                                                                                  <  2:83329/251‑1 (MQ=255)                                                                                                                                                                                  tGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCAtttt                                                             <  1:7731/250‑1 (MQ=255)                                                                                                                                                                                        aGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTc                                                      >  1:244333/1‑251 (MQ=255)                                                                                                                                                                                               atcatCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCCGGTGGAATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTACATGCGTCTGACAGCAAACCAGGTGACATGGCTGAACGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGCCGATCAATCACGTTGGACTGGGTTCCGATCCGGTGAAGGATATCCAGCATTTTACTGGTCgcactgc                                               >  2:126443/5‑251 (MQ=255)                                                                                                                                                                                                     gATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGa                                          >  1:155726/1‑250 (MQ=255)                                                                                                                                                                                                              cgcCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCt                                <  1:266348/251‑1 (MQ=255)                                                                                                                                                                                                                   gccaACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCg                            >  1:195194/1‑250 (MQ=255)                                                                                                                                                                                                                      aaCGATACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGGCTGGGCGCTCAATCACGGTGGACTGGGGTCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCATTGAAATTTTTGCTGGCTACCt                        >  2:80272/1‑251 (MQ=255)                                                                                                                                                                                                                            aCTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGtt                    <  2:57951/249‑1 (MQ=255)                                                                                                                                                                                                                               gtatAATGGTTACGCCGACCAGGTTGCATCGCTGTATCGGGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTGGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCTATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATc                <  2:141622/247‑1 (MQ=255)                                                                                                                                                                                                                                 ttAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCAc              <  1:147257/250‑1 (MQ=255)                                                                                                                                                                                                                                      ggTTACGCAGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCctct        >  1:192517/1‑251 (MQ=255)                                                                                                                                                                                                                                          aCGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTTAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCGCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGTTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCa    <  1:166884/251‑1 (MQ=255)                                                                                                                                                                                                                                          aCGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGc     >  1:218900/1‑250 (MQ=255)                                                                                                                                                                                                                                            gccgacCAGGAGGCATCGCTGTATCGAGGCATGGTTTTGCGGGAGAATCTCTAAATGCGCCTGACAGCAAAACAGGTGACCTGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGCGGACTGGGTGCAGATCCGGTGAAAGATATTAAGCATTTTACCGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACg  >  2:99042/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       GATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACG  >  NC_000913/2040194‑2040679

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