Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,829,873	T→G	20.8%	H298P (CAC→CCC)	astA ←	arginine succinyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,829,873	0	T	G	20.8%	141.3 / 3.4	72	H298P (CAC→CCC)	astA	arginine succinyltransferase
Reads supporting (aligned to +/- strand):  ref base T (10/43);  new base G (14/4);  total (25/47)
Fisher's exact test for biased strand distribution p-value = 1.20e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTC  >  NC_000913/1829639‑1830108                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               atTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CCAA‑AGCCCGGTTGCCGGGTCGGGACGCAC‑CAGGCCCCCCC‑GG‑CAATGGGGATa                                                                                                                                                                                                                                            >  1:257571/1‑239 (MQ=255)       ttgccATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑AGCTCGGTTGCCGGATAGGTAAGCCC‑CAGCACCACGC‑GG‑AAATGGGGATAAatt                                                                                                                                                                                                                                        >  1:549050/1‑237 (MQ=255)       ttgccATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAAT‑CATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTTAAAT‑CCAA‑AGCCCGGTTGCCGGATCGGGACGCAC‑CCGCCCCACCC‑GG‑AAAAGGGGAAAAatttta                                                                                                                                                                                                                                     >  2:438609/1‑237 (MQ=255)         gccATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCCCCGTTGGCGGATCGGTACGCCC‑CAGCACCCCGC‑GG‑AAATGGGGGAAATTTTCAt                                                                                                                                                                                                                                    >  1:266212/1‑239 (MQ=255)           cATAACACCTCGCCCGATACCGGATTACGCTTCACCCGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAACCCATAAAGTCATGCTGTTTGCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTg                                                                                                                                                                                                                                  <  1:547448/239‑1 (MQ=255)                    tCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTTAAAT‑CAAA‑CACCCGGTTGGCGGGCCGGGGCGCCC‑CAG‑GCCACACGGG‑ACACGGGGAGAAATTTCCTTGGCG‑ACcagg                                                                                                                                                                                                                        >  2:398689/1‑239 (MQ=255)                    tCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTACAATCCAAA‑CGC‑CGGGTGCCGGGACGGAACCGAC‑CACCCCCCCCC‑GG‑AAAAGGGGGAAAATTTATTTGAGG‑ACcagg                                                                                                                                                                                                                        >  2:454125/1‑239 (MQ=255)                    tCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGGTAAAAT‑CAAA‑CCCTCGGTTGCCGGGTCGGTACGCCC‑CAGCACCCCCC‑CG‑AAATGGGGGAAATTTTCATTGGCC‑ACcagg                                                                                                                                                                                                                        >  2:1017891/1‑239 (MQ=255)                    tCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCACCGGCGTGGCATTT‑GAGGGCATCCAGGTGTGCGGCGGTTACAAT‑CAAA‑CCACCGGGTGCCGGACTAGGAGGCAC‑CAACACCACGC‑GG‑AAAAGGTGAAAAATTTCATTGGCG‑aaaaaa                                                                                                                                                                                                                        >  2:530630/1‑234 (MQ=255)                    ccgaccGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTAACCAGTCCACGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACcagg                                                                                                                                                                                                                        <  1:522850/235‑1 (MQ=255)                    ccccccGATACCCGTTTGCCCTTCCCCCCCAATTCGGCCTGGGCCGTTATCCAGTCACCGTTAACCCATAAAGGCATGCGGTTTTCTCCCCTGCGCACAGGCGCACAAGACGAACGCGATCCCCGGCGTGGCATGT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCGCGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACcagg                                                                                                                                                                                                                        <  2:632515/236‑1 (MQ=255)                         cacatccGGTTTACGTCTCACAACCGGTGCCCCCTGGCCCCTTATCCGGGCACCGTTAATCCATAAAGTCATTCCGGTTTCCCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGGGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCaggcaggc                                                                                                                                                                                                                   <  2:317808/234‑1 (MQ=255)                           aTACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGGACGGAC‑CAGCACCACCC‑GGCAAAAGGGGA‑AAATTTCATTGGCG‑ACCCCGCAGGCTg                                                                                                                                                                                                                 >  1:692737/1‑239 (MQ=255)                            tACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑a                                                                                                                                                                                                                             <  1:898967/226‑1 (MQ=255)                                    tACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAAT‑CATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGcc                                                                                                                                                                                                      <  1:438609/239‑1 (MQ=255)                                     tcccttcccccgccttGCCCCCTGGCCCTTTTACCAGTCCCCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGcc                                                                                                                                                                                                      <  2:451564/224‑1 (MQ=255)                                      cGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCCCGGTTTGCGGGTCGGGACGGCC‑CAGCACCACCC‑GG‑AAATGGGGATAATTTTCATTTGCGCCCCAGGCGGGCGGGGaaacacc                                                                                                                                                                                                        >  2:764023/1‑234 (MQ=255)                                          ctcacgccgattgtccttGGCCCGTTTTCCCGTCACCGTTTATCCATAAAGTCATGCTGTTTTCCCCTCTGCGCACCGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTgcg                                                                                                                                                                                                 <  2:222245/222‑1 (MQ=255)                                            acccgcgATGCGCCCTGGCCCGTTACCCAGTCACCGTTAAACCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGAACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCa                                                                                                                                                                                               <  1:384371/236‑1 (MQ=255)                                                   aTGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACCC‑GG‑GAATGGGGATAATTTTCATTGGCG‑GCCAGGCAGGCTGGGGAAA‑CGCCCTGCCCAGGCTGc                                                                                                                                                                                         >  1:15759/1‑238 (MQ=255)                                                                      ccAGTCACCGTTAACCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAg                                                                                                                                                                     <  1:398689/239‑1 (MQ=255)                                                                        aGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CCGCACCACCC‑GG‑GAATGGGGATAATTTTCATTGGCG‑CCCAGGCAGGCTGGGGAAA‑CCCCCCCCCCACGCTGCCCTTCTTCCCCCCCCCCCAGcc                                                                                                                                                                   >  2:918281/1‑239 (MQ=255)                                                                          tccccttcacgcctaaaaGTCTGGCTGTTTTCGCCGCGGGGCACAGGGGCACCAGAAGAACGCGGTCCCCGGGGGGGCATTT‑GAGGGCATACAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑a                                                                                                                                                                <  2:654170/224‑1 (MQ=255)                                                                           ccccGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGcc                                                                                                                                                                   <  1:554840/234‑1 (MQ=255)                                                                           cACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGGCGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑Ac                                                                                                                                                               <  1:12869/239‑1 (MQ=255)                                                                            aCCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTTAAAT‑AAAA‑CCCCCGGTTGCCGGATCGGTACGCCC‑CCGCCCCCCCC‑GG‑CAAAGGTGATAATTTTCCTTGGCG‑ACCCGGCAGGCTCGGAAAA‑ACCCCCGCCCAGGCTGCCCTTCTTCCACTTccaccccccccct                                                                                                                                                               >  2:952515/1‑231 (MQ=255)                                                                            aCCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCCCGGTTGCCGGATCGGGACGCAC‑CAGCACCACGC‑GG‑AAATGGGGGTAATTTTCATTGGCG‑ACCAGGCCGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAACCG‑a                                                                                                                                                                >  1:198916/1‑237 (MQ=255)                                                                                  aaTCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCCC‑CAGCACCACGC‑GG‑AAATGGTGGTAATTTTGATTGGCA‑ACCAAGAGAGATGGGAAAT‑CACCCCGCCCAGGCTGCCCTCTTTCCACTCCCACCACCCGCACtttaag                                                                                                                                                         >  2:533138/1‑236 (MQ=255)                                                                                   aacccaaaaggcatctttttttCCCCTCTGCGCCGCGGCCCCCCGGACGAACGGGACCCCCGGCGTGGC‑TGT‑GGGGGCATCCGGTTGGGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCGCCACGC‑GG‑AAATGGGGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGa                                                                                                                                                      <  2:940401/221‑1 (MQ=255)                                                                                           aGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTgc                                                                                                                                                                                                  <  1:610650/189‑1 (MQ=255)                                                                                               aTGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGGTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑GAATGGGGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAAGCTGCCCTTTTGCCCCTTCCACCCCCCG‑CCTTTTACGGagtgggcgcccc                                                                                                                                           >  1:754072/1‑228 (MQ=255)                                                                                               aTGCTGTTTTCCCCTCTGCGCCAAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGGGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCAcg                                                                                                                                           <  1:518506/239‑1 (MQ=255)                                                                                                   tttttccccctcTGCGCCCAGGCCCACCAAACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGt                                                                                                                                       <  1:454125/231‑1 (MQ=255)                                                                                                       ttCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCCCGGTTGCCGGCTCGGTACGCAC‑CCGCACCACGC‑GG‑GAATGGGGGGAATTTTCATTGGGG‑GACAGGCAGGATGGGAAAT‑CGCCCTGCGCAGGCGCCCCTTCTGCCACTTACAACACCCG‑CCCTTTAAAGAAGGCACACACACGGGCggg                                                                                                                                   >  2:429778/1‑237 (MQ=255)                                                                                                        tcccctcTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGGGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATg                                                                                                                                  <  1:421747/236‑1 (MQ=255)                                                                                                                         cgaaCCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGt                                                                                                                    <  1:255301/233‑1 (MQ=255)                                                                                                                           ccccAGACGAACGCGTTCCCCGGCGTGGCTTTT‑GGGGGCACCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGcc                                                                                                               <  2:177944/237‑1 (MQ=255)                                                                                                                           cACCAGCCGAACGCGATCCCCGGGGTGGCAGTT‑GAGGGCATCCAGTTGTGCGGAGGTTAAAAA‑CAAA‑CGCTCGGTTGCCGGAGCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGcc                                                                                                               <  1:789794/239‑1 (MQ=255)                                                                                                                           cACCAGACGAACCGTAGCCCCGGCGTGGCATTG‑GGGGGCATCCAGTTGGGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCCCCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCCCCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGcc                                                                                                               <  2:692737/239‑1 (MQ=255)                                                                                                                                   gaagcggacccccGGCGTGGCATTT‑GGGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCa                                                                                                       <  1:732632/230‑1 (MQ=255)                                                                                                                                     accgctcccccGGGGTTGCATTTGGAGGACA‑CCCGGTGGGGGGCGGTTAAAAT‑CAAACCGCTGGGTT‑CCGGATCGGGCCGCCC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCaaa                                                                                                     <  2:254912/232‑1 (MQ=255)                                                                                                                                      cgcgATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGGACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCCCTTCCACCAGCCG‑ACTTTTTAGGATGGCGGGCACCCGGGCGATGTCACCCTCAAACGTCGGCCCCCCGCCAAAg                                                                                                    >  1:698693/1‑239 (MQ=255)                                                                                                                                           tCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGGGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGt                                                                                                <  1:533138/238‑1 (MQ=255)                                                                                                                                             cccGGGGTGGATTTG‑GGGGGTACCCAGTTGTGGGGCGGTTAAAAT‑CAAA‑GGCTCGGTTGCCGGGCCGGTACGAAC‑CAGCACCACGG‑GG‑AAATGGTGATAATTTTCATTGCCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCATTTCCACCAGCAG‑ACTTTTACGGATGGCGCGCACGCGGTAGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGa                                                                                             <  2:991958/239‑1 (MQ=255)                                                                                                                                              ccGGGGTGGCTTTT‑GGGGGCATCCAGTTTGGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTCCGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGat                                                                                            <  2:698693/239‑1 (MQ=255)                                                                                                                                              ccGGCGTGGCATTT‑GGGGGGATCCAGTTGTGCCGCGGTTAAAAT‑CAAC‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGATCCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGat                                                                                            <  1:530630/239‑1 (MQ=255)                                                                                                                                              ccGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGat                                                                                            <  1:875448/239‑1 (MQ=255)                                                                                                                                              ccGGCGGGGCATTT‑GAGGGCATCCAGTTGGGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGat                                                                                            <  1:1003933/239‑1 (MQ=255)                                                                                                                                                                      gcggttggaggggtttAAAAT‑CAAG‑CGCGCGGTTCGGCGATGGGTACCCCC‑CAGCACCAAGG‑GG‑AAATGGGGATAATTTTCATGGGCG‑ACCAGGAAGGCAGGGAAAG‑CGCCCTGCGCAGGCTGCCCTTCGGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACctt                                                                     <  2:32761/225‑1 (MQ=255)                                                                                                                                                                        gtttggcggggTTTAAAAT‑GAAC‑TGCTCCGTTGCCGGGCCGGTACCGAC‑CAGCCCCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACcttct                                                                   <  2:15759/234‑1 (MQ=255)                                                                                                                                                                        gttgggggggtgTTAAAAG‑AAAA‑CGCGGCGTTGCCGGGTGGGGACCCACACGGCCCAACG‑‑GG‑AAAGTGTGGTAATTTTCATTGGGC‑GCCAGGCCGGGTGGGAAATGCGCCTGGCCGAGG‑TGGCCTTCCGCCCCCTCCCCCCGGCG‑GCCTTTACCGATGGCGCGCACGGGGTCGATGGCACAATCAAGCGGCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACcttct                                                                   <  2:988433/228‑1 (MQ=255)                                                                                                                                                                        ggtgggcggcgGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCGC‑CAGCCCCACGC‑GG‑AAATGGTGATAATTTTCATTTGCG‑ACCAGGCAGGCTGGGAAAG‑CGCCCTGGGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACcttct                                                                   <  2:823743/234‑1 (MQ=255)                                                                                                                                                                        gTTGTGCGGGGGTTAAAAT‑CAAA‑CGCTTGGTTGCCGGAACGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACcttct                                                                   <  2:266212/239‑1 (MQ=255)                                                                                                                                                                         ttGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCCCGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACcttct                                                                   <  1:377821/238‑1 (MQ=255)                                                                                                                                                                                  ggTGAAAAA‑CAAA‑CGCTCGGTGGCCGGACCGGTACGCAC‑CAGCCCCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGACAT‑CGCCCTGGGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGc                                                          <  2:44759/238‑1 (MQ=255)                                                                                                                                                                                                    tCGGTTGCCGGATCGGTACGCAC‑CAGCCCCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGgcg                                          <  1:764023/238‑1 (MQ=255)                                                                                                                                                                                                      ggTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGt                                                                                                <  1:406499/182‑1 (MQ=255)                                                                                                                                                                                                      ggTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTAGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATgc                                                                                                >  2:406499/1‑181 (MQ=255)                                                                                                                                                                                                          gCCGGATCGGTACGCAC‑CAGCCCCACGC‑GG‑AAATGGTGATAATTTGCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTg                                   <  1:93265/239‑1 (MQ=255)                                                                                                                                                                                                            cGGATGGGTCCGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATGGGCG‑ACCAGGCAGGCTGGGAAAC‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGc                                  <  2:788316/238‑1 (MQ=255)                                                                                                                                                                                                              ggcgccgtccgcccccaccccccgc‑gG‑AACTGGGGATTATTTTCCTTGGGG‑GCCAGGCCGGCCGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCCCCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCgg                                <  2:667066/217‑1 (MQ=255)                                                                                                                                                                                                                  ggTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCCCGGGCAGGCGCGGTTTGCGGATGt                            >  1:44636/1‑238 (MQ=255)                                                                                                                                                                                                                         ac‑cagcaccaCGC‑GG‑AAATGGCGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGCCGATGTCACACTCAAGCGTCGGCCCACCGCCAAAGATTGCGACATAATTACGGGAGCGAAAACCTACTTTCCCCAGCACCCCCCGGGCAGGCGCGGCTTGCgagagaccccaacc                    >  2:456200/1‑226 (MQ=255)                                                                                                                                                                                                                              gccccccgc‑gG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATg                <  1:102012/233‑1 (MQ=255)                                                                                                                                                                                                                              gcaccaCGC‑GG‑AAATGGGGATAATTTTCAGTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATg                <  1:366760/239‑1 (MQ=255)                                                                                                                                                                                                                                    cgc‑gG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCCCACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTCCTTTCTCCCGCACCCCCCAGGAAGACGCGGTTTTCGGGGTTACCACGCCGATGAGGTcc          >  2:620713/1‑239 (MQ=255)                                                                                                                                                                                                                                    cgc‑gG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGCTATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTcc          >  1:39196/1‑239 (MQ=255)                                                                                                                                                                                                                                    cgc‑gG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGGGAAAACCTTCTTTCTCCACCACCGCGCGGGCAGGCGCGGGTTGCGGGTGTACCTGAcccaagaagccc          >  2:102012/1‑229 (MQ=255)                                                                                                                                                                                                                                           aaaTGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTggg      <  1:952515/238‑1 (MQ=255)                                                                                                                                                                                                                                           aaaTGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCCCGGGCAGGCGCGGGTTGCGGATGTACCGGACCGATGACGTCCTggg      >  1:583461/1‑238 (MQ=255)                                                                                                                                                                                                                                              tGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTc  <  1:918281/239‑1 (MQ=255)                                                                                                                                                                                                                                               ggTGATAATTTTCATTGGCG‑CCCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCCGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGACACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATg                <  1:533563/224‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               ATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCCCGTTATCCAGTCACCGTTAATCCATAAAGTCATGCTGTTTTCTCCTCTGCGCACAGGCGCACCAGACGAACGCGATCCCCGGCGTGGCATTT‑GAGGGCATCCAGTTGTGCGGCGGTTAAAAT‑CAAA‑CGCTCGGTTGCCGGATCGGTACGCAC‑CAGCACCACGC‑GG‑AAATGGTGATAATTTTCATTGGCG‑ACCAGGCAGGCTGGGAAAT‑CGCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCG‑ACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTC  >  NC_000913/1829639‑1830108

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