Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,728,754	C→T	A136A (GCC→GCT)	rnt →	ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,728,754	0	C	T	97.3%	115.3 / ‑4.0	38	A136A (GCC→GCT)	rnt	ribonuclease T
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (23/13);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01

CCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGCCGCCGAAGAGGTGTAATCG  >  NC_000913/1728513‑1728997                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    ccacATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGCATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGAAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGaa                                                                                                                                                                                                                                             >  2:220611/1‑250 (MQ=255)   ccatttcatttcacagccGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCACGCCTTCAACGGCATTTACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAAGACGAGGCGCTGCACGAAATTTGTAAATTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATGATGGTGGCGCACAATGCCAAGTTTGATCACAGCTTTATGATG‑GCCGCCTCAGAACGCTCTTCCCTGAAACg                                                                                                                                                                                                                                          <  2:253582/234‑1 (MQ=255)                                    gcAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGacac                                                                                                                                                                                                         >  1:185683/1‑251 (MQ=255)                                      aaaTTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGTTAGCGGCGCGGTCAGCGAATCCGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGAAG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGATCCACCCTTTCGCCACTTTTGACActg                                                                                                                                                                                                       >  1:188915/1‑251 (MQ=255)                                        aTTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACActgct                                                                                                                                                                                                     <  1:163215/251‑1 (MQ=255)                                             cAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACt                                                                                                                                                                                                <  1:243445/251‑1 (MQ=255)                                              aaCCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTAGCCACTTTTGACACTGCTGCACTg                                                                                                                                                                                               >  1:282855/1‑251 (MQ=255)                                              aaCCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATAAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTg                                                                                                                                                                                               >  1:169135/1‑251 (MQ=255)                                                              cTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCCCTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCg                                                                                                                                                                               >  1:341451/1‑251 (MQ=255)                                                                tCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCgg                                                                                                                                                                              <  1:163375/250‑1 (MQ=255)                                                                    cGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACaaa                                                                                                                                                                         <  1:182490/251‑1 (MQ=255)                                                                          tGACCTGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGCTGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGACTGTAACCGCACCATTATGGTGGCGCACAATGCCAATTTTCATCACAGCTTTATGATG‑GCCGACGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACCCTGCTGCACTGGCCGGGCTGGCACCCGGACCAACCGTAt                                                                                                                                                                   >  2:130564/1‑251 (MQ=255)                                                                                gAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCa                                                                                                                                                              >  1:369324/1‑250 (MQ=255)                                                                                                  cgcggTCAGAGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGGGCGGTTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACCGCTTTATGATG‑GCCGCCGCAGAACGCGCATCCCTGAAACGTATCCCGTTCCACCCGTTCGCCACTTTTGACACTGCTGCACTGGCCGGCCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTg                                                                                                                                           <  1:11661/251‑1 (MQ=255)                                                                                                    cggTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGCCAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGc                                                                                                                                         >  1:177311/1‑251 (MQ=255)                                                                                                    cggTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGc                                                                                                                                         >  1:353662/1‑251 (MQ=255)                                                                                                                         tGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCAccc                                                                                                                      >  1:1090/1‑249 (MQ=255)                                                                                                                                              tACGTAAAGGAATTAAAACGAGAGGCTGAAACCGTGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAGCGCGCTTCACTGAACCGTAACCCATTCCAACCTTTCGCCACTTTTGTCACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTATACAGCACCCCGGCGCACCCCGAGCTGTACGac                                                                                               >  2:331629/1‑251 (MQ=255)                                                                                                                                                 gTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAACCGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGCCAAACAGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGAGCCCTCCGCGCTGTACGACAcc                                                                                            >  2:72615/1‑251 (MQ=255)                                                                                                                                                  tAAAGGTATTAAAGCGAGCGGCTGTAATCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCg                                                                                           >  1:178212/1‑251 (MQ=255)                                                                                                                                                               gcgagcgGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGCCAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCActgctgg                                                                              >  1:304595/1‑250 (MQ=255)                                                                                                                                                               gcgagcgGCTGTAACCGCGACATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCAGCGCTGTACGACACCGCACGCACtgctgt                                                                              >  1:304616/1‑251 (MQ=255)                                                                                                                                                                 gagcgGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACtgctgtgc                                                                            >  1:113702/1‑251 (MQ=255)                                                                                                                                                                          tAACCGCGCCATTATGATGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTAAGACACCGAACGCACTGCTGTGCTGTTTtgtg                                                                   >  1:89778/1‑251 (MQ=255)                                                                                                                                                                                  ccATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTGCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTc                                                           <  1:28070/251‑1 (MQ=255)                                                                                                                                                                                  ccATTAGGGGGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTc                                                           <  2:290556/251‑1 (MQ=255)                                                                                                                                                                                       atgtgggCGCAAAATCCCATTTTTGATCACAGCTTTATGATG‑GCCGCTGCAGAACGCGCTTCATTGAAACGGAACCCGTTCCACCCTGTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCg                                                      <  2:249383/246‑1 (MQ=255)                                                                                                                                                                                         ggtggCGCACAAGGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCt                                                    <  1:177677/251‑1 (MQ=255)                                                                                                                                                                                           gggCGCACAATGCCAATGTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTgg                                                  <  1:42054/250‑1 (MQ=255)                                                                                                                                                                                               gcACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGAGCTTCAATGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGACTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCATCAACCGCTGGAAAc                                              >  1:326331/1‑251 (MQ=255)                                                                                                                                                                                                     tGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAACGGCTTGCCAGACCGCTGGCATGGGCTTCGACAGCACCCAGGCGCACTCCGCGCCGTACGACCCCGAACGCACTGCTGGGCTGTTTTGTGAAACCGTCAACCGCTGGAAACGTCTgg                                        >  2:237638/1‑251 (MQ=255)                                                                                                                                                                                                          aTTTTGATCACA‑TGTTATGATGTGCCGCCGCAGAACTCTCTTCATTGAAACGTAACCCGTTCCACCCTTTCGCCACGTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGc                                   <  2:151908/251‑1 (MQ=255)                                                                                                                                                                                                               gATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGcc                              >  1:330855/1‑251 (MQ=255)                                                                                                                                                                                                                   acaGCTTTATGATG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCt                           >  1:151059/1‑250 (MQ=255)                                                                                                                                                                                                                      gCTTTATGAGG‑GCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCt                       <  2:185683/251‑1 (MQ=255)                                                                                                                                                                                                                                  gccgccGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGCCGCCGa               >  1:152363/1‑248 (MQ=255)                                                                                                                                                                                                                                      ccgcAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGCCGCCGAAGAGgtg        <  1:63191/251‑1 (MQ=255)                                                                                                                                                                                                                                              cgcgCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACAGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGCCGCCGAAGAGGTGAAATCg  >  1:239144/1‑249 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    CCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATG‑GCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGCCGCCGAAGAGGTGTAATCG  >  NC_000913/1728513‑1728997

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