Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,151,328	T→G	28.1%	A590A (GCA→GCC)	ppc ←	phosphoenolpyruvate carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,151,328	0	T	G	28.1%	87.1 / 2.7	71	A590A (GCA→GCC)	ppc	phosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base T (30/12);  new base G (8/14);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 8.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTT  >  NC_000913/4151093‑4151564                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                ggTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTgcg                                                                                                                                                                                                                                              <  1:1188548/239‑1 (MQ=255)          cGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGCGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGG‑TGTGACCACAGCCCC‑GCATGAGCCGGGGCG‑CCGC‑cgcg                                                                                                                                                                                                                                    >  1:93919/1‑238 (MQ=255)             aGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCTTTTCAGGCTTCCTGGCGGG‑TGTGACAGCAGCGCC‑GCATGAGCAGGGGCG‑CCGC‑CGCGAccc                                                                                                                                                                                                                                >  2:1079390/1‑238 (MQ=255)             aGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CCCGACCa                                                                                                                                                                                                                                >  1:103604/1‑239 (MQ=255)               gagaTGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TTTGACAGCAGCGCC‑GCATGAGCAGGGGCG‑CCGC‑CGCGACCA‑At                                                                                                                                                                                                                             >  1:1157237/1‑239 (MQ=255)                    gactgacAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGGGCG‑CCGC‑CGCGACCA‑ATGGaa                                                                                                                                                                                                                         >  1:385862/1‑238 (MQ=255)                         acAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGGTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑TCCGTTCAGGCTTCCTGGCCGG‑CGTGACCACCACGCC‑CCATGAGCAGGTGcg‑ccgc‑ccccgcccctggtaccccc                                                                                                                                                                                                                     >  2:271147/1‑220 (MQ=255)                          taGTTCCGCCATAGGGCGACGCCGGCTTTTTTTCGGCGGCGGGGGGGGCGGCAGGTTGGCTCCCAGACCCGCCCCGGTAAAAGGCGACAGGCTGCGGGGGGTTTTTTCT‑GGCAGACCATATTTAAAGCGGATCATCCCGCCCTTTT‑CGTTTA‑CGCGCGGCCC‑GCTTTTCAGGCTTCCGGGCGTT‑TGTGACAGCAGCGCC‑GCAGGAGCGGGGGCG‑CCGC‑CGCGCCAA‑AGGGAACcgccgc                                                                                                                                                                                                                  <  2:914727/238‑1 (MQ=255)                                        tGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGCCCGTGGaa                                                                                                                                                                                                        >  1:583073/1‑235 (MQ=255)                                            aCGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGACTCGCCCCCGGATAAAGCGACAGGCTGCTGACAGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCCCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGA‑TGTGGCAACCGCGCC‑CCATGAGCAGGGGGG‑CCGC‑CGCCAGCA‑ACGGAAACGCCCCCACCgctgggcacacaca                                                                                                                                                                                                >  2:118712/1‑226 (MQ=255)                                            aCGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACCGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TTGTACAGCAGCCCC‑GCCTCAGCAGGCGCC‑GCGC‑CCGCACCA‑ATGGGACCGCCCCGCCAGCTGGACAACGtgc                                                                                                                                                                                                >  2:870513/1‑237 (MQ=255)                                            aCGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCAAGAGCCGGTGCG‑CCGC‑CGCGGCCC‑ATGGGACCGCCGCGCCCGGGGGACAACg                                                                                                                                                                                                   >  1:521087/1‑236 (MQ=255)                                               cgacttcgctgccgccTCGGGCGGTGGCAGGTGTTTTGCTCCCGGACTCCCCTCGTTATAAAGCCACAGCCTGCGGAC‑GTTTTTTCTGGGCAGACATTATTAAAAGCGGATCTTCTCGCCTGTTT‑CGTCTC‑CCCGCAGCCC‑TCTTTACAGGTTTCTGGGCGGT‑TGTCCCACCGGCCCC‑GCAGCACCGGGCGCC‑CCCC‑CGCGACAA‑ATGGATCCGCCGCGACCGTGGAACAACGTCAGCt                                                                                                                                                                                             <  2:803485/224‑1 (MQ=255)                                                   ctctctTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCCACCA‑ATGGGACCGCCCCGCCCGTGGAACAACGTTAGCTcgac                                                                                                                                                                                         >  1:564336/1‑236 (MQ=255)                                                   ctctctTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGGGGCAGCAACCCC‑GCATGAGCAGGGGCG‑CCGC‑CGCGACCA‑AAGGCACCGCCGCGACCGTGGAGCAACGGCAGCTCAAt                                                                                                                                                                                         >  2:744715/1‑239 (MQ=255)                                                                  ggcggTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGGGGAACAACGTCACGTCAATACCCGCCTTTTCGCa                                                                                                                                                                          >  1:722216/1‑239 (MQ=255)                                                                      ggggcggcGGGTGGGCTTCCAGACCCGCCCGGGTAGAAAGCGACGGGTTGCTGGCGTTGTTTTCT‑GGCAGCCTATTTTTAAGGGGTATCCTCCCGCCTTGTT‑GGGTTC‑CGCGCGGGCC‑GCCTTTCGGGCTTCCTGGGGGT‑GGGGCCGGCGGGCCC‑GCAGGAGCGGGCGCG‑CCGC‑CGCGACAA‑TTGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGt                                                                                                                                                                       <  2:103604/236‑1 (MQ=255)                                                                             agggTTGGCTCACAGAATCGCCCGGTTAAAAAGCGACGGGCTGCTGACGGTGATTTCG‑GGCAGACAATATTTAAGGGGGATCTTCTCGCCCTGTT‑GGGTTC‑CGCGCGGGCC‑CCTTTCCAGCTTTCGTGGGGGT‑TGTGACAGCGGCGCG‑GCATGGGCGGGGGCG‑CCGC‑CGCGCCAA‑GTGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGAtt                                                                                                                                                                <  2:1237554/236‑1 (MQ=255)                                                                                 gtGGCCTCCAGACTCCCCCGGGTATAAGGCGACAGGTTGCGGGGGGTGTTTCCT‑GCCAGACTATATTTAAAGCGACTCCTCCCCCTTTGTG‑CGGTCC‑CGCGCGGGCC‑GCTTTTCGGGCTTCCGGGGGGT‑TGTGCCAGCAGCGCC‑GCAGGAGCGGGGGCG‑CCGC‑CGCGACCA‑AGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATg                                                                                                                                                            <  2:1157237/237‑1 (MQ=255)                                                                                  gtgctcccaactccGCCCGGGTATAAAGCGCCAGGCGGCGGAGGGTGTTTCCC‑GGCAAACCATATTAAAAGCGGACCACCCGCCCCTGTT‑CGGTAC‑CCCGCAGGCC‑GCCTTTCAGTCTCCGTGGGGGT‑TGCGACGCCGGCCCC‑CCATGAGCGGGTGCG‑CCGC‑CGCCACCA‑TTGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCa                                                                                                                                                          <  2:1169161/226‑1 (MQ=255)                                                                                          aGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCtgtg                                                                                                                                                   >  1:366494/1‑238 (MQ=255)                                                                                              tCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGa                                                                                                                                              >  2:950754/1‑239 (MQ=255)                                                                                                                aCAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATg                                                                                                                                                            >  1:720331/1‑207 (MQ=255)                                                                                                                aCAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGaa                                                                                                                               >  1:25981/1‑236 (MQ=255)                                                                                                                aCAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAggt                                                                                                                            >  1:754438/1‑237 (MQ=255)                                                                                                                 cAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCC‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCt                                                                                                                            >  1:147317/1‑238 (MQ=255)                                                                                                                 cAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGAG‑CAGC‑CGCGACCA‑ATGGAACCGCCGCGACCGGGGAACAACGTAAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATACATCCTGTGCCTGAGTTTGGCCCCAGGAAGCt                                                                                                                            >  2:1085197/1‑238 (MQ=255)                                                                                                                 cAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGC‑TGTGACAGCCGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CCCGACCA‑ATGGAACCGCCGCGACCGTGGAAAAACGTCAGCTCACTCCCCGCTTTTTCGCAGGGTTTGGTTAATGCATCCTGTGCCTGACATTCCGCCCAGGAAGCt                                                                                                                            >  2:7899/1‑238 (MQ=255)                                                                                                                 cAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTGAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACCGCAGCGCC‑GCATGAGCAGGGGAG‑CCGC‑AGAGACCA‑ATGGTAACGCCGCGACCGTGGGACAACGTCAGCACAATACCCCATTTTTCTCTGGATTTTATTTATGCGTCCTGGGCCTGATATAGTGCCCAGAAAGCt                                                                                                                            >  2:248115/1‑238 (MQ=255)                                                                                                                                             aCCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCCGCCCAATACCCGCCTTTTCGCCGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGGAGCTGCCATCACTCCCTCATCTTTTGCttgg                                                                                                 >  2:1150789/1‑235 (MQ=255)                                                                                                                                                 ttttaaagccggtccctcccGCCCTTTT‑CGTTCC‑CGCGCGGGCC‑CCTTTTCGGGCTTCCGGGCGGT‑GGGGACAGCGGCGCC‑CCATGAGCGGGTGCG‑CCGC‑CGCGCCCA‑ATGAATCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGa                                                                                            <  2:454460/220‑1 (MQ=255)                                                                                                                                                        aGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAACGCC‑GCCTGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGGGGAACAACGTCAGCTCAATCCCCGCTTTTTCGCAGGCTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCACGAAGCTGCCATCACTCCCGCATCTTTTTCTGAGTCGGAATAGcgc                                                                                     >  2:1159379/1‑237 (MQ=255)                                                                                                                                                            gATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGCTATTGCGCCCAGGAAGCTGCCATCACTCCCGCCTCTTTTGCTGAGTCGGAATAGCCAat                                                                                   >  1:271007/1‑237 (MQ=255)                                                                                                                                                            gATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAatca                                                                                 >  1:316002/1‑239 (MQ=255)                                                                                                                                                             atcatcTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGc                                                                                                                                                                                              >  2:860045/1‑129 (MQ=255)                                                                                                                                                             atcatcTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGCGGCACAACGTCAGGTCAATACCCGCTTTTTCGCAGGGTTTTATTAATGTATCCTGGGCCTGATATTGCGCCCAGGAACCTGCTAACACTCCCGCATCTTTTCCTGAGGCGGAATAGCCAataa                                                                                 >  2:997699/1‑236 (MQ=255)                                                                                                                                                              cccccccccccTG‑T‑CGGTTACCGCGCGGCCC‑CCCTTCCGGCCTCCCGGGCGGT‑GGTGACAGCAGCGCC‑GCGTGAGCGGGCGCC‑CCCC‑CGCGCCAA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAatcat                                                                                <  2:25981/230‑1 (MQ=255)                                                                                                                                                                  tccGCCCTG‑T‑CGGTGA‑CGGGCCGGCCGGCCGTCCGGGCTTCTGGGGGGT‑TGTGAC‑GCAGGGCCAGAGGGGGAAGGGGCG‑CGGC‑CGCGCCCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCa                                                                           <  2:719938/237‑1 (MQ=255)                                                                                                                                                                   cgcgccTGTTCCGGTGC‑CGCGCGGGCG‑GCCTTGCAGCCTTCCGGGGGGT‑TGTGACAGCAGCGCC‑GCAGGGGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAGCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCa                                                                           <  1:121550/235‑1 (MQ=255)                                                                                                                                                                       ccTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGGC‑GCATGAGCGGGGGCG‑CCGC‑CGCGGCCA‑AAGGAACCGGCGTGACCGTGGAACAACGTGAGCTCAATAACCGCTTTTTCGCAGGTTTTGGTTAAGGCCTGCTGTGCCTGATATTGCGCCCAGGAAGCTGCCCTCACTCCCGCATCTTTTGCTGATTCGGAATAGCCCATCATCAACATCTGt                                                                      >  2:89008/1‑239 (MQ=255)                                                                                                                                                                           tt‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCCTCTTTTGCTGGGTCGGAATAGCCAATCATCACCATCTGTTTg                                                                   >  1:484350/1‑238 (MQ=255)                                                                                                                                                                              gcgttaccgcgcgGGCC‑GCTTTACGGGCTCCCGGGCGGT‑TGTGACAGCAGCGCC‑ACATGAGCGGGTGCG‑CCGC‑CGCGACAA‑AGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGcc                                                                 <  1:930179/232‑1 (MQ=255)                                                                                                                                                                                 ttA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGa                                                             >  1:843362/1‑239 (MQ=255)                                                                                                                                                                                      ggccgggCC‑GCCTTCCAGGCTCCGGGGCGGT‑GTTGACACCCGCCGC‑GCGTGAGCGGGGGCC‑CCCC‑CGCACCCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCa                                                         <  2:883544/234‑1 (MQ=255)                                                                                                                                                                                       cgcAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCC‑ATGGAACCGCCGCGACCGTGGAACAACGTCCGCTCAATACCCGCTTTTTCGCCGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCCCTCCCGCCTCTTTTGTTGAGTCGGAATAACCAATCATCCCCATCTGTTTGCCCTGAATccg                                                        >  2:706456/1‑237 (MQ=255)                                                                                                                                                                                       cgcAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTTCTGAGTCGGAATAGCCAATCATCCCCCTCTGTTTGCCCTGAATCAg                                                        >  2:282279/1‑239 (MQ=255)                                                                                                                                                                                          aGGCC‑GCCTTTCGGGCTTCCTGGCGGT‑TGTGACGGCAGCGCC‑GCAGGAGCGGGGGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGc                                                      <  1:1068068/238‑1 (MQ=255)                                                                                                                                                                                            gcc‑gccTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCAc                                                   >  2:741318/1‑239 (MQ=255)                                                                                                                                                                                            gcc‑gccTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCAc                                                   >  1:883544/1‑239 (MQ=255)                                                                                                                                                                                             gtccccTTTCCTGCGTCCCGGTGGTGTGTGACGGGAGCGCC‑GCATGGGCGGCTCCG‑CCCC‑CGCGCACA‑TAGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCa                                                    <  2:564336/235‑1 (MQ=255)                                                                                                                                                                                                 tcttcagggCTTCCTGGGGGT‑GGGCACAGCGGCGCC‑GGAGGAGCGGGGGCG‑CCGC‑CGCGACAA‑AGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa                                               <  2:271007/232‑1 (MQ=255)                                                                                                                                                                                                 cgttgcggtctgcggggCGGT‑TGTGACGGCGGCGCC‑GAAGGGGAAGGTCGG‑CCGC‑GGCGACCA‑AGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa                                               <  1:280343/224‑1 (MQ=255)                                                                                                                                                                                                 ccTTTCAGGCTTTCGGGCGGT‑TGTGACGACAGCGGC‑GTGGGGGCAGGTGCG‑CCCC‑CGCAACCA‑GTGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa                                               <  1:283325/239‑1 (MQ=255)                                                                                                                                                                                                 ccTTTCAGGCTTCCTGGCGGT‑TGTGACAGCCGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTTCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGAt                                                >  1:969262/1‑238 (MQ=255)                                                                                                                                                                                                  cgtccAGGCTTCCGGGGGGT‑TGCAGCAGCCGCCCC‑GCAGGAGCAGGGGCC‑CCCC‑CGCGCCAA‑ATGGACCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa                                               <  2:366494/234‑1 (MQ=255)                                                                                                                                                                                                      cggcctgccgggggtt‑tGGGACAGCAGCCCC‑GCATGGGCAGGTGCC‑CCGC‑CGCAACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGt                                          <  1:744715/224‑1 (MQ=255)                                                                                                                                                                                                      cAGCCTTCCTGGTGGT‑TGTGCCAGGAGCGCC‑GCAGGAGCAGGGCCC‑CCCC‑CGCGCCCA‑AGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGt                                          <  2:316002/239‑1 (MQ=255)                                                                                                                                                                                                           cgccctgccggtttgagcaggagCGCC‑GAAGGCGCAGGCGCG‑CCGC‑CGCGACAG‑ATGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAAGAGCCAATCATCACCAGCTGTTTGCCCTGAATCAGGCCACGATACCAGt                                          <  1:1057511/218‑1 (MQ=255)                                                                                                                                                                                                                      tgtggACAGCAGCCCC‑GCATGAGCAGGGCCC‑CCGC‑CGCAACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGagcagc                            <  1:747851/234‑1 (MQ=255)                                                                                                                                                                                                                       gtggcccacccgccccggaggAGCAGGTGCC‑GCCC‑CGCGACAA‑AGGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGTTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGagcagc                            <  2:654857/218‑1 (MQ=255)                                                                                                                                                                                                                                ggcgcC‑ACAGGAGCGGGGGCG‑CCGC‑GGCGACAA‑AGGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGAc                 <  2:266304/238‑1 (MQ=255)                                                                                                                                                                                                                                cgcgcG‑GCAGGAGCAGG‑GCG‑CCCCACGGGGCCATATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATCCCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCTCCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGa                  <  1:749273/238‑1 (MQ=255)                                                                                                                                                                                                                                acggcc‑gcGTGAGC‑GGGGCGCCCGC‑CGCGACCA‑ATGCAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGAc                 <  1:672314/236‑1 (MQ=255)                                                                                                                                                                                                                                  cgcc‑gcATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACa                >  2:454951/1‑238 (MQ=255)                                                                                                                                                                                                                                    cc‑gcATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCATTATTAAGCAGCTGGGTCATGACATCg             >  2:1138433/1‑239 (MQ=255)                                                                                                                                                                                                                                     c‑gcATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGt            >  1:56143/1‑239 (MQ=255)                                                                                                                                                                                                                                      ggcATGGGCAGGAGCG‑CCGC‑CGCGACCA‑AGGGAGCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCGTGACATCGt            <  1:1159379/238‑1 (MQ=255)                                                                                                                                                                                                                                      ggcAGGAGCAGGGGCG‑CCGC‑CGCGCCAA‑GTGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGt            <  1:282279/238‑1 (MQ=255)                                                                                                                                                                                                                                            ggCAGGGGCG‑CCGC‑CCCAACCA‑ATGGACCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCg       <  1:947579/237‑1 (MQ=255)                                                                                                                                                                                                                                             ggaGGTCCC‑CCCC‑CGCGCCAA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCg       <  1:141469/235‑1 (MQ=255)                                                                                                                                                                                                                                                 gTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCgttgtt  <  1:741318/238‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                GGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCT‑GGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTT  >  NC_000913/4151093‑4151564

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