Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,877,783	C→T	T779T (ACG→ACA)	gyrB ←	DNA gyrase subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,877,783	0	C	T	95.6%	140.3 / NA	45	T779T (ACG→ACA)	gyrB	DNA gyrase subunit B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/14);  minor base A (2/0);  total (31/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

AGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGCGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCACGTTCGATAAACGCAT  >  NC_000913/3877547‑3878009                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   aGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTg                                                                                                                                                                                                                      <  1:73136/251‑1 (MQ=255)     cGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGACCAACTGGTCg                                                                                                                                                                                                                  >  2:19329/1‑251 (MQ=255)                  cGCCCTTTTCTTATAGCTTCTTGCCGGAAGCGGCGTGAAAGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAATTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCCGGCATGGGTAGCGCAATTAAATATCGATATTGGCCGCTTTCAGGGCGGTCTATTCAATAAACGCACGGAGCGGTTCAACGGCGTCGCCCATAAGAGTGGTGAACAACTGGTCGGCAGCGATCGCAt                                                                                                                                                                                                     >  2:159321/1‑251 (MQ=37)                    cccTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCt                                                                                                                                                                                                   <  1:19329/251‑1 (MQ=255)                              ataGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTa                                                                                                                                                                                          >  1:217858/1‑250 (MQ=255)                                     cttGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACgcgc                                                                                                                                                                                    >  1:167026/1‑249 (MQ=255)                                        gCCGGATGAGGCGTGACCCCCTTATCCGGCCTACAAAATCATGAACATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGCACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCAt                                                                                                                                                                               >  2:406029/1‑251 (MQ=255)                                             aTGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTAAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATacgac                                                                                                                                                                          >  1:261568/1‑251 (MQ=255)                                                           ccTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTcc                                                                                                                                                            <  2:365073/251‑1 (MQ=255)                                                                 caggCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCTCATGGTTAGCGCCATTAAATATCGATATGCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGGCGGCAGCAAGCGCATCTTTAACGGTAACGCGCAGCATACGACGACGTTCCGGGTCCATAGt                                                                                                                                                       <  2:181583/248‑1 (MQ=255)                                                                    gCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCAGGTTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCTTTCTCGTCAAGAAACGCGCGGCGCGGTTCAACGGCGTCTCCCGTCAGTGTGGTGAACAAGTGGTCGGCAGCAATCGCAGCTGTAACTGTAACGCGCAGGATACGACGACTTTCCGGGTCCATAGTGGtt                                                                                                                                                   <  1:388670/251‑1 (MQ=255)                                                                     acTACAAAAGCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGAGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATAAGACGACTTTCCGGGTCCATAGTGGttt                                                                                                                                                  >  2:264033/2‑251 (MQ=255)                                                                       tACAAAATCATGAAAATTCAATACATTGCAAGCTTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCCTGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGCACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATACTGGATTcc                                                                                                                                                >  1:364387/1‑251 (MQ=255)                                                                           aaaTCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGCACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCcac                                                                                                                                             >  2:120481/1‑250 (MQ=255)                                                                           aaaTCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCcac                                                                                                                                             >  1:76750/1‑250 (MQ=255)                                                                             aTCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGc                                                                                                                                          >  2:318120/1‑251 (MQ=255)                                                                                     aaTTCCATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACg                                                                                                                                                                                             >  2:263478/1‑192 (MQ=255)                                                                                     aaTTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACg                                                                                                                                                                                             <  1:263478/192‑1 (MQ=255)                                                                                     aaTTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCAGCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCgg                                                                                                                                  <  2:219844/251‑1 (MQ=255)                                                                                           aTACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCAt                                                                                                                            >  1:428718/1‑251 (MQ=255)                                                                                               aTTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCg                                                                                                                        >  1:124444/1‑251 (MQ=255)                                                                                                       aTTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGAcc                                                                                                                >  2:365368/1‑251 (MQ=255)                                                                                                        ttttCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGTGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCt                                                                                                               >  2:279715/1‑251 (MQ=255)                                                                                                             gTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTtata                                                                                                          <  1:318120/251‑1 (MQ=255)                                                                                                                   cTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCtg                                                                                                    >  1:64160/1‑251 (MQ=255)                                                                                                                       tAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATCTCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCtggatg                                                                                                >  1:57807/1‑251 (MQ=255)                                                                                                                            gTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGgag                                                                                             >  1:143854/1‑249 (MQ=255)                                                                                                                            gTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGagag                                                                                           <  1:150171/251‑1 (MQ=255)                                                                                                                                         cACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTACCGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGAc                                                                              >  1:386040/1‑251 (MQ=255)                                                                                                                                            gCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGActct                                                                           >  1:134336/1‑251 (MQ=255)                                                                                                                                               cGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTc                                                                        >  2:370133/1‑251 (MQ=255)                                                                                                                                                             cATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTACCGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTcc                                                            >  1:270611/1‑249 (MQ=255)                                                                                                                                                                       cGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCg                                                <  2:134336/251‑1 (MQ=255)                                                                                                                                                                                    tttCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTAGCCCCTCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCACAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCg                                   >  1:254284/1‑251 (MQ=255)                                                                                                                                                                                        aGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTg                               >  2:160421/1‑251 (MQ=255)                                                                                                                                                                                        aGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTg                               >  2:129275/1‑251 (MQ=255)                                                                                                                                                                                                    tCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTGCCGGGTTCATCTCGCCCAGACCTGTATAACGCGGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCGCGAAGCTGGCTACCGGCTGACGACGCTCg                     <  2:428718/249‑1 (MQ=255)                                                                                                                                                                                                      aaTAAGCGCACGGCACGGTTCAACGGCGTAGCCAATCAGTGTGGTTAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTCACGCGCAGCATACGACGACTTTACGGGTCCATAGTGGTTTCCCACAGATGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGCTGGAGAGGCCGCGACGGGACTCTTTCACCAGACAGTCCAGCGCCTGATCGAAGCTGGCTACCGGCTGACGACGCTCGacac                 >  2:203679/1‑251 (MQ=255)                                                                                                                                                                                                       aTAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTACGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATATCACACAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGCAGCTGGCTACCGGCTGACGCCGCTCGccccg                >  1:272352/1‑248 (MQ=255)                                                                                                                                                                                                                cGGCGCGGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGCTGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGCCGCTCGCCACGTTCGATaa        >  1:383987/1‑250 (MQ=255)                                                                                                                                                                                                                     cgGTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCGGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCACGTTCGATAAACGCAt  <  2:143854/251‑1 (MQ=255)                                                                                                                                                                                                                      ggTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGCCTCTTTCACCAGCCAGTCCAgc                                                         >  2:240552/1‑195 (MQ=255)                                                                                                                                                                                                                      ggTTCAACGGCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAgc                                                         <  1:240552/195‑1 (MQ=255)                                                                                                                                                                                                                              ggCGTCGCTCATCAGTGTGGTGAACAACGGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCGGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCAccagc                                                                  <  2:230353/178‑1 (MQ=255)                                                                                                                                                                                                                              ggCGTCGCCCATCAGTGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCAccagc                                                                  >  1:230353/1‑178 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   AGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAACGGCGTCGCCCATCAGCGTGGTGAACAACTGGTCGGCAGCAATCGCATCTTTAACGGTAACGCGCAGCATACGACGACTTTCCGGGTCCATAGTGGTTTCCCACAGCTGTTCCGGGTTCATCTCGCCCAGACCTTTATAACGCTGGATGGAGAGGCCGCGACGGGACTCTTTCACCAGCCAGTCCAGCGCCTGCTCGAAGCTGGCTACCGGCTGACGACGCTCGCCACGTTCGATAAACGCAT  >  NC_000913/3877547‑3878009

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: ‑