Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,520,062	T→A	15.5%	*346K (TAA→AAA)	curA →	curcumin/dihydrocurcumin reductase, NADPH‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,520,062	0	T	A	15.5%	135.2 / 2.7	60	*346K (TAA→AAA)	curA	curcumin/dihydrocurcumin reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base T (30/19);  new base A (2/8);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 3.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTA‑TTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCA‑GGGTTGGACAGAAAACACGCGAAG‑GAATT‑AAAATGCCTGGAACGGGAAAAATGAAACATGTC  >  NC_000913/1519831‑1520291                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           aTGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTAACGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑GGGCGGGGGATGC‑TTAAAgat                                                                                                                                                                                                                                      >  2:200532/1‑236 (MQ=255)   ggCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAACCT‑TCGGCAAAGGGGTGAT‑CCCCG‑TGGCGGGTGATGA‑TTAAAGtttta                                                                                                                                                                                                                                   >  2:528689/1‑237 (MQ=255)      tACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGGGATGA‑TTAAAGTTTCAAcg                                                                                                                                                                                                                                >  2:150759/1‑239 (MQ=255)       tcAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTAATATCGGTCAGGATTAAGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAAcgg                                                                                                                                                                                                                               <  2:351665/238‑1 (MQ=255)        cAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTTAAcggc                                                                                                                                                                                                                              >  2:1148908/1‑239 (MQ=255)        cAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGGGAT‑CCGCG‑GGGCGGGTGTGGA‑TAAAAGTTTCAAcgcc                                                                                                                                                                                                                              >  2:526930/1‑237 (MQ=255)        cAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGGAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGGGGAGG‑TTAAAGTTTCAAcggc                                                                                                                                                                                                                              >  2:1019104/1‑239 (MQ=255)             cTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑GTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGCGGGGAT‑CCGCG‑GGGCGGGGGATGA‑TTAAAGGTTCAACGgcgtcgg                                                                                                                                                                                                                         >  2:555237/1‑238 (MQ=255)                     aCGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGCTTATGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGGGATGA‑TTAAAGTTTCAACGGCGGCGTAAgc‑cccc                                                                                                                                                                                                                >  2:167417/1‑239 (MQ=255)                     aCGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGGGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGCGGGGAT‑CCGCG‑TGGAGGGTGATGA‑TTAAAGATTAAAaggagagggtaagcgcg                                                                                                                                                                                                                 >  1:243765/1‑225 (MQ=255)                                                 gCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGCTAAAATCCACTACCGCG‑AAGAAATTACTGAAGGTTTAGAGAATGCGCAACAGACG‑ATTATCGGCCTGCTAACGGGGAAAAACT‑GCCGCAACATGGTGAT‑ACCCG‑TGGCGGGGGTGGA‑TAAAAATTTCAACAGCGGCGCAAGCACCACAAACAAAACCAGACaaaaaacaatttt                                                                                                                                                                                     >  1:622321/1‑231 (MQ=255)                                                 gCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGG‑TTAAAGTTTCAACGGCGGCGCAAAC‑CGCCATAAAAATCCAGGCCAAAAGCTTAttt                                                                                                                                                                                     >  2:638971/1‑238 (MQ=255)                                                        gtttgtgtcaccccacttgatggTTTAAGGGGGGGATGGGGGCAAGTGGAGAAGGAGGAAAAAACCCCCGCCGGAG‑AAGAAATCACTGGCGTTTGAGGCAAGGCGCCACAGCCT‑TGTGGCGGTCGGGTGAGGGATAAAATCT‑TCGAAAAGGGGTTGCT‑CCGCG‑TGGGTGGTGATGA‑TAAAACTTTCAAGGGCGGATAAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAg                                                                                                                                                                                 <  1:295871/214‑1 (MQ=255)                                                          tATGGTCACCGCATCCATGAGTTTCAGAGGGAGCT‑GGGGCAATGTGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGAGCCACAGACG‑TTTCTCGGCCTGCTGAAGGGGAAAACCT‑TCGGCACAGTGTGCAT‑CCGCG‑TGCCGGGAGGTGA‑GGAAAGGTTTACCGCCGGCTTAACC‑CCACTTAAAATTCCCCGGAAAAAAATTATTTTCaacgtttt                                                                                                                                                                           >  1:216091/1‑232 (MQ=255)                                                            tGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGGGATGA‑TTAAAGTTTCAACGGCGGCGTAAAC‑CGCCAACAAAAACCCGGCCAAAACGTTTTTTTTTGCTttaatt                                                                                                                                                                         >  2:52139/1‑239 (MQ=255)                                                                 aCCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAACGTGGTGAT‑CAGCG‑TGGCGGGAGATGT‑TTAAAGTTTCAACGGCTGCGTCAGT‑CTCAGATAAAATCCGGGCAAAAAACGTCTTTTTAGCCTTAATtaaata                                                                                                                                                                    >  2:759678/1‑239 (MQ=255)                                                                 aCCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAGA‑CGCCATAAAAATCCAGGGAAAAAGCTCATTTTTAACTTTAATtaaata                                                                                                                                                                    >  2:699943/1‑239 (MQ=255)                                                                 aCCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAAC‑CGCCCTAAAAAACCAGGCAAAAAACTTATTTTTAGCTTTAATtaacta                                                                                                                                                                    >  2:63549/1‑239 (MQ=255)                                                                 aCCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGGAAGC‑CGCCATAAAAAACCAGGCAAAAAACTTATTTTTAGCTTTAATttactt                                                                                                                                                                    >  2:275736/1‑238 (MQ=255)                                                                 aCCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACCGCGGGCGAAGC‑CGCCCTAAAAATTCAGGGAAAAAACTTATTTTCAACTTTAATtaacta                                                                                                                                                                    >  2:96243/1‑239 (MQ=255)                                                                 aCCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGGGATGA‑TTAAAGTTTCAACGGCGGCCGAAAC‑CGCCATAAAAATTCAGGGAAAAAAGTATTTTTCCGCTTTAATtaacta                                                                                                                                                                    >  2:1212944/1‑239 (MQ=255)                                                                 aCCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGATGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGGGGGTGAGGA‑TTAAAGTTTCAACGCCGGCGCAAGC‑GCCCGAAAAAAACCAAGCAAAAAACATATTTTAAGCAttatttaaaaa                                                                                                                                                                    >  1:868688/1‑236 (MQ=255)                                                                   cGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAGAGTGGTGAT‑CCGCG‑TGGCGGGGGATGA‑TTAAAGTTTCAACGGCGGGGTAAGC‑CGCCCTCAAAATCCCGGCAAAAACATTATTTTCAGATTTAATtaactaa                                                                                                                                                                   >  2:981489/1‑238 (MQ=255)                                                                                  cAGAGGGAGGT‑GGGGCAATGGGTGAAAGAGGAAAAAACCCCAGACCGCA‑AAGAAATTACTGACGGTTTAGAGAAGGCGCCACAGCCG‑TTTACCGGCCTGCAGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGGG‑TGGCGGGTGATGA‑TTAAAGTTCCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctg                                                                                                                                                   <  2:719066/239‑1 (MQ=255)                                                                                         agaT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGGAAGC‑CGCCCTAAAAATACCAGCAAAAAACTTTTTTTCACCTTTAATTTAATAAAcaacacaatgtttttatctcca                                                                                                                                            >  2:612144/1‑219 (MQ=255)                                                                                         agaT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCACCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGGGATGA‑TTAAAGTTTCAACGGCGGCGGAAGC‑CCCCAGAAAAAACCCAGGAAAAAAATTATTTTTGTCTTTAATTAAATAAacaccacatatttgtgctcacaa                                                                                                                                            >  1:480804/1‑220 (MQ=255)                                                                                           aT‑GGGGAAGTGTGTGAAAAAGGAAAAAGTCCCCT‑CCGCGAAAGAATTTGCGGAGGGTTTAAGGGATG‑GCCACGGGCGTTTTAGCGGCCGGCTGAAGGGTAAAAGCT‑TCGGCAAAGTGGAGAC‑CCGGG‑TGGCGGGTGATGA‑TTAAAGTTTAAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                                                                                                                                           <  1:150759/238‑1 (MQ=255)                                                                                              ggggggggggggGAAAGAGTATAAACACTTCCAGCGAG‑AAGAAATTGCGGGAGTTTTGGGGAGGGCGCCACAGTCG‑TTTAGCGGCCTG‑TGGGAGGAAAAAACTCTGGGCGAAGTGGTAAT‑CCGCG‑GGGCGGGTGATGA‑TAAATGTTTCAAGGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaataa                                                                                                                                         <  1:686288/230‑1 (MQ=255)                                                                                                          gAAAGAGGAAAAAACCCACTCCCGCG‑AAGAAATTACGGACGGTTGAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAC‑CCGCG‑GGGCGGGTGATGA‑TTAAAGTTACAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                            <  2:378165/239‑1 (MQ=255)                                                                                                          gAAAGAGAATAACATCCACCTCCAGT‑AAAAAATTCCTGACGATTCAGAGAATGCCCCACAGGTG‑TTTCTCGGCGTGGTGAAGGATAAAATCT‑TCGGAAAAGGTGTTAT‑CCGCG‑TGGCGGGTTATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                            <  1:1112725/239‑1 (MQ=255)                                                                                                          aaaaGAGGATAACACCCCTTCCGGAA‑AAAAAATTACTGGAGGTTGAGAGAATGCCCCACAGGCG‑TTTAGCGGCCGGGGAAAGGAAAAAACCC‑TGGACAAGGGGGGTAC‑CCGGG‑GGGCGGGGAATAA‑TTAAATTTTCATGGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                            <  1:1148908/238‑1 (MQ=255)                                                                                                           aaaGGGGATAAACTCCTCTCCCGAG‑AAGAATTTACTGGCGGTTTAGAGAAGGCGCCACAGACG‑TTTATCGGCCGGCTGAAGGGTAAAACTT‑GCGACAAAGTGGTGAT‑CCGGG‑TGGGGGGTGATGA‑TTAAAGTTTCAACGGGGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                            <  2:1211606/238‑1 (MQ=255)                                                                                                           aaaGGGGAAAAAACCCCTTGCCGAA‑AAAAAATTACGGACGGTTTAGAGATGGCGCCACAGGCG‑TTTAGCGGCCTGCTGAAGGTTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                            <  2:844205/238‑1 (MQ=255)                                                                                                           aaaGAGTAAAAAGTCCACTCCGGGA‑AAGAAATTACTGACGGTTGAGAGAAGGCGCCACAGACG‑TTTAGCGGCCTGTGGAAGGGTAAAACCT‑GCGGAAAAGTGGTGTT‑CCGCG‑TGGGGGGTGATGA‑TAAAAGTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                            <  2:334689/238‑1 (MQ=255)                                                                                                           aaaGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑tt                                                                                                                                                                                                                                            >  2:425861/1‑127 (MQ=255)                                                                                                           aaaGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑TGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                            <  2:130533/238‑1 (MQ=255)                                                                                                           aaaGAGGAAAAAATCCACTACGGAA‑AAAAAATTGCTGGCGGTTTAGAGAAGGCCACAGAGACT‑TTTATCGGCCTGCAGAAGGGTAAAAACT‑TCGGAAAAGTGGTGAC‑CCGCG‑TGGCGGGTAATGA‑TTAAAGTTTCACCGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                            <  2:1061309/238‑1 (MQ=255)                                                                                                           aaaGAGAATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑tt                                                                                                                                                                                                                                            <  1:425861/127‑1 (MQ=255)                                                                                                               aggaaaaactccctcacggag‑aagaaTTTACTGAGGTTTGAAAAAATGGCCCACAGAGG‑TTTATCGGCCTGCTGAAGGGTAAAACTT‑TCGGAAAAGTGGTGAT‑CCGCG‑TGGGGGGTGATGA‑TTAAAGTTCCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGtt                                                                                                                       <  2:67392/219‑1 (MQ=255)                                                                                                                 gATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCCCTCCACATCCAGGCAAAAAGCTTATTTTTAGCCTTAATTAAATAACCACATCATTGTTCTGCTGCTAATAACAACCAAGGAGAATGGGTTa                                                                                                                     >  2:189151/1‑239 (MQ=255)                                                                                                                 gATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCCCTAAAAATTCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACAACATTTTTCTGCTGCCAATAACAACAAAGGTGAGTGGTTTa                                                                                                                     >  2:110878/1‑239 (MQ=255)                                                                                                                       cttaccgacttag‑aaaaaATTCCTGATGGTTTAGAGAGTCCGCAACAGTCT‑TTTGTCGTCCTGATGAAGGGTAAAAACT‑CCGGCATGGTGG‑GAA‑CCGCGATGGGGGGCGATAA‑AAAAATTTTCAACGGCCGCGTAAGC‑CGCCATAAAAATCCATGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAATGTGAGTGGTTTATGatta                                                                                                               <  1:1011704/227‑1 (MQ=255)                                                                                                                         ccACTACCGAG‑AAAAAATTGCTGAGGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCAGAAGGATAAAACCT‑GCGGCAAAGTGGTGAG‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGatta‑tt                                                                                                            <  2:243765/239‑1 (MQ=255)                                                                                                                           acacccGCA‑AAGAAATGACTGACGTTTGAGAAAATGCGCCCGAGTCG‑TTTACTGGCGGTCTGAAGGTCAAAACCT‑CCGGCAAAGTGGGAAC‑CCGCG‑TGGCGTGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTA‑TTCa                                                                                                          <  2:589496/235‑1 (MQ=255)                                                                                                                            catcccgaagagaaATTCCTGGCGGTTTAGAAAAGGCGCCACGGTCG‑TTTATCGCCCTTCTGAAGGGTAAAACCG‑TCGGAAATGGGGTGAT‑CCGGG‑TGGCGGGTGATGA‑TTAAATTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTA‑TTCa                                                                                                          <  2:344873/229‑1 (MQ=255)                                                                                                                                cgcg‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAAATCATTGTTTTGCTGCTAATAAAAACAAAGGTGAGTGGTTTTTGATTA‑TTCACTTTa                                                                                                     >  2:1170703/1‑239 (MQ=255)                                                                                                                                  ggaaaaaaaTTACGGACGGTTGAGGGAATGCGCCAGAGCCG‑TTTACCGGCCGGAGGAAGGTAAAAAACG‑TCGGAAAGGTGGTGAC‑CGGGG‑TGGCGGGTGATGA‑TTAAATTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTA‑TTCACTTTaa                                                                                                    <  2:480804/233‑1 (MQ=255)                                                                                                                                       aaatgttgggggtgTTTAGAAAATGCGCCACAGGGT‑TTTGACGGGTTG‑AGAAGGGTAAAAACT‑TTGGGAAGGGGGTGTTCCCGGG‑TGGGGGTTTATGT‑TAAAAGTTGAATGTGCGTCGTAACC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCTCATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTA‑TTCACTTTAATacac                                                                                               <  1:275736/226‑1 (MQ=255)                                                                                                                                          aTTACTGACGGTTTAGAGAATGCGCCACAGATG‑TTTATCGGCCTGCGGAGGGGTAAAAACT‑TCGGAAAAGTGGTGAT‑CCGCG‑TGGGGGGTAATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTA‑TTCACTTTAATACACCAg                                                                                            <  2:868688/239‑1 (MQ=255)                                                                                                                                                              tGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAAAAACAAAGGTGAGAGGGTTATGATTA‑TTCACTTTAATACACCAGGTGAATTCCTTCTTCCAttc                                                                        >  2:74614/1‑237 (MQ=255)                                                                                                                                                               gcgcCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCCCTAAAAATCCCGGCAAAAAGCCTATTTTCAGCTTTTATTAACTAACCACATCATTGTTCTGCTGCTAATAAAAAAAAAGGGGAATGGTTTTTGATTT‑TTTACTTTTATACACCAGGGGAATTACtttttgcctgcc                                                                       >  2:1237354/1‑229 (MQ=255)                                                                                                                                                                  ccacaatggtttgtcggggcgggtgATGGATAAAATGT‑TGGGCAAGGTGGTGAC‑CCG‑G‑TGGGGGGGTAAGAGTAAAAGTTGCACCGGGGGTGACAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTA‑TTCACTTTAATACACCAGGTGAATTCCTTCTGCCATgcagg                                                                     <  1:83724/216‑1 (MQ=255)                                                                                                                                                                   atcaGACT‑TGTACCGGCCTGCGGAAGGGTAAAAACT‑TGGGCATAGTGGTGCT‑CCGCG‑TGGCTGGTGATGA‑TAAAAGGTTCAGCGGCGGCGTAAGC‑CGCCATAAAAAGCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTA‑TTCACTTTAATACACCAGGTGAATTCCTTCTGCCATgcaggca                                                                   <  2:677788/237‑1 (MQ=255)                                                                                                                                                                       agggttttgccGGCCTGCTGAAGGGTAAAAGCT‑TCGGAAAAGGGGTGAC‑CCGGG‑TGGCGGGGTATGA‑TTAAAGTTCAAACGGGGGCGAAAGC‑CGGCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTGAATTAACTAACCACATCATTGTTCGGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTA‑TTCACTGTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCA‑ggg                                                               <  2:929823/229‑1 (MQ=255)                                                                                                                                                                                     ctgctgAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTTAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTTTGATTA‑TTAACTTTAATACACCAGGGGAATTCCTTCTGCCATGCAGGCA‑ggggtctgtctctttta                                                 >  2:461673/1‑225 (MQ=255)                                                                                                                                                                                           gaagggaaatatctgtggaagggtggtgAT‑CCGGG‑GGGAGGGGTAGGA‑AAAAAGGTTAACGGGCGACTAAAGC‑CGACATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTA‑TTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCA‑GGGTTGGACAGAAAACACGCGa                                            <  1:96243/217‑1 (MQ=255)                                                                                                                                                                                               gTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGGGGGTGATGA‑TTAAAGTTACAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGatta‑tt                                                                                                            <  2:1186036/171‑1 (MQ=255)                                                                                                                                                                                                          ggtcaaagggggact‑CCGCG‑TGGGGGGTTATTA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCAATAAAAATCCAGGCAAAAAGCTTAGTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTA‑TTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCA‑GGGTTGGACAGAAAACACGCGAAG‑GAATT‑AAAATGcc                           <  2:1201520/225‑1 (MQ=255)                                                                                                                                                                                                                             cg‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTGAGC‑CGCCATAAAAAACCAGGGAAAAAGCTTATTTTTAGCCTTAAATTACTAACCAACTCCTTTTTTTGCTGCTAATAAAAACAAAAGGGAGTGGTTTATTATTATTTCATTTAAATACA‑CAGGTGAATTTCTTTTGCCATTGAGGAC‑GGGTTTGACCAAAAACAACCGGAGAGAATT‑AAAAAG‑CTGGAAAGGGaaaaaaaaa         >  1:1066777/1‑235 (MQ=255)                                                                                                                                                                                                                             cg‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAAC‑CGCCCTAAAAATCCAGGCCAAAAGCTTATTTTCCGCTTTTATTAACTAAACACAACATTTTTCTGCTGCTAATAACAACAAAGGTGAGGGGGTTATAATTT‑TTTACTTTAATACACCCAGGGAAATCCTTTTGCCCTGTAGGCAGGGGGTGGA‑AGAAAAAACACGAAG‑GAATTAAAAAACCCTGGAAAGGGAAAAAtaa          >  1:830600/1‑237 (MQ=255)                                                                                                                                                                                                                                   cGGGTGATGA‑TTAAAGTTTCAACGGCGGCGGAAGC‑CGCCATAAAAATCCAGGCCAAAAGCTTATTTTCAACTTTAATTAACTAACCAAATCATTGTTCTTCTGCTAATAACCACAAAAGTGAGTGGTTTATGATTT‑TTTACTTTAATTAACCCGGTGAATTCCTTTTGCCATTTCGGAC‑CGGTTTGACATAAAACACCCGCAG‑AAATT‑AAAAAGCCTGGGACGGGgaaaaaaaaaaaaa    >  1:961064/1‑230 (MQ=255)                                                                                                                                                                                                                                     ggggatga‑tTAAAGTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATTCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAAATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTA‑TTCACTTTAATAACCCAGGTGAAATCCTTCTGCCATGCCGGCA‑GGGGTGGACCGAAAACACGGGAAG‑GAATT‑AAAATGCCCGGAACGGGGAAAATGAAACATGta  >  2:326391/4‑238 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           ATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGAT‑GGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCG‑AAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACG‑TTTATCGGCCTGCTGAAGGGTAAAAACT‑TCGGCAAAGTGGTGAT‑CCGCG‑TGGCGGGTGATGA‑TTAAAGTTTCAACGGCGGCGTAAGC‑CGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTA‑TTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCA‑GGGTTGGACAGAAAACACGCGAAG‑GAATT‑AAAATGCCTGGAACGGGAAAAATGAAACATGTC  >  NC_000913/1519831‑1520291

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