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.5%	138.6 / NA	40	A136A (GCC→GCT)	rnt	ribonuclease T
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/15);  minor base G (0/1);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGCCGCCGAAGAGG  >  NC_000913/1728508‑1728989                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            ggACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCAc                                                                                                                                                                                                                                         <  1:464275/251‑1 (MQ=255)  gACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACt                                                                                                                                                                                                                                        >  1:304347/1‑251 (MQ=255)    caccacATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTg                                                                                                                                                                                                                                       >  1:351089/1‑250 (MQ=255)                  tCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGCGGCGCTGCACGAAATTTTTAAAGTTGTACGTAACGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCTCTGAAACGTAACCCGTTc                                                                                                                                                                                                                        >  2:64407/1‑251 (MQ=255)                        tCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCt                                                                                                                                                                                                                  >  1:453684/1‑251 (MQ=255)                           aaCCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTc                                                                                                                                                                                                               >  2:168855/1‑251 (MQ=255)                                  tGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACtt                                                                                                                                                                                                        >  1:465189/1‑251 (MQ=255)                                   gTCGGCACAAATTTGCACCAAGAAGCCCTGGCCTTCAACGGCATTGTCCCGAACGATCCCGATCGCGGCGCGGTCAGAGAATACGAGGCGCCGCAGGAAATTTTTAAAGTTGTACGTAAAG‑TATTTAAACCGAGCGGCTGTAACCGCGCCATTATGGGGGTGCACAAGGCCAATTTTGATCACAGCTTTATGAGGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACttt                                                                                                                                                                                                       <  2:501905/251‑1 (MQ=255)                                          cAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTc                                                                                                                                                                                                               <  2:177199/236‑1 (MQ=255)                                          cAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTc                                                                                                                                                                                                               >  1:177199/1‑236 (MQ=255)                                               ttGCAACCAGAAGCCCCCGCATTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTGGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAAGGCCATTTTTGATCACAGCTTTATGAGGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACActgctg                                                                                                                                                                                            <  2:369782/250‑1 (MQ=255)                                                       aGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCg                                                                                                                                                                                   >  2:184449/1‑251 (MQ=255)                                                           gCCCTCGCCTTCAACGGCATTGCCCCTAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGTCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCt                                                                                                                                                                               <  1:2079/251‑1 (MQ=255)                                                                      cAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGACGCTGCACGAAATTTTTAAAGTTGTACGTACAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAGACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCAATGGCCCGGCGGGCACTCGGAc                                                                                                                                                                    >  1:260937/1‑251 (MQ=255)                                                                           gCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACTAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCTCTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAAcc                                                                                                                                                               <  1:64407/251‑1 (MQ=255)                                                                                     gAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGAATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCa                                                                                                                                                      >  2:297632/1‑250 (MQ=255)                                                                                                      gcgcggTCAGAGAATCCGAGGCGTGGCACGAAATTTTTAAAGTTGTACGTAAAGGTTTT‑AAAGCGTGCGGCGGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCt                                                                                                                                    <  2:281717/251‑1 (MQ=255)                                                                                                          ggTCAGCGAATACGAGGCTCTGCTCGAAATTTTTAAAGTTTTACGTAAAGGTATT‑AAAGCGAGCGTCTGTAACCGCTCCATTATGGTGGCGCCCAATTCCATTTTTGATCACAGCTTTATGATGTCCGCCGCAGAACGCTCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCa                                                                                                                                <  1:329223/251‑1 (MQ=255)                                                                                                                   aTACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGCTGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCg                                                                                                                       >  1:494144/1‑251 (MQ=255)                                                                                                                        gggCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGc                                                                                                                  <  1:406740/250‑1 (MQ=255)                                                                                                                             cTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCCAAGGCTTGCCAGACCGCTGGCTTGGACTTCGACAGCAccc                                                                                                              >  2:183857/1‑250 (MQ=255)                                                                                                                                cACGAAATTTTTAAAGTGGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTGACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGgc                                                                                                          <  1:84620/251‑1 (MQ=255)                                                                                                                                  cGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTACACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGgcgc                                                                                                        >  1:201025/1‑251 (MQ=255)                                                                                                                                                   tACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGac                                                                                       >  2:154631/1‑251 (MQ=255)                                                                                                                                                    aCGTAAAGTTATT‑AAAGCGAGCGGCTGTACCCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTGTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGac                                                                                       <  1:205185/250‑1 (MQ=255)                                                                                                                                                    aCGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGaa                                                                                                                                                                                                                                     >  1:185820/1‑108 (MQ=255)                                                                                                                                                    aCGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGac                                                                                       >  2:113009/1‑250 (MQ=255)                                                                                                                                                    aCGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTGGATCACAGCTTTATGATGGCCGCCGCAGAACTCGCTTCACTGaa                                                                                                                                                                                                                                     <  2:185820/108‑1 (MQ=255)                                                                                                                                                        aaaGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCg                                                                                   <  1:417449/250‑1 (MQ=255)                                                                                                                                                          aGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAAc                                                                                >  1:375750/1‑251 (MQ=255)                                                                                                                                                               tt‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCAct                                                                           >  1:367204/1‑251 (MQ=255)                                                                                                                                                                t‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCActg                                                                          >  2:215884/1‑251 (MQ=255)                                                                                                                                                                t‑AAAGCGAGCGGATGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCGTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGTGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCTCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCAct                                                                           <  1:65678/250‑1 (MQ=255)                                                                                                                                                                           ggCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGtt                                                                >  2:397741/1‑251 (MQ=255)                                                                                                                                                                            gCTGTACCCGCGCCATTATGGTGGCGCACAATGCCAATTTTCATCTCAGCTTTATGAAGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCGTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGttt                                                               <  2:294672/251‑1 (MQ=255)                                                                                                                                                                                                               aaTTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAgg                            >  1:319781/1‑251 (MQ=255)                                                                                                                                                                                                                aTTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGc                           >  1:181343/1‑251 (MQ=255)                                                                                                                                                                                                                aTTTTGATCACAGATTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAAACCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCCAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGc                           >  2:257410/1‑251 (MQ=255)                                                                                                                                                                                                                          caggtttATGATGGCCGCCGCAGAACGCGCTTCACAGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACGCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAAGCGTCAACCGCTGGAAACGTCTGGGAGGCTTGCCGCTAt                 <  2:435311/247‑1 (MQ=255)                                                                                                                                                                                                                                         cgccgcAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGCCGCCGAagggg  >  2:11420/1‑248 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            GGACACCACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATT‑AAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGCCGCCGAAGAGG  >  NC_000913/1728508‑1728989

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