Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,504,667	T→G	19.1%	C317W (TGT→TGG)	insQ →	IS609 transposase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,504,667	0	T	G	19.1%	114.9 / 3.4	58	C317W (TGT→TGG)	insQ	IS609 transposase B
Reads supporting (aligned to +/- strand):  ref base T (29/17);  new base G (3/9);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 2.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGG  >  NC_000913/1504436‑1504901                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            ccACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATa                                                                                                                                                                                                                                      <  2:436368/239‑1 (MQ=255)    cGCAATGATAGTCATGGAGGATTTGAAGGGCGGCAACATGGCAAAGTCAGAAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTGTAAA‑CCGTTCGATACTGGATCAGGGCTGGTAGGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAg                                                                                                                                                                                                                                   <  1:1688504/239‑1 (MQ=255)    cGCAATGATAGACATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAg                                                                                                                                                                                                                                   <  1:167840/239‑1 (MQ=255)    cGCAATGAGAGTCATTAAGGATTAGAAGGTCAGGAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAg                                                                                                                                                                                                                                   <  1:798640/239‑1 (MQ=255)                   tGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACCAGCCAG‑CGTTGCGCGTACTTTGGGCATACAGCGAAAGAGAACCGc                                                                                                                                                                                                                     >  1:879037/1‑238 (MQ=255)                          ttGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGGCCTACAGCGAAAGAGAACCGCCTGTcaca                                                                                                                                                                                                             >  2:984881/1‑239 (MQ=255)                                gTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTaa                                                                                                                                                                                                       >  2:891234/1‑239 (MQ=255)                                       atacggcaaaatgaagcacGGGTGCGG‑GCAGCCAGGCGGGGGGCAAGTGCCGGGCAA‑AAATAAGGTTGAAC‑CGGTACGATAGTGGATCAGGGGTGGTTGAAAAGGCGCCCCAAGTGTGCGTAAAAGATGCTTGGGGGTGGGGGTGGGGTGTTTGCGGTCCCGCAAGCGTATAAAAGCCGG‑CGTTGCGGGTACTGTGGTCATACAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAg                                                                                                                                                                                                 <  1:629313/221‑1 (MQ=255)                                        cATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGGTCCGCCAGCGTATACAAGCCAG‑CGTTGTGCGTACTGGGGTCATACAGCGAAAGAGAAACGCCTGTCACAAAAGAAATTCAGAt                                                                                                                                                                                               >  2:912330/1‑239 (MQ=255)                                          tGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGTGGCGTGCTGTTCCGCCCGCGTATACAAGCCAG‑CGTTGCGCGTACTGGGGTCATCAAGCAAAAGGGAAACGCCTGTCCCCAAGAAAAATCAAATg                                                                                                                                                                                              >  1:721603/1‑238 (MQ=255)                                            aaaaagtagggggggggTACGG‑TCAGCAGGCGGGGGCGAAATGCCGGGGA‑A‑AAATAAGGTTAAAC‑CGGTGGGAAACTGGTTCGGGGTGGGTAGAAAATGCGCCGCCGGTTGGCGTAAAAGCAGCTGGGGGGGGGGGGTCAGGTGCTGGCTGTCCCGCAGGGGAAAAAAAGCCGG‑CGTTGCGCGTACGGGGGTAATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGc                                                                                                                                                                                             <  1:1248227/223‑1 (MQ=255)                                             cAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCTTGGCGGTCAGGGGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑GTTTGCGCGTACTGGGGTCATACAGCGAAAGGGAACCGCCTGTCCCAAAGTAAATTCAGATGCCAg                                                                                                                                                                                          >  2:735576/1‑239 (MQ=255)                                              aaaGGCAGCGGCGGGTAGGG‑TCAGCCGGCCGGGGGGAAAGGGCGGGGC‑A‑AAAGCGGGTTTAAA‑CGGTTGGATACGGGATCGGGGCGGGGATGAAGTGCGCCGCCAGCTGGCGTATAAGCAGCTCGGGGGGGGGGGTGGGGTGCTTGCTGTTCCGCCAGCATATACAAGCCGG‑TGTTGCGAGTATTGTGTTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAgg                                                                                                                                                                                         <  1:1562383/239‑1 (MQ=255)                                               aaagtcggcgcgggttcgggcccgccaccgGGGGGGAAATG‑CCGGGG‑A‑AAAACCGGTTTAAACCCGTTCAATACTGGACCGGGGCTGGTAGGAACGGCGCCGCCAGCTTGCGTAAAAGCAGCTCGGGCGTGGCGGTCGGGTGCTTGCTGTCCCGCCAGCGTATACAAGCCAG‑CTTTGCGCTTACTGTGGTCATACAGCGAAAGAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAgg                                                                                                                                                                                         <  1:346192/216‑1 (MQ=255)                                                aGTCAGGAGGGGGTACGG‑TCAGCCCGCCGGGGCGCAGTGTCGGGGA‑A‑AATACGGGTTTAAC‑CGGTGCGATTGGGGATCGGGGCGGGTAGGAAATGCGCCGCCGGTTGGCGTATAAGCAGCTGTGGCGGGGGGGTCGGGTGTTGGCTGTCCCGCCAGGGTATACAAGCCGG‑CGTGGCGCGTATTGGGGTCATACAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTa                                                                                                                                                                                       <  1:472553/239‑1 (MQ=255)                                                    cagcgggggTACGGCACAGCGCGCGGGGGCGCATTG‑CCGGGC‑ACAAAGCGGTTTAAAA‑CCTTTCGAGACGGGACGAGGGGTGGTATGAAGTGGGCCGCCGGCTTGCGTATAAGAGGCTGTGGCTTGGGGGTCGGGTGTTTGCTTTTCCGCCAGCGTATACAAGCCAG‑TGTTGCGCGTACTGTGGTAATAAAGCGAAAGAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTAtg                                                                                                                                                                                     <  1:657711/232‑1 (MQ=255)                                                       cagcggtgacggcgaccagcccGGGGCCCAAGTCCCGGAA‑A‑AAATCAGGTTAAAC‑CAGTCCAATCCGGGATCGGGGAGGGTATAAAAGGCGCCACCAGCTTGCGTATAAGCAGCTCGGGCGTGGCGTTCCGGTTTTTGCTGTCCGACCAGCGTAAAAAAGCCAG‑CGTTGCGAGTACTGTGGTCATACAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTgg                                                                                                                                                                                  <  1:248882/218‑1 (MQ=255)                                                                    acgccgcccGGGGAGCAATGGCCGGAC‑A‑AATCCGGGTTTAAC‑TCGTTCGATGGTGGAGTGGGGTGGTTATGAAGTCCGCCGCAGGCTTTCTTTAAAGAACCTCTGGCGTGGCGGTAGGGGGCTGGCTGTCCCGCCAGGGTAAACAAGCAAG‑CTTTGGGCGATTTGTGGTCACACAGGAAAAAAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACg                                                                                                                                                                    <  1:1075939/232‑1 (MQ=255)                                                                      gccagccGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGGCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGAAGCCCGGTATGTGGATATACAGCGAACGcc                                                                                                                                                                  >  1:1491885/1‑239 (MQ=255)                                                                           ccGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGTAAGAAAACCGCCTGGCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCCAACGCCGATGt                                                                                                                                                             >  2:841948/1‑239 (MQ=255)                                                                               ggCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGGAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAAt                                                                                                                                                         >  2:864979/1‑239 (MQ=255)                                                                                 cgcAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCTGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGGCATACAGCGAAAGAGAACCGCCTGTTACAAAGTAAATTCAGATGCCAGGTATGTGGATCTACAGCGAACGCCGATGTAAATg                                                                                                                                                        >  2:171491/1‑238 (MQ=255)                                                                                    aaTGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCGGTTCCGCCAGCGTATACAAGCCCG‑CGTTGCGCGGACGGTGGTCCTACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGgcg                                                                                                                                                     >  1:698552/1‑238 (MQ=255)                                                                                          cGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGGGCTTGCTGTTCCGCCAGCGTATACAAGACAG‑CGTTGCGCGTACTGTGGTTATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCCGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCa                                                                                                                                               >  2:9325/1‑238 (MQ=255)                                                                                                                    ttCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATa                                                                                                                                                                                                                                      >  1:398693/1‑128 (MQ=255)                                                                                                                    ttCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATa                                                                                                                                                                                                                                      <  2:398693/128‑1 (MQ=255)                                                                                                                    ttCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGGACGcc                                                                                                                       >  2:391900/1‑238 (MQ=255)                                                                                                                    ttCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGGGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGGAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAAAATTTTAGCGGC‑GGGGCACGcc                                                                                                                       >  1:845539/1‑238 (MQ=255)                                                                                                                            tGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGGACCCCgggccttc                                                                                                               >  2:167746/1‑231 (MQ=255)                                                                                                                              gATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGGGGATATACAGCGAACCCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCGCGCCGCTCTTGc                                                                                                               >  1:780420/1‑236 (MQ=255)                                                                                                                                   gggCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATCCAGCGAAAGAGAAACGCCTTTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGCACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGccggcccctactgctcgtggag                                                                                                       >  1:404728/1‑218 (MQ=255)                                                                                                                                   gggCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑gggggccacccgtctttcctctgggg                                                                                                       >  2:107513/1‑217 (MQ=255)                                                                                                                                   gggCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGcgcccccctcttttccttgggg                                                                                                       >  1:629078/1‑218 (MQ=255)                                                                                                                                   gggCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGCAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGGGCGccccgctttcccggttga                                                                                                       >  1:1296763/1‑223 (MQ=255)                                                                                                                                   gggCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGCGCTTGCTGTTCCGCCAACGTATACCAACCAG‑CGTTGCGCGTACTGTGGTCATACAGCGCAAGAGAACCGCCTGTCCCAAAGGAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGGTGTAAAAGGCGCTCCGAACATTTTAGCGGC‑CGGGCACGCCGCTCttcgcctcggag                                                                                                       >  1:995808/1‑229 (MQ=255)                                                                                                                                           agaaaaTGCGCCGCCGGTTGGGGTATAAGCAGCTCGGGGGGGGGGGTCGGGTGGTGGCTGTCCCGCAGGCGAAAAAAAGCCAG‑GGTTGCGCGGACGGGGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTg                                                                                               <  1:107513/236‑1 (MQ=255)                                                                                                                                                   cgccgcCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGGCCGCCGTGCTTTCCTGTGgggagagagggcggtcgcg                                                                                       >  2:1657660/1‑226 (MQ=255)                                                                                                                                                   cgccgcCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGGACGCCGCTCTTTGCTCTGGTGAGGCGGgtcagtcaa                                                                                        >  1:166238/1‑237 (MQ=255)                                                                                                                                                   cgccgcCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCCTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCCCAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACCTTTTAGCGGCGGGGGCCCCCCCTCCTTGCCTGTgggggggggggggggtcag                                                                                        >  1:962587/1‑222 (MQ=255)                                                                                                                                                        ccAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGGTCTTGCCCGTGGAGAGGTGGGGGAGACAAgccgc                                                                                   >  2:1722442/1‑238 (MQ=255)                                                                                                                                                                     aaGCAGCTCTGGCGTGGCGGTCGGGTGCTTGCTGTTCCGCCAGCTTATAAAAGCCAG‑GGTTGCGCGTACTGGGGTCATCCAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGa                                                                     <  2:1201933/239‑1 (MQ=255)                                                                                                                                                                      agcggcTCGGGCGGGGCGTTAGGGTTCTGGCTGTCCCGCAGGCGTATAAAAGCCGG‑CGTGGGCCGTACGGGGGTCAAACGGCGAAAAAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGaa                                                                    <  2:1336824/239‑1 (MQ=255)                                                                                                                                                                                gcgtggcgGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATg                                                                                                                                                              >  2:854833/1‑140 (MQ=255)                                                                                                                                                                                     gcgGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGGCCGCCGGTCTTTTCTGGGGGGGGATGGTGTGGTCAGGGCGCCCGCTGGAGCAGGGAcccccccccaaagtt                                                     >  2:1108176/1‑230 (MQ=255)                                                                                                                                                                                             gtgtTTGCTGTCCCGCCAGCGTATACAAGCCAG‑AGTTGCGCGTACTGGGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGa                                              <  2:857263/238‑1 (MQ=255)                                                                                                                                                                                             gTGCTGGCGGTCCCGCCGGCGTAAACAAGCCGG‑CGTGGCGCTTACTGTGGCCATACAGCAAAAAAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGa                                              <  2:962587/238‑1 (MQ=255)                                                                                                                                                                                              tgcttg‑tgTTCCCCCAGTGTTTACAAGCCCGATGTTGCGGGACCTGGGGCCATAGACAGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCg                                               <  2:177260/236‑1 (MQ=255)                                                                                                                                                                                                  tcgttttCCGCCAGCGTATACAAGCAAG‑CGTTGCGCGTATTGTGGTCAAACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCt                                        <  2:1717341/234‑1 (MQ=255)                                                                                                                                                                                                     tGTCCCGCCGGCGTATAAAAGCCCG‑GTTTGCGGGTATTGTGGTAATACAGCAAAAGAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGaa                                     <  2:1718510/239‑1 (MQ=255)                                                                                                                                                                                                     tGTCCCGCCAGCGTATCCAAGCGAG‑GGTTGCGCGTATTGGGGTAATACAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGaa                                     <  2:1111438/239‑1 (MQ=255)                                                                                                                                                                                                      taataggatgaagataAAAGCCGG‑GGTTGGGCGAACGGTGGCAAAACAGCAAAGGAGACCCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGaa                                     <  2:1419784/225‑1 (MQ=255)                                                                                                                                                                                                                   aaaCAAGCCGG‑GGTTGCGCGTCCTGTGTACAAACGGGGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGACCCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATcc                       <  1:1720600/237‑1 (MQ=255)                                                                                                                                                                                                                     aCAAGCCAG‑CGTTGCGCGTACTGTGGTCATAAAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCAc                     <  2:879037/239‑1 (MQ=255)                                                                                                                                                                                                                       aaGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTca                                                                                                                                                                                                               >  2:918180/1‑52 (MQ=255)                                                                                                                                                                                                                          ccAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGGTCTTTCCTTTGGAGAGATGGGGCAGGCACGCCGCCCGCTGGAGCAGGAACCCCCCCAAATAATTCAAGCGGCCGCCTGAAAGCAGCAGGGCTCCCCGTCCt                >  2:674935/1‑239 (MQ=255)                                                                                                                                                                                                                          ccAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCCGTCTTGCCTGTGGGGAGATGGGGCCGTCAAGCCGCCCGCTTAAGCAGGAACCCCCCCAAATGATTTCGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCt                >  2:173481/1‑239 (MQ=255)                                                                                                                                                                                                                               cGTTGCGCTTCCTGGGGACTTACAGGAAAGGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCa             <  2:721603/238‑1 (MQ=255)                                                                                                                                                                                                                                     gcgaattgtgGTCATACAGAGAAAGAGACCCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTg      <  1:167746/233‑1 (MQ=255)                                                                                                                                                                                                                                         aCTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGGTCTTGCCTGTGGAGAGATGGGGCAGGCAGGCCCCCCGTTGCAGCAGGAACCCACCGAAATGATTCAGGGGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTggggg  >  2:1498591/1‑239 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            CCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGC‑A‑AAATCAGGTTTAAA‑CCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGC‑GGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGG  >  NC_000913/1504436‑1504901

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