Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,471,154	A→G	15.6%	A58A (GCT→GCC)	tufA ←	translation elongation factor EF‑Tu 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,471,154	0	A	G	15.6%	133.3 / 2.3	67	A58A (GCT→GCC)	tufA	translation elongation factor EF‑Tu 1
Reads supporting (aligned to +/- strand):  ref base A (44/8);  new base G (6/6);  total (50/17)
Fisher's exact test for biased strand distribution p-value = 1.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGTTCAGGAACACGATGATGTACGGAACGCCTACCTGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCT  >  NC_000913/3470917‑3471392                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 tttGTTCAGGAACACGATGATGTACGGAACGCCTACCTGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGG‑GTCGTATTCAACGT‑GAGAAGGGT‑TGATGGTGATACCACGa                                                                                                                                                                                                                                                  >  2:271750/1‑238 (MQ=14)   ttttCAGGAACACGATGATGTACGGAACGCCTCCCTGACGACCCCGCAGGATGTGCTCCCGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCt                                                                                                                                                                                                                                               <  1:573238/237‑1 (MQ=14)      tCAGGAACACGATGATGTACGGAACGCCTACCTGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTCGGCGCGCGTAGTGCCGGGTCGGGGG‑GTCGGCATTAACGT‑GGGAAGTGT‑TGAGGGGGATACCACGGGCtttt                                                                                                                                                                                                                                            >  1:532117/1‑239 (MQ=11)               cGATGATGTACGGAACGCCTACCTGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGC‑GTCGTATTCAACGT‑GAGAAGGGT‑TGTTGGGGATACCACGAACTTTtttttctcgc                                                                                                                                                                                                                                   >  1:129917/1‑235 (MQ=255)               cGATGATGTACGGAACGCCTACCTGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCGGAGCAGCACCGGGGATCATGTTTTTAACATATGCGGCGTGCCCCGGGCAGCCTCCGTGGGCGAAGTGCCGGGTGGGGGG‑GAGGTCTTGACCGG‑GGGAAAGGG‑TGTTGGGGTAACCACGAGGTTTTTTTTCgggg                                                                                                                                                                                                                                   >  1:281776/1‑238 (MQ=11)                    atgTACGGAACGCCTACCTGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGgcg                                                                                                                                                                                                                                  <  2:336606/235‑1 (MQ=31)                            aaCGCCTACCTGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGG‑GCCGTCTTCAACTT‑GGGAAGGGT‑TGATGGTGAAACCACGAGCTTTTTCTTCCgggggggtagtgatcg                                                                                                                                                                                                                      >  1:436140/1‑225 (MQ=11)                            aaCGCCTACCTGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGGGCCCCGGGCAGTCTACGTGTGCGTAGTGCCGGGTCGGGGG‑GTGGGGTTCAACGG‑GGCACGAGT‑TGTGGGGGGTACCCCGGGCttttgtttcgggggggtttctggatg                                                                                                                                                                                                                      >  1:688967/1‑217 (MQ=11)                              cGCCTACCTGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTgg                                                                                                                                                                                                                   >  2:449672/1‑239 (MQ=31)                              cGCCTACCTGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGGGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTgg                                                                                                                                                                                                                   >  2:609311/1‑239 (MQ=14)                              cGCCTACCTGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACAGGTGATCATGTTTTTAACATATTCGGCGAGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGG‑GTGGGATTCAACGG‑GGGAAGGTC‑TGATGGGTATACACCAGGCTTTTTTTTCTGgccggtttaggaa                                                                                                                                                                                                                        >  1:336606/1‑224 (MQ=11)                                 cTACCTGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCg                                                                                                                                                                                                                >  2:713909/1‑239 (MQ=31)                                   cccTGACGGCCCAGGAGGATGTGCTCGCGCGTCTGCGGCATCGGGCCGTCAGTGGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGaa                                                                                                                                                                                                              <  2:839926/238‑1 (MQ=9)                                   accctacgacCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGCCGCAGCAACTACCAGGAGCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGaa                                                                                                                                                                                                              <  2:347222/234‑1 (MQ=14)                                      tGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGACGTGT‑TGATGGGGATACCACGAGCTTTTTTTTCCGGCGCGT‑TATCGATCTGGTCGaagg                                                                                                                                                                                                            >  1:780958/1‑236 (MQ=14)                                        acgacCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCACCGT‑GGAAAGGTT‑TGATGGGGTTAACACAGGCTTTTTCTTCCGGCGCGT‑TTACGATCTGGTCGAATGCAc                                                                                                                                                                                                         >  1:569955/1‑239 (MQ=11)                                               agcagGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGag‑cag                                                                                                                                                                                                  >  2:756722/1‑238 (MQ=32)                                                      tgtgCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCAccgccg                                                                                                                                                                                          >  2:639959/1‑239 (MQ=32)                                                      tgtgCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGGTT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAAGTCCCGGG‑CAGCCccgccg                                                                                                                                                                                          >  2:439333/1‑239 (MQ=255)                                                          cTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAgg                                                                                                                                                                                      >  2:650912/1‑239 (MQ=32)                                                          cTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAgg                                                                                                                                                                                      >  2:170096/1‑239 (MQ=32)                                                            cACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTCTTCAACGT‑GAGAAGTGT‑TGATGGTGATACCGCGGGCTTTTTTTTCCGGCGCGG‑TTTCGGATTGGTCGAATGCACGAA‑CAGCACCCCCGGAGGtt                                                                                                                                                                                    >  1:463190/1‑239 (MQ=11)                                                               gAGTCTGCGGCATCGGGCCGTCAGACGCAGCAACTACCAGGATCGCGCCGTCCCTCTGAGCAGCACCGGTGACCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCAACGTGTGCGTAGTGACGGGTCGGGGG‑GTCGTAGTCACCGT‑GAGAAGGGT‑GGATGGGGATACAACGGTCTTTTTCTTCCGCCGCGT‑TTTGGCTCGGGTGGAATGCA‑GAGACAGACCCGCCGTAGGgttaa                                                                                                                                                                                 >  1:543909/1‑237 (MQ=11)                                                                           tCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GGGAAGGGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAA‑CAGCACCGCCGGGGGTTTTAGCCAGGACgggg                                                                                                                                                                     >  1:24126/1‑237 (MQ=255)                                                                                         cctccatctcccAGGATTGCCCCGTCCTTCTGGGCAGTGCGGGTGATTATTTTTTAAACATGGGGGGGGTCCCCGGGGGTGTCTGCGGGGGCGTAGGTGCGGGTGGGGGG‑GTCGTTTTCAAGGT‑GAGAAGTTT‑TGATGTTGATCCCCCGGCTTTTTTCTCCCGGCGCTT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGt                                                                                                                                                       <  1:183856/229‑1 (MQ=2)                                                                                                 tACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCACGTTTTTAACATAGTCGGCGTGCCCCGAGCCGTCTACGTGTGCGAGGTGACGGGGCGGGGT‑GTAGTATTCAACGT‑GAGCAGTGT‑TGATGGTGATACCACGAGCCTTTTCTTCGGGCGCGT‑TATCGGTCGGGTCGCAGGCACGAG‑CAGCACCGCCGACGGGTTTAGCCAGTCGGGTGGGGGTTggacggggcaggggga                                                                                                                                               >  1:264463/1‑224 (MQ=255)                                                                                                            cgcCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTg                                                                                                                                     >  2:730956/1‑238 (MQ=38)                                                                                                            cgcCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTg                                                                                                                                     >  2:142555/1‑238 (MQ=38)                                                                                                             gcCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCCACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCTGAGGTTTTAGCCAGTACGGGGGTGATTGCAGCGGTCAGAGGAGTTTTAcctgggg                                                                                                                                   >  1:315327/1‑234 (MQ=255)                                                                                                                 tCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAAc                                                                                                                               >  2:753825/1‑239 (MQ=38)                                                                                                                 tCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGAAGTTTTCCCTGGGTCAAc                                                                                                                               >  2:768281/1‑239 (MQ=255)                                                                                                                 tCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCATGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTCCGGGGGGGATTGGAGCGGTCGGAGGAGTTTTACCGTGGTCacg                                                                                                                               >  1:858613/1‑237 (MQ=255)                                                                                                                           cagcaCCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTGGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCCCCGCCGTAGGTTTTAGCCAGTACGGTGGGGATTGCAGCGGCCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATa                                                                                                                     >  2:721168/1‑239 (MQ=255)                                                                                                                           cagcaCCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATa                                                                                                                     >  2:825685/1‑239 (MQ=37)                                                                                                                           cagcaCCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGAGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGCCGCGT‑TATCGATCTGGTCGAATGCACGAC‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGGAGTTTTACCGTGGTCAACGTGGccacaa                                                                                                                     >  2:427048/1‑235 (MQ=255)                                                                                                                                   gTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGGGGTGTTTGCAGCGGTCAAGGGAGTTTTACCGGGGTCAAAGGGGGCGAGAGAaacaac                                                                                                             >  1:400448/1‑239 (MQ=11)                                                                                                                                        cATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGCTTTAGCCAGTACGGTGTGGGTTGCGGCGGTCAGAGTGGTTTTACCGTGGTCAACGGGGCCGATAGTACCCACGTTaa                                                                                                        >  2:253416/1‑239 (MQ=11)                                                                                                                                        cATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGGGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGGGATTGCAGCGGTCAGAGAAGTTTTACCGTGGTAACCGTGGCCGATAGTACCAACGtccg                                                                                                        >  2:192638/1‑236 (MQ=255)                                                                                                                                            ttttgaagtttcggcgggTGCGCGGGGCACTTTACGTGTGGGTGGGTGCGGGGTGGGTT‑GTGTTATTCATGGG‑GAAAATGGT‑TGTGTGTGAAGCCACGAGTTTTTTCTT‑CGGGGTGTGTATGGATTTGGGCGAAGGCAGGAC‑CACCCCCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTTGCCGATAGTACCAACGTTAACGTg                                                                                                    <  1:127885/222‑1 (MQ=255)                                                                                                                                              tttAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACGCCGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTCGCCAGTACGGGGGTGATTGCAGCGGTCAGAGACGTTTTACCGGGGTCACAGGGGACGGAAGGACCAACGTAAAGGTggg                                                                                                  >  1:664024/1‑237 (MQ=11)                                                                                                                                                      aGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCGGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGGACCACCGTTAACGTCCGGTTTTGTa                                                                                          >  2:550087/1‑239 (MQ=255)                                                                                                                                                      aGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGGGATTGCAGCGGTCAGAGTAGTTTTCCCGGGGTCAACGGGGCCGATGGTACCAACGTTAACGTGCGGTTTTGt                                                                                           >  2:4057/1‑238 (MQ=255)                                                                                                                                                        tCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGGGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACg                                                                                        >  2:13693/1‑239 (MQ=38)                                                                                                                                                                     ggCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGGTCGGTTTTGTACGTTCAAATTTTTCt                                                                           >  2:332206/1‑239 (MQ=255)                                                                                                                                                                          tCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGa                                                                      >  2:391225/1‑239 (MQ=40)                                                                                                                                                                          tCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGGGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGGGCGGTTTTGTACGTTCAAATTTTTCTTTAg                                                                       >  2:861783/1‑238 (MQ=255)                                                                                                                                                                            ttggtgtgtgttggggggggTGGGGGT‑GTTTTATTGAAGGT‑GAGGTGTGT‑GGGTGGAGCACCCAGAGGTTTTTTCGTGCGGGGGGT‑TTGCGATTTGGTCAAGTGCACGAG‑AAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGaca                                                                    <  1:817051/222‑1 (MQ=9)                                                                                                                                                                            gaggtgggggTAGTGAGGGGTCGGGGT‑TTTGTATTCACCGT‑GAGAAGTGT‑TGATGGTGATACCACGAGTTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGaca                                                                    <  2:463190/230‑1 (MQ=14)                                                                                                                                                                                  tgCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGGGGCCGATAGTACCAACGTTAACGCGGGGTTTTGTAGGTTCAAATTTTTTTTTAGACACGGGta                                                              >  1:894553/1‑239 (MQ=255)                                                                                                                                                                                    cGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGGTCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGGGTAGTTTTACCTTGGTCACCGTGGCCGAAAGTACCAATGTTAAGGGGCGGGTTTGTGCGTTCAAATTTTTCTTTAGACACGGCaata                                                            >  1:350403/1‑239 (MQ=14)                                                                                                                                                                                    cGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTACCGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCtata                                                            >  2:462071/1‑239 (MQ=255)                                                                                                                                                                                    cGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCttttaaccgggctaaa                                                            >  2:448182/1‑226 (MQ=255)                                                                                                                                                                                    cGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCtata                                                            >  2:679083/1‑239 (MQ=255)                                                                                                                                                                                    cGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTAAACGTGGCCGATAGTACCACCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCtata                                                            >  2:114119/1‑239 (MQ=255)                                                                                                                                                                                    cGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGGTGTCAACGTGCCCGATAGTGCCAACGGTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACCGCtata                                                            >  2:850445/1‑239 (MQ=255)                                                                                                                                                                                     gTAGTGGCGGGTCGGGGG‑TTGTTTTTAAAGGG‑GAGAGGGTT‑TGATGGTTAACCCACGGGTTTTTCCTTCCGGCGCGT‑TACCGATGGGGTCGAATGCAGGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCtata                                                            <  2:894553/238‑1 (MQ=14)                                                                                                                                                                                          gcgggggCGGGGT‑GTGTTTTTCAAGTG‑GAGAGGTGT‑GGATGGTGATCCCCGGGGCTTTTCCTCCCGGCGCGT‑TACGGTTTGGGCCGAATGACCGAG‑CAGCTCCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCt                                                       <  2:858613/235‑1 (MQ=14)                                                                                                                                                                                          gcggggTGGGGGT‑GTGGTTTTCAACGG‑GAGAATTTT‑TGATGGTAATCCCCGGGGTTTTTTCTTCGGGCGCTT‑TATCGATCGGGTCGAATCCACGAG‑CGGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCt                                                       <  2:24126/235‑1 (MQ=14)                                                                                                                                                                                            gcgggtggggggtcttgtttCAAAGTGGAGAAGTTT‑TGATGGTGATCCCACGAGCTTTTTCTCCGGGCGGTT‑TATCGCTTGGGTCGAAGGCACGAG‑AAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCt                                                       <  2:436140/220‑1 (MQ=14)                                                                                                                                                                                               tgcGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATCCCAGGCTCTTTTTCTGCCGGCGCGT‑TATCGATCGGGTCGAATGCACGGG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACta                                                  <  2:129917/236‑1 (MQ=255)                                                                                                                                                                                                gtggggTGGTCGTTTTCAAGGG‑GAGAA‑TGTGTGAGGGTGATACCACGCGTTTTTTCCCCCTGCGCTT‑TATCGATCTGGTCGATTGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACta                                                  <  2:342256/236‑1 (MQ=17)                                                                                                                                                                                                          tgTATACAACGT‑GGAAGGTGT‑TGAGTGTAATCCCCCGGTTTTTTTCTCCCGGCGCGT‑TATCGACGTGGTCGAATGCACGAG‑CCGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCtctc                                        <  2:664024/237‑1 (MQ=17)                                                                                                                                                                                                           ttttttcacgttgggatgttt‑tGATGGTGATCCCGAGCGTTTTTTTTCCGGGCGGTT‑TAGCGCTTGCGGCGAATGCAGGGG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCtctc                                        <  1:730956/218‑1 (MQ=17)                                                                                                                                                                                                                      tagtaggtgt‑TGATGGTGATCCCTCGAGTTTTTTCTTCCGGCGGGT‑TAGCGATTTGGCCGAATGCAGGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGga                              <  2:780958/231‑1 (MQ=18)                                                                                                                                                                                                                       gagaAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagaa                           >  2:490502/1‑238 (MQ=255)                                                                                                                                                                                                                       gagaAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTAATTTTCTTACTATAGGGCTCTCCCCTTCCGgagag                           >  1:756612/1‑239 (MQ=255)                                                                                                                                                                                                                                        gATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGCTAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTAAAGTGCTCTCCCCTTCAGGCGAGAGCACggggattgggt           >  2:771496/1‑231 (MQ=255)                                                                                                                                                                                                                                        gATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTACGAGAAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCAAAATTCCTTACTATAGGGGTCCCCCCTTCAGGGGAAGTCACGGGCCTTggga           >  1:833899/1‑238 (MQ=17)                                                                                                                                                                                                                                                 aGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGGGGCCGCTAGTACCAACGTTAACGTGCGGTTTTGTCCGTTCAAATTTTTCTTTAGCAACGGCTAAATTCCTTACTATAGGGCTTTCCCCTTCAGGGGAGGGACCGGGGTTTTGGTTTTAACCCt  >  1:806493/1‑239 (MQ=17)                                                                                                                                                                                                                                                 aGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGGGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGTTAGTACCACCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTg                                             >  2:860718/1‑196 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 TTTGTTCAGGAACACGATGATGTACGGAACGCCTACCTGACGACCCAGCAGGATGTGCTCACGAGTCTGCGGCATCGGGCCGTCAGTCGCAGCAACTACCAGGATCGCGCCGTCCATCTGAGCAGCACCGGTGATCATGTTTTTAACATAGTCGGCGTGCCCCGGGCAGTCTACGTGTGCGTAGTGACGGGTCGGGGT‑GTCGTATTCAACGT‑GAGAAGTGT‑TGATGGTGATACCACGAGCTTTTTCTTCCGGCGCGT‑TATCGATCTGGTCGAATGCACGAG‑CAGCACCGCCGTAGGTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCT  >  NC_000913/3470917‑3471392

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