Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,239,855	A→G	18.0%	E326G (GAG→GGG)	dadX →	alanine racemase, catabolic, PLP‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,239,855	0	A	G	18.0%	128.6 / 3.0	65	E326G (GAG→GGG)	dadX	alanine racemase, catabolic, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base A (18/34);  new base G (9/4);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 3.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGCGATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGG  >  NC_000913/1239653‑1240091                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   gcgcgcGATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTCAGCCGGGTGGGCCGT‑GGGGCAAGGAGGTCACAAATGGATG‑TGTcggcgcccgcgccggcaa                                                                                                                                                                                                               >  1:894877/1‑228 (MQ=255)  cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGT‑GGGGGAAGGAGGTCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACg                                                                                                                                                                                                              >  1:365544/1‑239 (MQ=255)  cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCCGT‑TGGGGCAGGGGATCA‑AAAATGTTGATTGTGCCGCCCCGGCCGGaaag                                                                                                                                                                                                              >  1:776991/1‑237 (MQ=255)  cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGGTTGGCTGT‑GGGGCAAGGAGGTCA‑AAAATGATGATGTCTCCGCCGCTGCCGGGACg                                                                                                                                                                                                              >  1:235213/1‑239 (MQ=255)  cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTAAGCCGGTTTGGCTGT‑GGGGGAAGGGGGTCA‑CAATTTTTGATTTGCCCGCCCCTGCCGGGACg                                                                                                                                                                                                              >  1:1228248/1‑239 (MQ=255)  cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGAAGCCGGTTGAGCCTG‑GGGGCACGGGGGTCA‑AAATTGATGATGTGGCCGCCGCGGCCGGaaac                                                                                                                                                                                                              >  1:158986/1‑237 (MQ=255)    aggagacagcAGGGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGt                                                                                                                                                                                                            <  1:1041484/230‑1 (MQ=255)     ggcagtaacaggtgaTCGGCATGGTGGCCGCGGGGTACGCCGACGGTTTGCCGCGCCACGGG‑CCTACCGTTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGt                                                                                                                                                                                                            <  1:1022113/224‑1 (MQ=255)     gcgATTAACCGGCAATCGGCATTGGCGCCGCAGGGTACGCCCGCGGTTATCCGCGCCACGCG‑GCTGCCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGt                                                                                                                                                                                                            <  1:1267717/238‑1 (MQ=255)       gATGAACCGCGAATGCGCATTTTCCCCCCGGGGTACGCCGACGGTTACCCGGCCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTggg                                                                                                                                                                                                         <  1:1301635/239‑1 (MQ=255)          gAACAGCGAATCGGCATGGCCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCt                                                                                                                                                                                                       <  1:967194/238‑1 (MQ=255)           aaCAGGGAATCGGCATTGGCGCCGGAGGGTACGCAGACGGTTAGCCGCGCCACGCG‑CCTACCGGTACCCCGGTTTTAGTGGACGGCGTGCGCCCCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCt                                                                                                                                                                                                       <  1:844603/237‑1 (MQ=255)            ccaccgatgcggcgttttgcccgcTGTGTCCGCCGCGGGTTAGCCGGGCCACGGG‑CTTACCGGTCCCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCGCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTAt                                                                                                                                                                                                     <  2:1307344/218‑1 (MQ=255)                  attCGGCATTGTCCCCGCTGGGTCCGGCGCGGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGGTGGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTg                                                                                                                                                                                               <  1:541813/236‑1 (MQ=255)                          tggTCGCCGCAGGGTTCGCCGGCGGTTACCCGCGCCACGGG‑CCGACCGGTACCCGTGTTTTAGGGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑gcgcg                                                                                                                                                                                    <  1:536608/237‑1 (MQ=255)                           tggcgccgcAGGGTACGGCGACGGTTATCCGGGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑gcgcg                                                                                                                                                                                    <  1:1173024/235‑1 (MQ=255)                                          gcgcGACGGTTACCGGCCGCAGGCG‑CGTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTc                                                                                                                                                                     <  1:298652/235‑1 (MQ=255)                                           gCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTCTTGGGACCGCGGTTGGGCGGTGGGGGCAGGGGGACAA‑AAATTG‑TTATGGCGCCGCCGCTCCCGCACCGGGGGGGTATTAAATGA‑TGTGgcgcgcgcgggaacagggcac                                                                                                                                                                    >  2:84774/1‑224 (MQ=255)                                             ccggggttgcccgcgcc‑cgcgcCCTACGGGTCCCCCTGGTTTAGGGGGCGGCGGGGCGACCATGACGGTGGGGACCGTCTCGGTGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTTTTGGTACGCCGGTTGAGCTGG‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCgg                                                                                                                                                                 <  2:22239/229‑1 (MQ=255)                                                    tacccGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGGGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGa                                                                                                                                                           <  1:842908/235‑1 (MQ=255)                                                      tcccccgcccccccctttccGTTACCCCTGTTTTGGTGGAGGGCGTGGGCGCCAGGGGGGGGGGGATGGTCTCGATGGATTGGCTAGCGGGCGATTTAACGCCTTGCCCGCAGGGGGTTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACg                                                                                                                                                         <  2:667535/221‑1 (MQ=255)                                                      tCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGAGGTGCCCGCCGCTGCCGGGACGGTGGGCGTTGAGGTTA‑TGT‑GGGCGCTGGGGCTACGCGTCCCGGGTGTGGCgg                                                                                                                                                        >  1:371455/1‑239 (MQ=255)                                                                                 cTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGGGATCA‑AAATTGATGAAGTCGCCGCCGCTGCCGGGACGGTGGGGTATGAGTTGA‑TGT‑GGGCGCTGGCGCTACCCGTCCCGGTTTGGACGGGGTTAC‑TTGTTGTTAGGCCGGTCggg                                                                                                                             >  1:547168/1‑238 (MQ=255)                                                                                      ttggTGGGCGGCGTGGGCCCCATGGCGGTGGGGGCCGGCTCGATGGATACGCTAGCGGTCGATTTACGGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGGGCGGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCaa                                                                                                                        <  2:1069939/236‑1 (MQ=255)                                                                                              cGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGGACGGGGGGCTATGAGTTGA‑TGT‑GCGGGCTGGCGCTACGCGTCCCGGTTGGGACGGGGTAAC‑TTGTTTTTAGACGGGTCGGCGGCAACGgct                                                                                                                   >  1:824201/1‑236 (MQ=255)                                                                                                           aTGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGATACGCCGGGTGAGCTGT‑GGGGCCAGGAGATCC‑AAATTGGTGGTTGCGCCGCCGCTGCCGGGACGGGGGGCTATGAGGTTA‑TGT‑TCGCGCTGGGGCTACGCGTCCCCGTTGTGACGGGGTAAC‑TTTTTTTAAACCGGATCGGAGGCAACACCATCTGGGGGGaaaaaaa                                                                                                   >  2:533299/1‑239 (MQ=255)                                                                                                                  tGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGG‑TGGGCAAGGGGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGGACGGGGGGGTATGTGGTGT‑TTT‑GCGGGCTGGCGCTACGCGTCCCGGGTGTGGCGGGGGAAA‑TTGTTTTTATGCGGGACGGGGGGCACGGCTTCTGGGTGCaaaaaaaaaaaccc                                                                                            >  2:23635/1‑232 (MQ=255)                                                                                                                  tGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCCGGTGAGCTGG‑GGGGCAAGGAGAGCA‑AAATTGTTGTTGTCGCCGCCGCTGCCGGGAAGGAGGGCTAGGAGTTTATTGT‑CCG‑GCTGGGCCTAGCCGGCCCGGGTGGGACGGGGGGAC‑TTTTTGTTAGCCGGGCTCGAGGCAACCGCTTTTGGGGGGGAAAAAAAaatcc                                                                                             >  2:957579/1‑238 (MQ=255)                                                                                                                   ggggCCCGTCGCGATGGATATGCTAGGGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATcatcca                                                                                            <  1:796644/238‑1 (MQ=255)                                                                                                                    gggACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACCCCGGGTGAGCTGT‑GGGGCAAAGGGATCC‑AAATTGTTGATGTCGCCGCCCCTGCCGGGAGGGGGGGCTACGAGTTTA‑TTG‑GCGCGCTGGCGTGACGCGGCCCGGTGGGTGCGGGGTTACTTTGTTGTACGCGGGGTGGGGGTCACCGCTCTCTGGGGACaaaaaaacaaaccac                                                                                           >  2:983549/1‑231 (MQ=255)                                                                                                                    gggACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTAAGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGGTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGGCCCGGGTGTGACGGGGTAAA‑TTTTTGTAAACCGGATCTGAGGCAACGTCTTCTTGGGGCAAAAAAAAcatccat                                                                                           >  1:500059/1‑238 (MQ=255)                                                                                                                    gggACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGACGCCGGTTGAGGTGT‑GGGGCCAGGAGATCA‑AAATTGGTGTTGTCGCCGCCGCCGCCGGAACGGGGGGGTATGGGGTGG‑TGG‑TCACGCTGGCGCCCCCGGCCCCGGTTGGTGCGGGGTAAC‑TTTTTTGAAACCGGGGCGGGGGCACGGCGTTTTGGGGGCaaaaaaaaaataccc                                                                                           >  2:1181027/1‑230 (MQ=255)                                                                                                                        ccGTCTCGATGGATTTGCTGGCGGTCGGTTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCgg                                                                                       <  1:15699/238‑1 (MQ=255)                                                                                                                           tctcGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCCAGGGGAACA‑AAATTGATGATGTCCCCGCCGCTGCCCGGAGGGGGGGCCTTTAGTTGA‑GGT‑GGGGGCTTGGGCTACGCCTCCCGGTGGTGTCGGGGTACC‑TTTTTGTAAGGCGGGACGGGGGGGAGGGCGTCTTGGGGCaaaaaaaaacacaaaccgggggg                                                                                   >  2:922797/1‑223 (MQ=255)                                                                                                                           tctcGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAAATGATTATTTCCCCGCCCCCGCCGGAACGGGGGGGGTTGGGTTTA‑TGT‑GCGCCCTGGCGCCACCCGGCCCCGTGGTGGCGGGGTAAA‑TTTTTTTAAGCCGGGCCGGGGGCCAAGTTTTTCGGGGGGaaaaaaaaaaaacaacagggggg                                                                                   >  2:1205255/1‑223 (MQ=255)                                                                                                                             tcGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTc                                                                                <  1:922797/239‑1 (MQ=255)                                                                                                                                     tctGCGAGCGGTCGATTTAACGCCTTGCCCGCAGGCGTGTATTGGGACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACt                                                                        <  1:772204/236‑1 (MQ=255)                                                                                                                                        gCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGGTGAGCTGT‑GGGGCAAGGGGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGGGGGCTGTGAGTTGA‑TGT‑GCGCGCTGGGGCTACGCGTCCCGGGTGGGCAGGGGTGAC‑TTTTTTTAAACGGGAGCGGGGGGCACGGCTCTTGGGGGCAAAAAAAACATACATACCGC‑TGGGGAGCA‑ACGGtgg                                                                    >  2:905932/1‑237 (MQ=255)                                                                                                                                            gCGGTCGATTTAACGCCTTGCCCGCAGGGGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAgttgt                                                                <  1:319857/239‑1 (MQ=255)                                                                                                                                                      tAACGGCTTGCCCGGAGGCGGGTATTGGTACGCCGGTTGAGCTGG‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACa                                                      <  1:14272/239‑1 (MQ=255)                                                                                                                                                            ctgggcgggaGGCGGGTATTGGGACGCCGGTTGAGCTGG‑GGGGCACGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGGTGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGGGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCa                                                <  1:1241375/230‑1 (MQ=255)                                                                                                                                                            cTTGCCCGGGGGGGGGTGGTGGTGCGGCGGTTGGGCTGG‑GGGGCAAGGAGGTCA‑AAATTGATGATGGCGCCGCCGCGGCCGGGACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCa                                                <  2:894877/239‑1 (MQ=255)                                                                                                                                                            cTTGCCCGCGGGGGGTTTTTGGTACGGGGTTTGGGCTGG‑GGGGAAGGGGGATAA‑AAATGGATGATGTCGCCGCCGCTGCCGGAACGGTGGGGTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGGGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCa                                                <  2:371455/239‑1 (MQ=255)                                                                                                                                                            cTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCa                                                <  1:905932/239‑1 (MQ=255)                                                                                                                                                            cTTGCCCCGAGGGGGGTTTTGGTACGCGGGTTGAGCTGG‑GGGGCAAGGGGATAA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCa                                                <  2:824201/239‑1 (MQ=255)                                                                                                                                                             ttGCCGGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGGCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCa                                                <  1:1205255/238‑1 (MQ=255)                                                                                                                                                              tGCCCGCGGGCGGGTTTTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCAtt                                              <  1:23635/239‑1 (MQ=255)                                                                                                                                                              tGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCCCCGCTGCCGGGACGGTGGGCTATGAGTTTA‑TGT‑GCGCGCTGTCGCCACGCGTCCCGGTGTTGGCGGGGGAAC‑TTGTTGTAAACCGGGTCGCGGGCAACGTCCTTTGGGTGCAAAAAAAACATCCATCCCGCGTGGGC‑GCA‑AACGTAGTTGTTATtggtaaaaagccaa                                               >  2:730382/1‑226 (MQ=255)                                                                                                                                                                 ccGCGGGGGGGTTTTGGGTCGCCGGTTGGGCGGG‑GGGGCAGGGGGATCA‑AAATTGATGGTGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGc                                            <  2:776991/238‑1 (MQ=255)                                                                                                                                                                 caggaGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGcc                                           <  1:324321/235‑1 (MQ=255)                                                                                                                                                                   gCAGGGGGGTATTGGTACGCCGGTTGAGGTGT‑GGGGCAAGTAGATCA‑AAATTGTTGCTGGCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCa                                         <  1:957579/239‑1 (MQ=255)                                                                                                                                                                                    cGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAATGTCCATTaaa                        >  2:223807/1‑239 (MQ=255)                                                                                                                                                                                     gCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTTTAAGCCGGATCGGAGGCACCGCCTTTGGGGGGCAAAAAAATCTTCAACCCGG‑‑TGGGCCGTACACGGAGGTTTTTAAGGTGAAAGAGCAATTGCCCATGGTATAGTCCATAaaaa                       >  2:1371628/1‑239 (MQ=255)                                                                                                                                                                                     gCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAAAGCCATTGCCCATGATAGTGTACATTaaaa                       >  1:529056/1‑239 (MQ=255)                                                                                                                                                                                     gCCGGTTGAGCTGT‑GGGGCAAGGAGATAA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTaaa                        <  1:1244113/238‑1 (MQ=255)                                                                                                                                                                                      ccGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATAATATGGTCCATTAAAAg                      >  2:1359484/1‑239 (MQ=255)                                                                                                                                                                                      ccGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATAATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGCAAGGGCCATTGCCCCTGATAGTGGCCATTAAAAg                      >  2:1244113/1‑239 (MQ=255)                                                                                                                                                                                         gTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGCGTTGA‑TAT‑GGGCGCTGGCGCTACGCGGCCCGGTTGTGACGGTGTAAC‑TTGTGGTAAGCCGGATCGGCGGCACCGTCTTCTGGGTGCAAAAAAACCTTAAATCCGGC‑TGTTCAGCA‑ACGGTAGTTGTTAATGTGCCACAACCATTGCCCAtaaagagctccttaaaaaggat                   >  2:1150135/1‑218 (MQ=255)                                                                                                                                                                                          tttggtgtt‑gGGGAAAGGAGTACA‑AAATTGGTGGGGTCGCCGCCGCTTCCGGAACGGTGGGTTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGGGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATTCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATg                  <  2:547168/231‑1 (MQ=255)                                                                                                                                                                                            ggggttt‑gGGGCAAGGGATTTA‑AATTTGATGGGTTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAA‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGa                <  2:235213/233‑1 (MQ=255)                                                                                                                                                                                            gAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGa                <  1:223807/239‑1 (MQ=255)                                                                                                                                                                                              gCTGT‑GGGGCAAGGAGATAA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGGCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCAGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGac               <  1:310620/238‑1 (MQ=255)                                                                                                                                                                                               cTGT‑GGGGCAAGGAGATCA‑AAATTGAGGATGGCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGaca              <  1:553599/238‑1 (MQ=255)                                                                                                                                                                                                  gtcgggCAACGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGGGCAAAAAAATCATCCATCCGGC‑TGGCGAGCA‑ACTGTAGTTGGTAATTGGACAGAACCATTGCCCATGACAAGGGCCATTAAAAGGAAGGgcacaa           >  2:654893/4‑237 (MQ=255)                                                                                                                                                                                                          aGGAAAACA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTcccc    <  1:654893/238‑1 (MQ=255)                                                                                                                                                                                                             agaTCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑CGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCgg  <  1:429192/237‑1 (MQ=255)                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   GCGCGCGATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGT‑GGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGT‑GCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGC‑TGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGG  >  NC_000913/1239653‑1240091

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