Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,340,431	C→T	18.6%	F63F (TTC→TTT)	yciS →	DUF1049 family inner membrane protein, function unknown

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,340,431	0	C	T	18.6%	143.5 / 5.7	70	F63F (TTC→TTT)	yciS	DUF1049 family inner membrane protein, function unknown
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base T (12/2);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 3.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TAAGTTGACCATAATTTATTCGCTCTAACCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCAAGATGAAGCCAACCGCTTGTCGCG  >  NC_000913/1340196‑1340667                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  tAAGTTGACCATAATTTATTCGCTCTAACCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTTTTG‑TTTTCTTCGGGGTTTGCTTTTGGTTGGTTGCTTTGTTGCCTGTTCTgg                                                                                                                                                                                                                                                >  1:704366/1‑239 (MQ=255)  aaGTTGACCATAATTTATTCGCTCTAACCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTgg                                                                                                                                                                                                                                                <  1:91376/238‑1 (MQ=255)  aaGTTGACCATAATTTATTCGCTCTAACCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGGATTT‑TTTGCTGCGGGGTTTTCTATCGGTTGGTTGTTTTTTGGCCTGTTTTggg                                                                                                                                                                                                                                               >  2:955348/1‑238 (MQ=255)        aCCATAATTTATTCGCTCTAACCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTTCTGCGGGGGTTGCTATCGGGTGGGTGATTTTTGGCCTGTTCTGGTTGCggg                                                                                                                                                                                                                                         >  1:586820/1‑237 (MQ=255)              aTTTATTCGCTTTTACCCAATAACGGGAAGTAATTTGAAATTTTTACTCATTTTCTTGGTGGTGTTAGCG‑ATCTTCGTTATTTCTGTTACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCgtg                                                                                                                                                                                                                                   <  2:989084/239‑1 (MQ=255)                  aTTCGCTCTAACCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGGTTGCTATCGGTTGGTTGATTTTTGGCCGGTTTTGGCTGCGGGTTCGGG‑ttt                                                                                                                                                                                                                               >  1:876897/1‑238 (MQ=255)                   ttCGCTCTAACCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTcc                                                                                                                                                                                                                             <  1:875556/239‑1 (MQ=255)                    tCGCTCTAACCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTT‑TTTTCTGCGGGGGTTTCTATCGGGTGGTTTATTTTGGGCCTGTTTTTGCTggggggtctggtttt                                                                                                                                                                                                                               >  2:626893/1‑223 (MQ=255)                    tCGCTCTAACCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTTCGGGGTTTGCTATCGGTTGGTTGAATTGGGGCCGGTTCTTGCTGCGAGTTCGTG‑TTTc                                                                                                                                                                                                                              >  1:316319/1‑237 (MQ=255)                       ctctAACCACATAACGGGAAGTAATGTGACATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCCCATTGCTGGCGGGCTTG‑TTTGCTGCGGGGTTTTCTATCGGTTGGTTGATTTGTGGCCTGTTTTTGTTGCGAGTTCTTG‑ttttctctg                                                                                                                                                                                                                         >  2:739335/1‑233 (MQ=255)                       ctctAACCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCCCATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTTTCGGTTGGTTGGTTTCTGGCCTGTTTTGGCTCGGGGTTCGTG‑TTTCCCTgt                                                                                                                                                                                                                         >  1:445256/1‑238 (MQ=255)                       ctctAACCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGGTTTGTATCGGGTGTGTGTTTTTTGGCCGGTTTTGGCTGCGGGtttggggttgccctg                                                                                                                                                                                                                          >  1:610510/1‑226 (MQ=255)                        tctAACCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGgc                                                                                                                                                                                                                        >  1:60138/1‑239 (MQ=255)                          tAACCACATAACGGGAAGTAATGTGAAATATTTACTCCTTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGCATTG‑TTTGCTCCGGGGTTGGCGATCGGTTGGTTGATTTTTGGGCTGTTTTGGCTGAGAGTTCCTG‑TTTTCCTgttgc                                                                                                                                                                                                                      >  1:565312/1‑235 (MQ=255)                            aCCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTTGCCTGTTCTGGTTGCGAGTTCGGG‑TTTCCCTGgggcgc                                                                                                                                                                                                                    >  1:965150/1‑239 (MQ=255)                                     cGGGAAGTAATGTGAAATTTTTACCTATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGt                                                                                                                                                                                                           <  1:686877/239‑1 (MQ=255)                                          aGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGGATTT‑TTTGCTTGGGGGTTTGCTTTCGGTTGGGTGTTTTGGGGCCTGTTTTGGGTGCGGGGTTGTG‑TTTCCCCGGGGGGCGCGGGACGAaaaaa                                                                                                                                                                                                      >  2:1006548/1‑238 (MQ=255)                                          aGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGGATTGTTTTACTG‑GGGGTTTGCTATCGGTTTGTTGTTTTTTGGCCTGTTTTGGCTGTGGGGTTCTG‑TTTTCCTgggcgggcgtggaagcgaaaa                                                                                                                                                                                                      >  2:332205/1‑220 (MQ=255)                                               tgtgTAATATTTACTTATTTTCTTTCTGGGGTTAGCG‑ATCTTGGTGATTTCGGTTACGGTGGGTGCGCAGAACGATCAACAGGTGAAGTTTAATTAGCTGTTAGCGGAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGc                                                                                                                                                                                                <  2:704366/239‑1 (MQ=255)                                                      tatTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTT‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTTGGGCCTGTTTTGTGTGCGGGTTCGGG‑TTTCCCCGGCGGGCCCTGAAAGTaaaaa‑aaaagaactggg                                                                                                                                                                                         >  2:491794/1‑230 (MQ=255)                                                                 tttCTTACTGGTGTTAGCGTATC‑TGGGGATTCCGGTCACGTGGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCttt                                                                                                                                                                             <  1:470446/239‑1 (MQ=255)                                                                       aCTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTTCTTCGGGGTTTGCTATCGGTTGGTTGATTTTTTGCCTGTTCTTGCTGCGGGTTCGTGTTTTCCCGGGGGCG‑GCCGGACGGAAAATAAAAGCG‑CCGGAAAA‑ACAGATTTTCCCcg                                                                                                                                                                       >  2:963844/1‑239 (MQ=255)                                                                           gtgtTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGTTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGCCGGTATTG‑TTTGGTTGGTGGTTTGCTATCGGTTGGTTGGTTTGGGGCCTGGTCTGGCTTCCAGTTCGTG‑TTTTCCTGGCGCGCGTTTACCGCGTAAA‑AAAACGGCTGGaaaaaaagacctacacaccgca                                                                                                                                                                    >  2:20207/1‑221 (MQ=255)                                                                            ggtTAGCG‑ATCTTCGGTATTTCGGTCACGGTTGGTGCGCAGAACGATCAACAGGTTACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGAc                                                                                                                                                                   <  1:1006548/237‑1 (MQ=255)                                                                               tttggcatttttgtgtTTTTGGTCCCGTTGGGTGCGCCGAACGTACAACAGGTGACGTTTGATTATCTGTTAGCCGAAGGGGAGTTCCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCCGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGa                                                                                                                                                                <  2:269807/227‑1 (MQ=255)                                                                               tAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGCATTG‑TTTGCTGCGGGGGTTGCTTTCGGGTGGTTGATTTTTGGCCTGTTTTTGCTGCGGGTTCGGG‑TTTCCCTGGCGCGCGCTGGACGGCAAAT‑AAAGCGGCGTGAAAA‑CACACTTTCACCCCCGCCTggg                                                                                                                                                               >  2:1054293/1‑237 (MQ=255)                                                                                       cTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGGTCGTG‑TTTCCCTGGCGCGCGCCGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTg                                                                                                                                                            >  1:738160/1‑235 (MQ=255)                                                                                       cTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGGTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGGGTTCGGG‑TTTCCCTGGCGCCCGCTGAACTGAAAAT‑AAAGCGACTGGAAAA‑CCAGCttttcccccggcccgactgggcgg                                                                                                                                                        >  2:1038285/1‑218 (MQ=255)                                                                                       cTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGAATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGCCTGACGTGGCTg                                                                                                                                                        >  1:342579/1‑239 (MQ=255)                                                                                          cGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTTGGGCCTGTTCTTGCTGCGGGTTCGGG‑TTTCCCCGGGGGGCGCTGGACGTAAAAA‑AAACGGACGGGGAAA‑CCAGCTTTCccccccccccggcggggggggag                                                                                                                                                     >  2:808188/1‑220 (MQ=255)                                                                                            ttttttCGGTTCCGTTGGGGGCGAAGAACGATCAACAGGGGACGTTTAATTTTCTTGTAGCGGAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTg                                                                                                                                                   <  2:708966/236‑1 (MQ=255)                                                                                            ttttttCGGTCACTGTGGGGGCGCAGAACCATCAACAGGGTACGTTTAATTATCTGTTAGCTCAAGGGGGGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTg                                                                                                                                                   <  2:316319/236‑1 (MQ=255)                                                                                                 tCGGTTGCGTTGGGGGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTTCCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGGGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCa                                                                                                                                             <  1:955348/239‑1 (MQ=255)                                                                                                 tCGGTCTAGTTGGGTGCGAAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCa                                                                                                                                             <  2:977283/239‑1 (MQ=255)                                                                                                 tCGGTCAGGTGGGGTGCGCAGCACTATCAACAGGTGACGTTTAATTTTTTGTTTGCGCAAGGGGAGTACCGTATTGCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCa                                                                                                                                             <  1:1054293/239‑1 (MQ=255)                                                                                                                       aCGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGGG‑TTTCCCTGGCGGGCGCTGAACGTAAAAT‑AAAACGCCTGGAAAA‑CCAGCTTTCACCCGCGACTGACGGGGCTGTAGGGC‑CGCACTCGTCAACGGCGAAGGAATaa                                                                                                                       >  1:729241/1‑239 (MQ=255)                                                                                                                           tCAACAGGTGGCGTTTAATTTGCTTTTATGGCAAGGGGTGTACCGTATTGCCAGATTGGTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACttt                                                                                                                   <  2:504584/239‑1 (MQ=255)                                                                                                                                aGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGCCTTACGTGGCTGTAGTGC‑CGCACTCGGCAGCGGCGAAGGAATAAAtttttttg                                                                                                              >  1:1018720/1‑234 (MQ=255)                                                                                                                                          tAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGTTGCGGGGGTTGCTATCGGGTGGGTGTTTTGTGGCCTGTTTTGGCTTGGGGTTTGTG‑TTTCCCTGGCGCGGGCTGAACGTAAAAT‑AACGCGCCTGGAGAA‑ACAGCGTTTACCCGCGACCGACGTGGCTGTAGTGC‑GGCACCCGCCAGCGGCGCAGGGAAAACTTTTTATggggggggttg                                                                                                    >  2:687514/1‑229 (MQ=255)                                                                                                                                                 cTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGgcgcgcg                                                                                                                                                                                                                   >  2:457597/1‑124 (MQ=255)                                                                                                                                                 cTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGgcgcgcg                                                                                                                                                                                                                   <  1:457597/124‑1 (MQ=255)                                                                                                                                                    ttAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGGTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCCGCGGGTTCGTG‑TTTTCCTGGCGCGCGCTGCACGGAAAAT‑AAAGCGCCCGGGAAA‑CCAGCTTTCACCCGCGCCTTCCCGGGCGGTAGGGC‑CCCACCCGCCATCGGCGAAGGCAAAACGTTCTTTCGT‑GGGGGTGtttttttttgt                                                                                         >  2:704019/1‑234 (MQ=255)                                                                                                                                                                  gTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTTGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGGCTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCGTTTGCCTGTTGcc                                                                            >  1:192651/1‑238 (MQ=255)                                                                                                                                                                   ttccgttttttcgattggCTTGCGGAAATG‑GTTGCTGGGGGGTTTTCTTTCGGGTTGTGTGTTGTTGTCCTGTTCGGGCTGCGAGTTTGTG‑TTTCCCTGGCGCGCGCTGAACGGAAAAT‑AAAGCGACGGGAAGA‑CTAGCTTTCACCCGCGACTGACGTAGCTGTAATGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCGTTTGCCTGTAGCCg                                                                           <  2:445256/221‑1 (MQ=255)                                                                                                                                                                   ttccGTTTTTCCACATTGTTGGCGGTTTTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCGGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCg                                                                           <  2:646416/236‑1 (MQ=255)                                                                                                                                                                   gtccGTATTTCCTCATTTCTGGCGTTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCg                                                                           <  1:381379/236‑1 (MQ=255)                                                                                                                                                                    caggTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCg                                                                           <  1:355819/234‑1 (MQ=255)                                                                                                                                                                      cGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCa                                                                                                                                             <  1:328350/170‑1 (MQ=255)                                                                                                                                                                         ttttCCCCTTTGGTGGCGGTTTTG‑TTTGGTGGGGGGTTTGCTATGGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTa                                                                    <  2:586820/238‑1 (MQ=255)                                                                                                                                                                          tttCCACATTGCTGGCGGCATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGTCTGCGAGGTCGTGTTTTCCCTGGCGCGCGCTGGACGGAAAAT‑AAAGCGACTGGAAAACCCAGC‑TTCCCCCGCGCCTTACGTCGCTGTAGAGG‑CGCCCCAGTTCACGGCGGAGGAGTAACATTCCACGCT‑GGAGTGGTTGTTTTGTCTTTTGCCCGAAGccccgtccca                                                                    >  2:603253/1‑231 (MQ=255)                                                                                                                                                                            tCCACATTGCTGGCGGTTTTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAC‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATgg                                                                 <  1:879837/239‑1 (MQ=255)                                                                                                                                                                            tCCACATTGCTGGCGGTATTG‑TTTGCTGGGGGGTTTGCTAGCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATgg                                                                 <  1:90737/239‑1 (MQ=255)                                                                                                                                                                            tCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATgg                                                                 <  1:808188/239‑1 (MQ=255)                                                                                                                                                                             ccACATTGCTGGCGGTATTG‑TTTGCTGCGGGGGTTGCTATCGGTTGGTTGATTTTTGGGCTGGTCTGGCTGCGAGTTCGTG‑TTTTCCTGGGGCGCGCTGAACGTAAAAT‑AAAGCGACTGGGAAA‑ACAGCTTTTACCCGCGCCCGATGTGGATGTGGGGCGCCCAC‑CCTAACCGGCGGAGGAAAAATTATTTATGTG‑GGAGGTTTTGTTTTTGCATTTGCCTGTAAACGCGGACTATggg                                                                >  2:487855/1‑238 (MQ=255)                                                                                                                                                                                caTTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGGATAACTTTCTATGCT‑GGGGTTGTTGTTTTTGCTTTTGCCTGTAGCCCGTGCCTTTGGCCgg                                                             >  1:418178/1‑239 (MQ=255)                                                                                                                                                                                  ttGCTGGCGGTATTT‑TTTGCTGCGGGGGTTGGTTTAGGGTGGTTGATTTTTGGCCTTTTCTGGCTGCGAGTTCGTG‑TTTTCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTTTCCCCCGACTGCCGGTGCTGGAGGTC‑CGCACTCGCTAGGGGGGACGGAATAACCTTTTATGTTGGGGGTTGTTG‑TTCTCTTTTTGCCTGttggacgttccctgtgcccg                                                             >  2:355819/1‑218 (MQ=255)                                                                                                                                                                                  ttGCTGGCGGTATTG‑TTTGCTGCGGGGGTTTCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGGAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CCCACTCGTCCGCGGCGAAGGGATAACTTTCTTTGCT‑GGGGTGTTTGTTTTTGCTTTTTCCCTTAGCCCCCGCCTATgggggggg                                                           >  2:594534/1‑233 (MQ=255)                                                                                                                                                                                   ttcgggCGGTTTTT‑TTTTCTGGGGGTTTTGTTTTCGGTTTGTTGATTTGGGGCCTGTTCGGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGTGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGtat                                                          <  2:418178/235‑1 (MQ=255)                                                                                                                                                                                            tttttttttttctggggtTTTGTTATGGGTTGGTTGATTTGGGGCCTGTTTTGGCTGCGAGTTCGTT‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTGACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCg                                                  <  2:342941/225‑1 (MQ=255)                                                                                                                                                                                                       tgcGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGCCTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCCTTAGCCGCTGCCTTTGGCTGGGAAATGGGCCGCGGAAgggggg                                       >  2:1077781/1‑234 (MQ=255)                                                                                                                                                                                                       tgcGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTAAGGCTGGGATTTGGGCCGCaaaatggcg                                        >  2:145613/1‑233 (MQ=255)                                                                                                                                                                                                         cGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGGGTATGGGgcgccagagggggcca                                     >  2:655436/1‑227 (MQ=255)                                                                                                                                                                                                         cGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGGTTCTGCTTTTGCCTGTAGCCGCGGCCCATGGCGGGGTTAGGGGCCGCAGAAGGgcggca                                     >  2:597059/1‑236 (MQ=255)                                                                                                                                                                                                                tGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGGGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGGAATTTGGGCGCCGAAGGGCGCaaaaaaaaa                              >  2:793321/1‑237 (MQ=255)                                                                                                                                                                                                                        gttggttgATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTTGTATATGGGCGGCAGAAGTGCGCACCAAACcaagcaa                         >  1:247959/1‑236 (MQ=255)                                                                                                                                                                                                                        gttggttgATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGGGCGCAACAAAACAAGCAAGAt                      >  1:269178/1‑239 (MQ=255)                                                                                                                                                                                                                                 tttGTGGGCTGTTTTGGCTGGGGGTTCGTG‑TTTCCCGGGCGCCGGCTGAACGTAAAAT‑AAAGCGACGGGAAAC‑CCAGCTTTCCCCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCAAGATGAAGCCAAc             <  2:1018720/239‑1 (MQ=255)                                                                                                                                                                                                                                            ttCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTTTATGGGCCGCCGAAGTCGGCAACAAAACAAGCAAGGTGAAGCCAACCGCTTGTcgcg  >  1:564266/1‑239 (MQ=255)                                                                                                                                                                                                                                            ttCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCAAGATGAAGCCAACCGCTTGTcg    <  1:1008694/237‑1 (MQ=255)                                                                                                                                                                                                                                             tCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCAAGATGAAGCCAACCGCTTGTcgcg  <  1:594534/238‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  TAAGTTGACCATAATTTATTCGCTCTAACCACATAACGGGAAGTAATGTGAAATATTTACTCATTTTCTTACTGGTGTTAGCG‑ATCTTCGTGATTTCGGTCACGTTGGGTGCGCAGAACGATCAACAGGTGACGTTTAATTATCTGTTAGCGCAAGGGGAGTACCGTATTTCCACATTGCTGGCGGTATTG‑TTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTG‑TTTCCCTGGCGCGCGCTGAACGTAAAAT‑AAAGCGACTGGAAAA‑CCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGC‑CGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCT‑GGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCAAGATGAAGCCAACCGCTTGTCGCG  >  NC_000913/1340196‑1340667

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