Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,559,621	C→T	P452S (CCC→TCC)	rtcR →	DNA‑binding transcriptional activator RtcR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,559,621	0	C	T	97.8%	144.2 / NA	46	P452S (CCC→TCC)	rtcR	DNA‑binding transcriptional activator RtcR
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (30/14);  total (32/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 = 1.00e+00

CGAACCTGGATTATGAAGTGGAGCGCCACGCCTCACTCACTGGCGACAGCGTGCGTTTTAACACCGAAGCGCGGCGCGCCTGGTTGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTG  >  NC_000913/3559376‑3559865                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      cgaaccTGGATTATGAAGTGGGGCGCCTGGCCTCACTCACTGGCGACAGCGGGCGTTTTAACACCGAAGCGCGGCGCGCCTGGTTGGCGTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAATTTTCGCGAACTTGCTGCCAGCGTCACGCGGAGGGCCCCCTTTGCCACGAGCGGAGGCATCACTCTGGACGTTGTTGAAGATGAGATAAACCGTCTGCTCGATAACTGTCAGGTGAGCCGCtcctcc                                                                                                                                                                                                                                                  <  2:274468/251‑1 (MQ=255)               gAAGTGGAGCGCCACGCCTCACTCACTGGCGACAGCGTGCGTTTTAACACCGAAGCGCGGCGCGCCTGGTTGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTTTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGtt                                                                                                                                                                                                                                    >  2:214995/1‑251 (MQ=255)               gAAGTGGAGCGCCACGCCTCACTCACTGGCGACAGCGTGCGTTTTAACACCGAAGCGCGGCGCGCCTGGTTGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTTGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGtt                                                                                                                                                                                                                                    <  2:301875/251‑1 (MQ=255)                   tGGAGCGCCACGCCTCACTCACTGGCGACAGCGTGCGTTTTAACACCGAAGCGCGGCGCGCCTGGTTGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTg                                                                                                                                                                                                                                >  1:247352/1‑251 (MQ=255)                    ggAGCGCCACGCCTCACTCACTGGCGACAGCGTGCGTTTTAACACCGAAGCGCGGCGCGCCTGGTTGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTgg                                                                                                                                                                                                                               <  2:57667/251‑1 (MQ=255)                         gcCACGCCTCACAGACGGGCCACGGCGTGCGTTTTACCACAGAAGCGAGGCGCGCCTGGTTGGCTTTTGCGACCTCACCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCCGGACGTGGTTGCAGCTGAGATCAACCGTCTGCGCTATAACTGGCTGGAGAGCCGCCCCTCCGCGCTTCCGGCGTTGCTGGgcgcc                                                                                                                                                                                                                          >  2:306918/1‑250 (MQ=255)                                          ggCGACAGCGTGCGTTTTAACACCGAAGCGCGGCGCGCCTGGTTGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATc                                                                                                                                                                                                           <  2:1109/249‑1 (MQ=255)                                            cGACAGCGTGCGTTTTAACCCCGAAGCGCGGCGCGCCTGGTTGACTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTACCTTTCGCGAACTTTCTGCCAGCGCCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGCGGCAGAAAACATCGAtc                                                                                                                                                                                                       >  1:177700/1‑251 (MQ=255)                                                      cgtTTTAACACCGAAGCGCGGCGCGCCTGGTTGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCg                                                                                                                                                                                             >  1:277739/1‑251 (MQ=255)                                                         tttAACACCGAAGCGCGGCGCGCCTGGTTGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTAt                                                                                                                                                                                          >  1:219887/1‑251 (MQ=255)                                                               acCGAAGCGCGGCGCGCCTGGTTGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACt                                                                                                                                                                                    >  2:166940/1‑251 (MQ=255)                                                                ccGAAGCGCGGCGCGCCTGGTTGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTg                                                                                                                                                                                   >  2:13202/1‑251 (MQ=255)                                                                 cGAAGCGCAGCGCGCCTGGTTGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACGTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACGCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTgg                                                                                                                                                                                  <  2:276402/251‑1 (MQ=255)                                                                                cTGGTTGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCt                                                                                                                                                                   >  1:217073/1‑250 (MQ=255)                                                                                    ttGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTg                                                                                                                                                              >  1:123940/1‑251 (MQ=255)                                                                                              cGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCa                                                                                                                                                    <  1:271349/251‑1 (MQ=255)                                                                                                  ctctCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGt                                                                                                                                                >  1:232227/1‑251 (MQ=255)                                                                                                                 aTGGCGGGGTAAGTGTTGCGAACTTTCTGCCATCCTCACGCAGATAGCCACCTTTGACACTAGCGGAGGCATCATTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGAGCTATAACTGGCAGGAGATCCGATCCTCCGCGCTTACGGCGTTGCTGCGCGCACAGGCAGAACACCTCGCGCTCTTCGACCGTATGCAACTGCAAATAGTTCTGCGC‑ATCTGCCGCACGGCAACGACGCTTTCCGCAGccc                                                                                                                                 >  2:28884/1‑250 (MQ=255)                                                                                                                   ggCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCgg                                                                                                                                >  1:186853/1‑250 (MQ=255)                                                                                                                     cgcgGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTAt                                                                                                                                                                                          >  2:208865/1‑191 (MQ=255)                                                                                                                     cgcgGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTAt                                                                                                                                                                                          <  1:208865/191‑1 (MQ=255)                                                                                                                               ttCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCtttt                                                                                                                    <  1:60659/250‑1 (MQ=255)                                                                                                                                   cgAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACg                                                                                                               >  1:226451/1‑251 (MQ=255)                                                                                                                                         ttCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTcgc                                                                                                         >  2:124356/1‑251 (MQ=255)                                                                                                                                           cTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCgcgc                                                                                                       >  1:162137/1‑251 (MQ=255)                                                                                                                                           cTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCgcgc                                                                                                       >  1:272472/1‑251 (MQ=255)                                                                                                                                             gCCAGCGTAACGCGGAAGGCCACCTTTGCAACTAGAGGACGCATCACTCTGCACGTGGTTGAAGCTGAGATAAACCATCTGCGCTCTCACTGGCAGCCGCGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAATCATCGATCTCTTCGACCGCATGCAACTGGAACACGTTAT‑CGCTATCTGCCGTCAGGCAAAGTCGCCTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGcc                                                                                                      >  2:153809/1‑250 (MQ=255)                                                                                                                                               cAGCGTCACGCGGATGGACACCTTTGCCACTAGAGGAAGCATCACTCTGGACGTGGTTGAAGATGAGCTAAACCGTCTGCGCTATAACTGGCAGGAGCGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAgg                                                                                                   >  2:304017/1‑251 (MQ=255)                                                                                                                                                   gTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCaa                                                                                               <  1:184348/251‑1 (MQ=255)                                                                                                                                                      aCGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGc                                                                                            >  2:140077/1‑251 (MQ=255)                                                                                                                                                                 cACCTTTGCCACTAGCGGACGCATCACTCTGGAAGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATg                                                                                 >  1:106733/1‑251 (MQ=255)                                                                                                                                                                            cTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGATAGCCGCACCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAGACATCGATCTCTTCGACCGTATGCAACTGAAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGCCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCATGCCAAAGCCAGCGTCATGGACGCGGATCgg                                                                      >  1:159985/1‑251 (MQ=255)                                                                                                                                                                                 ggACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACg                                                                 >  1:233576/1‑251 (MQ=255)                                                                                                                                                                                        acACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGCGCCGCTCCTCCGCGCTTACGGCGATGCTGGGCGCAGAGGCAGCAAACATCGATCTCTTCCACCAGATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCCGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGACAATGACGCGGATCGGATACGCACATAc                                                          >  2:260403/2‑251 (MQ=255)                                                                                                                                                                                             cgggACGTGGTGGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGGCGCTTTCCGCAGCCGGACGCCAGCTGTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTgg                                                      <  2:270733/248‑1 (MQ=255)                                                                                                                                                                                                  cGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGtt                                                <  2:217073/251‑1 (MQ=255)                                                                                                                                                                                                      gTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTgg                                            >  2:28744/1‑251 (MQ=255)                                                                                                                                                                                                      gTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAAGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTgg                                            >  2:105372/1‑251 (MQ=255)                                                                                                                                                                                                               gagaTAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCCACTGGAACACGATAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCCGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTg                                   >  1:198707/1‑251 (MQ=255)                                                                                                                                                                                                                    aaaCCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGaa                               <  2:233576/250‑1 (MQ=255)                                                                                                                                                                                                                    aaaCCGTCTGCGCTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGaa                               <  2:272472/250‑1 (MQ=255)                                                                                                                                                                                                                               gcTATAACTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGAt                   >  1:279956/1‑251 (MQ=255)                                                                                                                                                                                                                                    aaCTGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGc               <  1:105309/250‑1 (MQ=255)                                                                                                                                                                                                                                       tGGCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGcac             >  2:218587/1‑249 (MQ=255)                                                                                                                                                                                                                                         gCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAg                                                                                             >  2:135990/1‑167 (MQ=255)                                                                                                                                                                                                                                         gCAGGAGAGCCGCTCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAg                                                                                             <  1:135990/167‑1 (MQ=255)                                                                                                                                                                                                                                           aGGAGAGACGCTCATCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTcc       >  2:199849/1‑251 (MQ=255)                                                                                                                                                                                                                                                 gccgctCCTCCGCGCTTACGGCGTTGCTGGGCGCCGCGGCAGAAAACATCGATCTCTTCGCCCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTCCGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTg  >  2:133655/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      CGAACCTGGATTATGAAGTGGAGCGCCACGCCTCACTCACTGGCGACAGCGTGCGTTTTAACACCGAAGCGCGGCGCGCCTGGTTGGCTTTTGCGACCTCTCCCCAGGCAACATGGCGCGGTAACTTTCGCGAACTTTCTGCCAGCGTCACGCGGATGGCCACCTTTGCCACTAGCGGACGCATCACTCTGGACGTGGTTGAAGATGAGATAAACCGTCTGCGCTATAACTGGCAGGAGAGCCGCCCCTCCGCGCTTACGGCGTTGCTGGGCGCAGAGGCAGAAAACATCGATCTCTTCGACCGTATGCAACTGGAACACGTTAT‑CGCTATCTGCCGCCAGGCAAAGTCGCTTTCCGCAGCCGGACGCCAGCTTTTTGACGTTTCGCGCCAGGGCAAAGCCAGCGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTG  >  NC_000913/3559376‑3559865

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