Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,151,323	A→G	14.7%	A592A (GCT→GCC)	ppc ←	phosphoenolpyruvate carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,151,323	0	A	G	14.7%	142.7 / 2.7	63	A592A (GCT→GCC)	ppc	phosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base A (41/11);  new base G (5/5);  total (47/16)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

GCCGCGGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGG  >  NC_000913/4151089‑4151558                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                gCCGCGGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑g                                                                                                                                                                                                                                           >  1:786226/1‑239 (MQ=255)          gACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTgcggcg                                                                                                                                                                                                                                   >  1:825669/1‑238 (MQ=255)               cGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGGCAGCAGCGGCGCATGGG‑CA‑GGTGCGCCGCcgcg                                                                                                                                                                                                                              >  1:937768/1‑238 (MQ=255)                         gactgacAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTTTGACAGCAGCGCCGCATGGG‑CA‑GGTGCGCCGCCGCCACCAATGGa‑a                                                                                                                                                                                                                   >  1:40720/1‑238 (MQ=255)                         gactgacAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGCTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGa‑a                                                                                                                                                                                                                   >  1:810727/1‑238 (MQ=255)                            tgacAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGGGACAGCAGCGCCGCATGAG‑CA‑GGGGCGCCGCCGCGACCAATGGA‑Accg                                                                                                                                                                                                                >  1:496399/1‑238 (MQ=255)                              acAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGACAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTCGCGGGTGTGGCAGCAGCCCCGCCTCGG‑CA‑AGATCGCCGCCCCGGCCAATGGacacacc                                                                                                                                                                                                               >  2:111395/1‑232 (MQ=255)                                                 aCGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTTTGACAGCAGCGCCGCATGAG‑CA‑GGGGCGCCGCCCCGACCAATGGG‑AACGCCGCGGCCGGGGACCAACGTCa                                                                                                                                                                                          >  1:268179/1‑239 (MQ=255)                                                 aCGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGGGGGACAACGTc                                                                                                                                                                                           >  1:124247/1‑238 (MQ=255)                                                  cGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGGG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGGACAACGTcc                                                                                                                                                                                          >  1:894230/1‑237 (MQ=255)                                                  cGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGCG‑CA‑GGTGCGCCGCCGCGACCAATGGG‑ACCGCCCCGACCGCGGAACAACGTCa                                                                                                                                                                                          >  1:378021/1‑238 (MQ=255)                                                              ttCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATAacc                                                                                                                                                                               >  1:70284/1‑237 (MQ=255)                                                                  gCTCCGGCGGTGGCAGCAGGTTGGATTCCAGCATCGCCCCGGTATCAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCCGTTCGGTTACTCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAACAGCCCATCAG‑CA‑CGCGCGCCGGCCGGAACAATGGG‑ACCGCCGCGCCCCTGGAACAACGGACGCTCGACACCCGCtgtt                                                                                                                                                                         >  2:575590/1‑239 (MQ=255)                                                                      gggcggTGGCGGGAGTTTGGCTGCAAAACCCCCCCGGGTAAAAGGCGACGGGCTGGTGGCGGTTTTTTGTGGCATA‑CATTTTTAAAGGCGGATCATCCCCCCCTTTGCGTTTCCGCGGGGGCCCCTTTTCGGGTTTCGGGGGGGTTGGAAAGGCGGCGCC‑CATGGG‑GAGGGGGCCCCCCCGCCACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCg                                                                                                                                                                      <  2:236658/237‑1 (MQ=255)                                                                      cggcggTGGCAGGGGGTTGGCTTCCAGAATCCCCCGGGTAAAACGCGACAGTCTGCGGCGGGTGTTTTCGGGCAGACCATTTTAAAG‑GCGGTTCTTCTCGCCTTGTTCGGTCACGCGCGGGCCGCCTTTCGGGCTTCCGGGGGGTTGTGACCGCAGCCCCGCAGGAG‑CG‑GGCGCCCCCCCGCGACCAATGAA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCg                                                                                                                                                                      <  2:769381/238‑1 (MQ=255)                                                                       ggcggTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGGCAGCACCGCCGCATGAG‑CA‑GGGGCGCCGCCGCGGCCAATGGA‑ACCGCCGCCGCCGTGGGACAACGTCAGCCCAAAACCCCCTTTTTCGCa                                                                                                                                                                    >  1:621207/1‑239 (MQ=255)                                                                       ggcggTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGGTTTGACCGCAGCGCCGCATGGGACA‑GG‑GCGCCGCCCCCACCCATGGA‑ACCCCCCCGGCCCTGGAACCACGTCAGCCCAAAATCCCCttttttccc                                                                                                                                                                    >  2:249472/1‑235 (MQ=255)                                                                                tggagtttGGCTGCCAAAGTCGCCCTGTTAAAAGGCGAGGGGCGGCGGCCGGTGTTTTCTGCCAGCCCATTTTTAAA‑CGGGTTCACCCCCCCTGGTTCGGTAGCGCCCGGCCCTCTTTCCGGGCTCCCTGGCGTTTGTGAGGGCGGCGCCGCATGAG‑AG‑GGTGCCCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGa                                                                                                                                                            <  2:172746/232‑1 (MQ=255)                                                                                  cagGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCCCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCCATACCCGCCTTTTCGCCGGTTTTGATTa                                                                                                                                                         >  1:663099/1‑239 (MQ=255)                                                                                        ggCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCa                                                                                                                                                    >  1:299960/1‑238 (MQ=255)                                                                                                      ccccGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGatat                                                                                                                                     >  1:543359/1‑239 (MQ=255)                                                                                                             aaaaaGCGAGAGGCTGCGGAGGGTGATCTCTGTCAGACATTTTTTAAA‑CCGGTCCATCCCGCTGTTTCCGTTCACCCGCAGGCCCCCTTTCAGGCTTCCTGGCGGTTGGGACCGCAGCGCCGCAGGAG‑AG‑GGGGCGCCCCCGCAACAAGTGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGccc                                                                                                                             <  2:360996/237‑1 (MQ=255)                                                                                                                     aCAGGCTGCTGACGGTGATTTCTTGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGCTGTGAAAGCAGCGCCGCATGAG‑CC‑AGTGCGCCGCCTCGCCCAATGGA‑AACGCCGCGACCGTGGAACAACGGCAGCTCAATACCCGCCCTTTTGCCGGGTTTGCTTAAATCCTCCTTCGCCTGATCTTCGCCCCCAgaa                                                                                                                         >  2:594373/1‑234 (MQ=255)                                                                                                                     aCAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGaa                                                                                                                        >  1:108804/1‑236 (MQ=255)                                                                                                                      cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGTCAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCt                                                                                                                     >  1:738891/1‑238 (MQ=255)                                                                                                                      cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCt                                                                                                                     >  1:820179/1‑238 (MQ=255)                                                                                                                      cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCCCATGAG‑CA‑GGGGCGCCGCCCCGACCAATGGA‑ACCGCCGCGGCCGGGGAACAAAGTCATCTCAATACCCGCTTTTTCTCAAGTTTTTATTAATGCATCCTTGTCCTGAGATT‑GCGCCCAGGAAgat                                                                                                                     >  2:959667/1‑236 (MQ=255)                                                                                                                      cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAACAGCGCCGCATGAG‑CA‑GGTGCCCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCACAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCt                                                                                                                     >  1:286352/1‑238 (MQ=255)                                                                                                                         gctgctGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCa                                                                                                                 >  1:448614/1‑239 (MQ=255)                                                                                                                                aCGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTc                                                                                                           >  1:828524/1‑238 (MQ=255)                                                                                                                                   gTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGAAAGCCGCGCCGCATGAG‑CG‑GGGGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGCGGAAAAACGTCAGCTCAATACCCGCTTTTTCCCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCg                                                                                                        >  2:733390/1‑238 (MQ=255)                                                                                                                                                   cATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATg                                                                                                                                                      >  1:1020930/1‑177 (MQ=255)                                                                                                                                                     tatTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGa                                                                                     >  1:837849/1‑239 (MQ=255)                                                                                                                                                          aaa‑GCGGCTCCTCTCCCCGTGTTGGTTGCCGCGCAGCCCCCCTTTCGGCTTTCGTGGCGTTTGGGGCGGGAGCGCCGAAGGGG‑CA‑GGTGCCCCCCCCCGCCCAAGGGA‑TCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGc                                                                                <  2:738891/239‑1 (MQ=255)                                                                                                                                                           aa‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGcc                                                                               >  1:284041/1‑239 (MQ=255)                                                                                                                                                            a‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGGGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTTCCATGACTCCCGCACCTTTTGCTGCGTCGGAATAGCCa                                                                              >  1:919220/1‑239 (MQ=255)                                                                                                                                                              gCGGATCATCTCGCCCTGTTCGGTTACGCTCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATAACCGCTTTTTCGCAGGTTTTGTTTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAACTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCaa                                                                             >  2:997947/1‑239 (MQ=255)                                                                                                                                                              gCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCCACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCCGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCaa                                                                             >  1:199068/1‑239 (MQ=255)                                                                                                                                                                ggATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCCGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCCGGAAGCTGCCATCACTCCCCCCTCTTTTGCTGAGTCGGAATAGCCAatc                                                                           >  2:107746/1‑239 (MQ=255)                                                                                                                                                                       cccccccTTTGCGTTTCCGCGCAGCCCCCCTTCCAGGCTCCCTGGGGTTTGGCACAGCAGCGCCGCTGGAG‑CA‑GGTCCCCCCCCGCCACAAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCa                                                                    <  2:512297/235‑1 (MQ=255)                                                                                                                                                                           cccTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCCCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGT‑ACCGCCGCGACCGTGGCACAACGTCAGCTCAATACCCGCTTTTTCGCAGGGTTTGATTAATGCTTCCTGTGCCTGATATT‑GCGCCCCGGACGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCa                                                                    >  2:606448/1‑235 (MQ=255)                                                                                                                                                                            ccTGTTGGGTCACCGGGAGCCCGCCTTTCGGGTCCCGGGGGGGTTGTGACAGCAGCGGCCAATGAG‑CA‑GGGGCGCCGCCGCGCCCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGt                                                               <  2:825669/239‑1 (MQ=255)                                                                                                                                                                            acTGTGCGGTTCCGCGCAGGCCGCTTTTGGGGTTCCGTGGGGGTGGGGACAGCAGCGCCGCAGGAG‑AA‑GGTGCGCCGCCGCAACCAAGGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGt                                                               <  1:111395/238‑1 (MQ=255)                                                                                                                                                                               gTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGttt                                                             >  1:933111/1‑238 (MQ=255)                                                                                                                                                                                ttCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTg                                                            >  1:1191112/1‑238 (MQ=255)                                                                                                                                                                                ttCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTg                                                            >  1:215190/1‑238 (MQ=255)                                                                                                                                                                                ttCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGGGCCTGATATT‑GCGCCCAGGAAGATGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTg                                                            >  1:256064/1‑238 (MQ=255)                                                                                                                                                                                        cgcgCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATc                                                   >  2:362766/1‑239 (MQ=255)                                                                                                                                                                                          cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAg                                                 >  1:890169/1‑239 (MQ=255)                                                                                                                                                                                          cgcAGGCCGCCTTTCAGGCTTCCTGGCGGGTGTGACAGCAGCGCCCCATGAG‑CA‑GGTGCGCCGCCGCGACCAAAGGA‑ACCGCCGCGACCGTGGAACAACGCCAGCTCCATACCCGCTTTTTCGCCGGTTTTTGTTAATGCATCCTGGGCCTGATATT‑GCGCCCAGCACGCTGCCATCACTCCCGCATCTTTTGCTCAGTTGGAGTAGCCCATCCTCACCCTCTGTTTGACCTGAATCa                                                  >  2:880596/1‑238 (MQ=255)                                                                                                                                                                                           ggaGGCCGCCTTTGAGGCTTCCTGGCGGTTGAGACAGCAGCGCCGCATGAG‑GA‑GGTGCGCCTCCGCGACAAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAg                                                 <  1:1025937/236‑1 (MQ=255)                                                                                                                                                                                               cccgccTTTCCGGCTCCCGGGCGTTTGGGCCGGCAGCGCGGAATGGG‑CA‑GGTGCGCCCCCGCGACCAATGGA‑CCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCa                                             <  2:215190/237‑1 (MQ=255)                                                                                                                                                                                                   ccTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGAt                                         >  1:1057350/1‑238 (MQ=255)                                                                                                                                                                                                       tCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCCATGGA‑ACCGCCCCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGGAGGTTTTGATTAATGCATCCTGTGCCCGATATT‑GCTCCCAGGAAACTGCCATCACTCCCGCATATTTTGCTGAGTCGGAAAAGCCCATCATCACCAACTGTTTGCCCTGAATCAGGCCCCGATAcac                                     >  2:704497/1‑236 (MQ=255)                                                                                                                                                                                                        ctggCTTCCGGGGGGTTGTGACAGAAGCGCCGCAGGAG‑AA‑GGTGCGCCGCCGCGACAAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGt                                   <  1:798146/237‑1 (MQ=255)                                                                                                                                                                                                        cgggCTTCGGGGCGTTGGGGACAGGCGCGCCGCATGAG‑GG‑GGGGCGCCGCCGCAACCAATGAA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGt                                   <  2:543359/237‑1 (MQ=255)                                                                                                                                                                                                              tCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAgg                                                >  1:117466/1‑220 (MQ=255)                                                                                                                                                                                                                       tgggaagagcggcGCCCAAGGAG‑CA‑GGTGCCCCGCCGCCCCAAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGagcagc                     <  1:704497/231‑1 (MQ=255)                                                                                                                                                                                                                        ggtgccggcagcCCCGCATAGG‑AG‑GGTGCCCCGCC‑CGACCAAAGGACACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCACAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCt                    <  2:806849/231‑1 (MQ=255)                                                                                                                                                                                                                          tgACCGCAGCCCCGCATGAG‑CA‑GGTGCGCCCCCGCGGCCAATGGA‑ACCGCCGGGCCCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGAGCCTGATATT‑GCGCCCAGGAAGCTGCCATCCATCGCGCATCCTTTGACGCGTCGGACTAGCCAATCATCATCCTCTGTTTGCCCTGAATCAGGCCACGATACTAGTCACTATTTAGCAGCCggg                 >  2:1010164/1‑239 (MQ=255)                                                                                                                                                                                                                                 agcGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCCTCACCATCTGTTTGCCCTGAATCAGGCCACGATCCCAGTCAATTTTGAGCAGCTTGGTCCTGAc          >  2:977495/1‑239 (MQ=255)                                                                                                                                                                                                                                      cgcAGGGG‑CA‑GGTGCCCCGCCGCGACCAATGAA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGt     <  1:107746/239‑1 (MQ=255)                                                                                                                                                                                                                                         aTGAG‑CG‑GGTGCGCCGCCGCGTCCAAGGGG‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTgg  <  1:575590/239‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                GCCGCGGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAA‑GCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAG‑CA‑GGTGCGCCGCCGCGACCAATGGA‑ACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATT‑GCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGG  >  NC_000913/4151089‑4151558

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