Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	609,147	(C)7→8	coding (1083/1113 nt)	insL1 →	IS186 transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	609,140	1	.	C	90.7%	141.4 / 8.8	54	coding (1076/1113 nt)	insL1	IS186 transposase
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base C (24/25);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAACCTGAACTCGCGAAAGCGTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TT‑C‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGG  >  NC_000913/609054‑609434                                                                                           |                                                                                                                                                                                                                                                                                                          acgtgtgctcttccgatctgCGTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAg                                                                                    <  2:757307/281‑1 (MQ=2)   acgtgtgctcttccgatctgTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAg                                                                                    <  2:606060/279‑1 (MQ=2)   aaCCTGAACTCGCGAAAGCGTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAAGCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTaa                                                                                 <  1:1091277/301‑1 (MQ=14)     ccTGAACTCGCGAAAGCGTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTTGCACTTATGGGATTAAt                                                                                >  1:1097737/1‑300 (MQ=11)      cTGAACTCGCGAAAGCGTGGATATTTGCTAATCTACTC‑GCCGCATTTTTCATTGACAACATGATCCAGCCAACGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTTCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACATGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTAAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATcc                                                                              <  2:403810/301‑1 (MQ=255)        aaaCACGCGAAAGCGTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATcccc                                                                            <  1:420162/300‑1 (MQ=17)         aaCTCGCGAAAGCGTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATcccc                                                                            <  2:1150693/300‑1 (MQ=17)          aCTCGCGAAAGCGTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCAGCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTCCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGc                                                                          <  1:285283/301‑1 (MQ=17)          aCTCGCGAAAGCGTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAAGTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGc                                                                          <  1:967361/301‑1 (MQ=17)             cgcgAAAGCGTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCata                                                                       >  2:253278/1‑301 (MQ=18)                  aaGGGTTGATATTTGCTAATCTACTA‑GCCGCATTTTTAATGGACGCCATAATCCAGCGAACGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTGGTTGTGGAGAATAACAAAAATGGTCATCTGGTGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCAGTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCgg                                                                  <  1:610092/301‑1 (MQ=11)                    gCGTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAgagatcggaagagcgtcgtgt                                                                >  1:756446/1‑281 (MQ=2)                     cGTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCACATATCCGGtt                                                                <  2:1153838/300‑1 (MQ=25)                      gTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAgagatcggaagagcgtcgtgtag                                                              >  1:605415/1‑279 (MQ=2)                         gATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtca                                                            <  2:906289/300‑1 (MQ=38)                         agtatTTGCAAATCT‑CTCGGCCGCATTTGTAATTGCGGTCAAATACCAGCCACCGCTGG‑ATGTC‑C‑C‑CCC‑CCGGAAGTGCCGGATCCGAAAAGAAGAACTATCTGGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtca                                                            <  1:399577/298‑1 (MQ=11)                             ctcaccggacgcgcttccgatctttgtAATTGACGACAGAATCCAGCCAGCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGAGCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtcaggtc                                                        <  1:848939/279‑1 (MQ=14)                               tGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGtcaggtcagg                                                     >  2:592701/1‑301 (MQ=255)                                 cccgatctcTC‑GCCGCATTTTTAATTGACGACATAAACCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGAt                                                   <  1:1095009/293‑1 (MQ=255)                                 cTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGAt                                                   >  2:1063661/1‑301 (MQ=255)                                 cTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGAt                                                   <  1:964833/301‑1 (MQ=255)                                  tAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTCGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATg                                                  >  1:755646/1‑301 (MQ=255)                                      cTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCGGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTaa                                              <  1:1194619/301‑1 (MQ=255)                                         cTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGAtagatcgga                                           >  2:1096313/1‑293 (MQ=255)                                           c‑gccgcATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCa                                          >  2:1004586/1‑300 (MQ=255)                                               cgcATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCTGAAGGGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTg                                      <  2:662825/301‑1 (MQ=255)                                               cgcATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTg                                      <  2:1231321/301‑1 (MQ=255)                                                 cATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTGAGGATGGGAAGGCACCTg                                      >  1:487510/1‑299 (MQ=255)                                                 cATTATTAATTGACGACAAAATCAAGGCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGc                                     <  1:176125/300‑1 (MQ=255)                                                  aTTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TT‑CAC‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACGCGTTGTGGAGAATAACTAGAATGGTCATCTGGAGCTGACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTAAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCTCCGCATATCCGGTTGTAAGATCAGGATGGTAAGGCACCTGCtt                                   >  2:1078446/1‑301 (MQ=255)                                                   tttttAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCttt                                  <  2:521137/301‑1 (MQ=255)                                                    ttttAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtcaggtcagatcggaagagcacacgtctga                                 >  2:849926/1‑280 (MQ=255)                                                    ttttAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTa                                 >  2:467428/1‑301 (MQ=255)                                                     tttAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTac                                <  1:466922/301‑1 (MQ=255)                                                       tAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTacac                              <  2:263171/301‑1 (MQ=255)                                                        aaTTGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTaca                               >  2:909829/1‑299 (MQ=255)                                                           tGACGACATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGAAAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTACCGCTTATGGGAGTAATCCCCGCATATCCTGTGTCCAGGTCAGGCTGGTAAGGCACCTGCTTTACACttt                           >  1:1152537/1‑300 (MQ=255)                                                             acgacATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCg                         <  1:1062399/300‑1 (MQ=255)                                                             acgacATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCg                         >  2:1043815/1‑300 (MQ=255)                                                             acgacATAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCGTACAGGTGGCCAGTCGTGGGACAGTATCCCTGACAGCCTACAAGAAGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCAACGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGc                        <  2:456193/301‑1 (MQ=255)                                                              cgacAGAATCAATCCAACGCTGG‑AT‑TT‑C‑CGTCT‑CCAGAAATGCCGGATCCGAAAAGAAGAACTAACGCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTTAAGACCGCGAGGCATCGTCTTAACGAGGCACCGATGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGcc                       <  2:1076523/301‑1 (MQ=17)                                                                  gaaaTCCAGCGATGGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTg                   >  1:71326/3‑301 (MQ=255)                                                                  ctAATCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGACGGATCCGAAAAAAAGAACTAACTCGTAGGGGAGATGAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTg                   <  2:1142808/300‑1 (MQ=255)                                                                  atcaTCCAGTCACCTCTAG‑AG‑TT‑C‑C‑CCCACCAGAAAAGCCGGATCCGAGAAGAAAAACTAACTCGTTGTGGAGAATAACAAAAATGGGAAGCGGGAGCTGACAGGTGACCATTCATGGGACAGTATCCCTGACAGCCTACAAAACGGAATAGAAGAACGCGAGGCATCGGCTTAGCGATGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCGTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTg                   <  1:813427/298‑1 (MQ=17)                                                                    aaTCCAGCGATGGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGCTGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGGTGGTAAGGCACATGTTTTACACTTTCGCCTGTGGt                 >  1:846807/1‑301 (MQ=255)                                                                    aaTCCAGCGATGGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGCTGGTAAGGCACCTGCTTTACACTTTCGCCAGTGGt                 >  1:139695/1‑301 (MQ=255)                                                                    aaTCCAGCGATGGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATAACCGGTTGGCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTggg                 >  2:727044/1‑300 (MQ=255)                                                                     aTCCAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTTGCCCGTGGt                 >  1:865897/1‑300 (MQ=255)                                                                       ccAGCCATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACGAACGCGTTGTGGAGAAGAGCAAAAAAGGTCATCTGCAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTACTCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAg              <  1:230235/301‑1 (MQ=255)                                                                            cATCGCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATg         >  2:747159/1‑301 (MQ=255)                                                                            aagcGCTGGAAT‑TC‑C‑C‑CCC‑CCAGAAGAGCCGGATACGAAAAGAAGAACTAACTCGTGGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCAGTCGTGGGACAGTATCCCTGACAGACTACAAAACGCAAGTGAAGAACGCGAGGCATCGTCGTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATg         <  1:1164536/298‑1 (MQ=18)                                                                                gCTGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTCCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATCACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGccg     >  2:11891/1‑299 (MQ=255)                                                                                  tGG‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGACATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAAATGGTTCAACCATTAAGTTAGCGCTTATGGTAGTCATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTATGGCCCCTGCTTTACACTTTCGCACGTCGTCAGTGATGGCTgtg   >  1:230810/1‑299 (MQ=17)                                                                                   gg‑AT‑TTCC‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCgg  >  2:1194210/1‑301 (MQ=255)                                                                                           |                                                                                                                                                                                                                                                                                                          GGAACCTGAACTCGCGAAAGCGTGGATATTTGCTAATCTACTC‑GCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGG‑AT‑TT‑C‑C‑CCC‑CCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGG  >  NC_000913/609054‑609434

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑