Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,730,870	T→C	noncoding (288/1542 nt)	rrsG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,730,870	0	T	C	97.6%	89.2 / NA	43	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base C (17/24);  minor base A (0/1);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTT  >  NC_000913/2730708‑2731108                                                                                                                                                                    |                                                                                                                                                                                                                                               cctccCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACCAGCTAATCCCATCTGCGCACATCCGATTGCAAGAGGCCCGCAGGTccccc                                                                                                                                                        >  2:164475/1‑251 (MQ=9)   cccccGTTGAAAGTACTTTACAACCCGAAGGCCTGCTTCATACCCGCGGCAGGGCTGCATCAGGCTGGCGCCCAGTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGG‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCCGGGCTCGTCGCCTTGGTGAGCCGTTCCCCCACCAACAAGCTAATCCCATCTTGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                                                                                                                      <  2:444377/250‑1 (MQ=11)    ccccGCTGAAAGTTCTTTACACCCCGAAGGCCTTCTGCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCct                                                                                                                                                       <  1:165463/249‑1 (MQ=9)       cGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg                                                                                                                                                  <  1:551345/251‑1 (MQ=14)        gCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTgg                                                                                                                                                 <  1:182796/251‑1 (MQ=14)             aaGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTg                                                                                                                                            <  1:564253/251‑1 (MQ=14)             aaGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGATGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCCCATCCGATGGCAAGAGGCCCGCAGGTCCCCCTCTTTGGACTTg                                                                                                                                            >  2:209449/1‑251 (MQ=9)                 aCTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGAc                                                                                                                                        >  1:162897/1‑251 (MQ=14)                           ccGAAGGCCTTCTTCATACACGCGGCATGGCTTCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGATTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCGGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGt                                                                                                                              <  1:265952/251‑1 (MQ=16)                           ccGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCCCCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACCTTATGCGGt                                                                                                                              >  2:313749/1‑251 (MQ=14)                                 gCCTTTTTCATACACGCGGCATGGCTGCATCAGGCTGGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                                                                                                                        <  2:304120/251‑1 (MQ=18)                                  ccgtcttcATACACGCGGCATTGCTGCATCAGTCTGGCGCACATTGCGCAATATTCCCCACTTCTTATTCCCGGAGGATTCTGGACCG‑GTGCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGAGCAGCCAGGGATCGTCGCCGTGGTGAGCCGTTACCCCACCAACAAGCTAAGCCCATCTGGTCATATCCGAGGGCAAGAGGCCTGAAGGGCCCCCTCTTTGGTCTTGCGACATTATGGGGTATTAGCt                                                                                                                       <  2:444304/248‑1 (MQ=2)                                  ccTTCTTCATTCACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGGGCAATAGTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAGCCCATCTGGGCACATCCGATTGCAAGAGTCCCGAAGGGCCCCCTCTTTGGTCTTGCGCCATTATGCGGTATTAGCt                                                                                                                       <  2:485040/251‑1 (MQ=11)                                         cATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGttt                                                                                                                >  1:390689/1‑251 (MQ=18)                                          aTACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCGCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAACCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGCGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAGGCGGAATTAGCTACCGTTTc                                                                                                               >  1:321046/1‑251 (MQ=255)                                             cacGCGGCAGTGCTCCAGCAGGATTGCGCCCCTTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCAGGACCGTGT‑CTCAGTTTCAGGGTGGCTGGTCATCCACTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAGGCGGTATTAGCTACCGTTTCCag                                                                                                            <  2:565096/251‑1 (MQ=11)                                                gcgGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtag                                                                                                         <  2:557975/251‑1 (MQ=18)                                                gcgGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtag                                                                                                         <  1:168154/251‑1 (MQ=18)                                                 cgGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtagt                                                                                                        >  1:302150/1‑251 (MQ=18)                                                  ggCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGtt                                                                                                       >  1:454716/1‑251 (MQ=18)                                                    cATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAt                                                                                                     <  1:560694/251‑1 (MQ=18)                                                     aTGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAT‑AGCTACCGTTTCCAGTAGTTATcc                                                                                                   >  1:82407/1‑251 (MQ=18)                                                       ggCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGCAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc                                                                                                  >  1:391176/1‑251 (MQ=18)                                                             aTCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTcc                                                                                             >  1:76117/1‑250 (MQ=18)                                                               cAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATc                                                                                          <  1:367979/251‑1 (MQ=18)                                                                          ccATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTcc                                                                               >  1:266422/1‑251 (MQ=21)                                                                          ccATTGTGCAAAATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTTGCTTGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGGGCCGTTACCCCACCAACAAGCCAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTcc                                                                               <  1:430091/251‑1 (MQ=255)                                                                            aTTGTGCCATATTCCCCCCTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCa                                                                             <  2:163368/251‑1 (MQ=21)                                                                              tgtgCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAg                                                                            <  2:162897/250‑1 (MQ=25)                                                                               gtgCAATAGTCCCCCCTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCGAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGAc                                                                          <  2:10704/251‑1 (MQ=255)                                                                                  cAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCCGGCAGCTTCCCAGACAtt                                                                       >  2:564601/1‑251 (MQ=21)                                                                                     tatTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCAAATCCCATCTGGGCACATCCGATGGCTAGAGGCCCGAAGGGCCCCCTCTTTGACCTTGCGACATTATGCGGTATTTGCTAACGTTACCAGTCGTTATCNCCCTCCATCCGGCAGCTTCTCAGACAttcat                                                                    >  2:264741/1‑248 (MQ=11)                                                                                      atTCCCCGCTGCTGCCTCCCGTAGGATTCTGGACCGTGT‑CTCAGTTCCAGTTGGGCTGGTCATCCGCGCAGACCAGCCAGGGATCGTCGCCTTGTTGAGCCGTTACCCCCCCAACAAGCTACTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCGTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTc                                                                   <  2:47289/251‑1 (MQ=11)                                                                                              ctgctgCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTc                                                                                              >  1:478633/1‑216 (MQ=18)                                                                                              ctgctgCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAGCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCGTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAGCCCCCTc                                                                                              <  2:478633/216‑1 (MQ=255)                                                                                                    ccTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACt                                                     >  1:352642/1‑251 (MQ=21)                                                                                                             ggAGTCTGGCCCGTGT‑CTCAGTTCCAGTGGGGCTGGTCATCCTCCCAGACCCGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAACCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCCCCATCAGGCAGCTTCCCAGACATTACCCACCCGTCCGCCACTCGTCAGCGa                                            <  1:522452/251‑1 (MQ=11)                                                                                                                    ggctcGTGT‑TTCAGTGCCAGTGTGGGGGGTCATCCTCGCAGCCCAGCCAGGGATCTTCCCCTTGGATTGCCGTTACCCCCCCAACAAGCTAATCCCATCGGGGCATCTCCGATGGCAAGAGGCCCGAAGGGCCCCCCCTTTGGTCTGCCGACATTATTCGGTATTGGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTACCCATACTTTACTCACCCGCCCGCCACTTGTCAGCGAagcagc                                      <  1:426305/246‑1 (MQ=2)                                                                                                                     gACCGGGT‑CTCAGTACCAGTGGGGCGGGCCCTCCTCCCAGCCCGGCAAGGAATCGTCGCCTTGGTGGGCCGTTCCCCCCCCAACAACCTAATCCCAGCTGGCCACATCCGAGGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAGGCGGTCTTAGCTACCGTTTCCAGTAGTTATCCCCCTCCACCAGGCAGCTTCCCAGACATTTCTCCCCCATCCGCCACTCGTCAGCGAAGCAGCaa                                    <  2:212610/251‑1 (MQ=11)                                                                                                                                     ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGCAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCCGACATTACTCACCCGTCCGCCACTCGTCAgaga                                            >  2:356535/1‑225 (MQ=255)                                                                                                                                     ccAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGa                                            <  1:356535/228‑1 (MQ=25)                                                                                                                                      cAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGCCATTATGCGGTATTAGCTCCCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTAcc                    <  1:136955/251‑1 (MQ=255)                                                                                                                                                         ctcAGACCAGCCAGGGATCGTCGCCTTGGTGCGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGtt  >  1:196558/1‑250 (MQ=255)                                                                                                                                                                    |                                                                                                                                                                                                                                               CCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGT‑CTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTT  >  NC_000913/2730708‑2731108

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