Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,406,283	T→G	21.5%	V102G (GTG→GGG)	smrA →	DNA endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,406,283	0	T	G	21.5%	128.4 / 6.9	63	V102G (GTG→GGG)	smrA	DNA endonuclease
Reads supporting (aligned to +/- strand):  ref base T (17/29);  new base G (12/4);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGGCGCAGGCGAAGCAAGAAAACTGGGAG  >  NC_000913/1406048‑1406518                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             gATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGtt                                                                                                                                                                                                                                           <  2:848839/239‑1 (MQ=255)             cATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGTAAGGGGTGCAACAATGGGTGCCGCGAAAG‑CCGCGCCGTGGGAAATATCCGCAACAGGCG‑GGCCCGAATCTTTTGCCCCCACCGGTAGAAAATGGACCAACAAAGGGGTTGAGTTTTAttga                                                                                                                                                                                                                               >  1:436456/1‑237 (MQ=255)              aTCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGGGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGGGCCGCAAAATGGGGGTACGTTTTATTca                                                                                                                                                                                                                               >  2:595424/1‑238 (MQ=255)               ttccaCGCGTAACCAACGTTCGCCGGAGCGTTTCAACCCGCGGCAGGTTGAAAGTTTCCTCCCCCGCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTcaac                                                                                                                                                                                                                             <  1:442730/235‑1 (MQ=255)                    cgcgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGGCTGAATCTTTTGCGCCAGCCGGGGGAAGAGGGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGc                                                                                                                                                                                                                         >  2:575012/1‑238 (MQ=255)                    actgTAACCACCGTTCCCCGCACGGTTCCGCCCCGCTGCAGTTTGACAATTCCCTCACCACCGGGTTTCTCG‑AAATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGcc                                                                                                                                                                                                                        <  1:155512/236‑1 (MQ=255)                           ccAACGTGCGCCGCAGCGTATCGACCCGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCgg                                                                                                                                                                                                                 <  2:496579/239‑1 (MQ=255)                             aaCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGAt                                                                                                                                                                                                               <  2:848698/239‑1 (MQ=255)                                gTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAAAATGGGGTGCTGGATAAA‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGAgggg                                                                                                                                                                                                            >  2:1154449/1‑238 (MQ=255)                                   cgccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGGGTGCCGCAAAATGGGGTTCAGTTTTATTTAACAAGCCCTGGCGGAGGGtttt                                                                                                                                                                                                         >  2:1098790/1‑238 (MQ=255)                                   cgccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGGAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGGGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTTAACAAGCCCTGGCGGATggcgg                                                                                                                                                                                                          >  2:155512/1‑235 (MQ=255)                                              atctccacGCTGCATGTTGATAATTTCCTCACCGCCGGATTTCGCG‑ACATCCTCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCt                                                                                                                                                                                             <  2:1016735/234‑1 (MQ=255)                                                     ccccggaggtttaataTTTCCTCCCCCCGTGATTTCTCG‑ACATCATCCCACTAAGTCCGCCGCGGGAGTTGCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGattatt                                                                                                                                                                                      <  1:575012/224‑1 (MQ=255)                                                             cTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGGTTGCG‑TAACGTGCTGATTATTCATGGTaa                                                                                                                                                                              >  2:880039/1‑239 (MQ=255)                                                                   aaTTTCCTCTCCCCCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCg                                                                                                                                                                        <  2:432322/239‑1 (MQ=255)                                                                         cTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGcc                                                                                                                                                                                                                        <  2:151472/186‑1 (MQ=255)                                                                              caccGGATTTGTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCa                                                                                                                                                                                    <  2:725800/216‑1 (MQ=255)                                                                                ccGTATTTCTCGTACA‑CATCCTACTTAGTCAGCGCCTGTGGTTTGGGCGGGAATGGTTATAGAAGCGGGGGGTGGATAAA‑CGAAGAAGGGAAAAACAGCCGCAAGAGACG‑GTCCTGACTCTTTTGCGCCAGCCGGTGGAAGAAAGCCGCAACATGGGGTTCAGTTTTATTCAACAAGCCATGGCGGATGGTGTGCG‑TAACTTGCTGATTATTCATGGTAAAGGGCGGGATGATaaagct                                                                                                                                                           <  1:1013431/239‑4 (MQ=255)                                                                                ccGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCCAAATGGGGTTCAGGTTTATTCAACAAGCCCTGGCGGATGGTTTGGG‑TAACGTGCTGGTTTTTCATGGTAAAGGGCGGGATGGTAAATc                                                                                                                                                            >  2:896413/1‑238 (MQ=255)                                                                                ccGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGGGGAAGAGAGCCGCAAAATGGGGTTCCGTTTTATTCCACAAGCCCTGGCGGATGGTTTTCG‑TAACGTTCTGATTTTTTATGGTAAAGGGCGGGGTGATAAAtg                                                                                                                                                            >  2:1070824/1‑237 (MQ=255)                                                                                ccGGAGTCCTCG‑AATTCACCCCACAACGGCAGCCTGTGGGGTTTCGGGGGGAAGGGTTGAAACATGGGGTGCTGGATAAG‑CTGCGCAGGGGTAAATATCCGCAACAGGCG‑AGGCGGAATCTTTTGCGCCAGCAGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCg                                                                                                                                                           <  1:809305/239‑1 (MQ=255)                                                                                acGGATTTCCCG‑GCAACATCCCACTAAGGCAGGCGCTTGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCg                                                                                                                                                           <  2:127949/238‑1 (MQ=255)                                                                                 cgcgagttttcgactgcgtcCCTCTAAGTCGGCCGCGGGGGTTACGGCGGGAAGGGTTGCGAATTGGTGGCGTACAAAAC‑CTGCCCGGTGGTAGATAGCGTCAA‑AGGTGAAGCCTGAACCTTTGGCGCCAGCCGGTGGAAGAGGGCCGCAAAATGGTGTTCAGTTTTATTGAACAAGCCCTGGCGGATGGTTTGCG‑TAATGTGCTGATTATTCATGATAAAGGGCGGGAGGATAAAGCg                                                                                                                                                           <  1:702131/221‑1 (MQ=255)                                                                                        ctcG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGATGCCACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGGAACAGGCG‑TGCCCGACTCTTTTGCGCCGGCCGGGGGTAGAGTGCCCCAAACTGGGGTTCAGCTTTATTCTAAAAGCCCTGGCGGCTAGTTTGCGCAAAGGTGCT‑ATTATTCATTGTAAAAAGCGGGATGGTAAAACgacagccaa                                                                                                                                                   >  1:921802/1‑231 (MQ=255)                                                                                               atcatcCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGGGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCCGGCGGATGGGTTTCG‑TAACGTGCTGATTTTTCCTGGTAAAGGGCGGGGTGATAAATCGCATGCCCATATTg                                                                                                                                              >  2:594154/1‑238 (MQ=255)                                                                                                tcatcCCAGGAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTTCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTc                                                                                                                                            <  2:93180/239‑1 (MQ=255)                                                                                                  gccccaccaattcggcgtgtGGGGTTGCGGCGGGAAGGGTTGCAACGTGGGGTGCTGGTTAAG‑CTGCGCAGGGGTAAATCCCGGCAACAGGCT‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCg                                                                                                                                          <  1:851246/220‑1 (MQ=255)                                                                                                  atcCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTTCGCCAGCCGGTGGAAGAGAGCCGCAAAATGGGGTTCAGTTTTATTCAAAAAGCCCTGCCGGAGGGTTTGTG‑GAACGTGCTGATTATTTCTTGTAAAAGGCGGGATGATAAATGGCATGCCAATAtt                                                                                                                                               >  2:661247/1‑234 (MQ=255)                                                                                                        cTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGGAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGGTTGCG‑GAACGTGCTGGTTATTCATGGTAAAGGTCGGGATGGTAAATGGCATGCCAAAATTTTccccgacta                                                                                                                                    >  2:410749/1‑232 (MQ=255)                                                                                                                 ccGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGGGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGGTTATTCATGGTAAAGGGCGGGATGATAAATc                                                                                                                                                            >  2:984068/1‑206 (MQ=255)                                                                                                                         gtgtCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCgctggctg                                                                                                                    <  2:921802/235‑1 (MQ=255)                                                                                                                           ttCGGGGGGAAGGGTTGCAACAGGGGGGGCTGGATAAG‑TTGGCGAGGGGTAAATAACCCGAACAGGCG‑GGGCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTTCGCAGCTATGTGGCGCGCTGGCTGAcc                                                                                                                 <  1:880039/239‑1 (MQ=255)                                                                                                                                       gattGCAACATGGGTTTC‑GGGTAAGCCGGCGACGGGGTAAATGCCCGAAACAGCGG‑AGACGTAACTTTTTGCGCCAGCCGGTGGAGGCGTGCCTCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAAGGTGCTGAGTATTCAGGGTACAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGGTGACCTAATTtgatg                                                                                                       <  1:595424/235‑1 (MQ=255)                                                                                                                                        gTTGCAACATGGGGTGGTGGATAAG‑CTGCGAAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTtgatga                                                                                                      <  2:757068/237‑1 (MQ=255)                                                                                                                                             aaacgggggTTCTGGATTAC‑GTGCGCAGGGGTAAATTTCGCCAACGGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAg                                                                                               <  1:1098790/234‑1 (MQ=255)                                                                                                                                              ccATGGGGGGCTGTATAAG‑CTGGGCAGGGGTAAATATCGCCAACAGGCG‑AGCCGTAATCTTTGGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAg                                                                                               <  1:594154/237‑1 (MQ=255)                                                                                                                                               cATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCggg                                                                                                                                                                      <  2:88524/166‑1 (MQ=255)                                                                                                                                                         tGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGGGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGGAAAGGGCGGGAGGGTAAATCGCAGGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCGGACGGAATTTGGTGATGTTCAGGGACAATGCAcg                                                                                   >  1:1269498/1‑238 (MQ=255)                                                                                                                                                         tGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCAAATTGCAc                                                                                    >  2:517646/1‑238 (MQ=255)                                                                                                                                                           gATCAG‑CTGCACGGTGTTAAATATCCGCAATCGGCG‑GGCTTTAAATTTTTGCGCCAGCCGGTTGAAGAGTTCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACcg                                                                                  <  1:1154449/238‑1 (MQ=255)                                                                                                                                                                g‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAata                                                                                                                                                 >  2:270843/1‑170 (MQ=255)                                                                                                                                                                g‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGTCGGGATGATAAATCGCCTGCCCATATTGTCCGCAGCCATGTGACGCGCTGGCTGACCGAATTTGATGAGGTTCAGGCATATTGCACCGCGCTa                                                                             >  2:1285202/1‑238 (MQ=255)                                                                                                                                                                          tGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTtga                                                                                                         <  2:870947/201‑1 (MQ=255)                                                                                                                                                                          tGGTAAATATCCGAAACAGGCG‑AGCCTGAATCTGTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTtgatg                                                                                                       <  2:1233550/203‑1 (MQ=255)                                                                                                                                                                            gTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGACTCATggc                                                                 >  2:44649/1‑239 (MQ=255)                                                                                                                                                                               atgtatccgaaatagtgctgCCTGTATTTTTTGCTCAAGCCGGTGTAAGAGTGCCACAAAATGGGGTACAGTGTTATTCAACAGGCTCTTGCGGATGGTTTGCG‑TAACGTGATTATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCACTATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAAGGTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATggcgg                                                               <  1:1035820/220‑1 (MQ=255)                                                                                                                                                                                  atCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGc                                                           <  2:913525/239‑1 (MQ=255)                                                                                                                                                                                    ccGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCCTGGCAAAGGGCGGGATGATCAATCGCATGCCAATATTGCCCGCAGCAATGTGTCGCGCTGGCTGACCGAATTTGATAATGTGCAGCCATATTGCACCACGCAAcggaacaaagccggaacggg                                                         >  1:344390/1‑220 (MQ=255)                                                                                                                                                                                    ccGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCGGCTATTTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGAACCGCGCTACCGCTCCATGGCGGCAACgg                                                         >  2:841565/1‑239 (MQ=255)                                                                                                                                                                                    ccGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCCCGCTACCGCATAATGGCGGCCGCgg                                                         >  2:1045581/1‑239 (MQ=255)                                                                                                                                                                                      ttttaagggcctgcttgatttttttGCGCCAGCCGGTGGAAAGGGGCCGGAAAATGGGTGTCAGTTTTATCCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTTCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCggg                                                        <  1:385547/218‑1 (MQ=255)                                                                                                                                                                                      gCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCgggg                                                       >  2:59949/1‑239 (MQ=255)                                                                                                                                                                                         aCAGGCG‑AGAATGAAACTTTTGAGCCAGCCGGTGGAAGAGTGCCGCAAAATGGAGTTCAATTTTATTCAACAAGCCCTGGCGGATGATTTGCC‑CAACGTGCTGATTATTCATGGTACAGGACGGGCTGATCAATCGCATGCCAATCTTATCCCTACGTCAGTGGCCCGATCGCGGACACAGTTTGATGATGTGCATGCCAATTGCCCGGCGCACCCGCCACAAGCCGGCAGCGGcgcgc                                                    >  1:137571/1‑238 (MQ=255)                                                                                                                                                                                               cgggCCTGAATTTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCGGGGCgtg                                                   <  2:1269498/232‑1 (MQ=255)                                                                                                                                                                                                   ccTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGGGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAACGGGGCGTGTTAAGTc                                            >  2:613985/1‑238 (MQ=255)                                                                                                                                                                                                   ccTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCGGGGCGGGGTACGTc                                            >  2:661920/1‑238 (MQ=255)                                                                                                                                                                                                            tttGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTTTCCGCAGCTATTTGGCGCGCTGGCTGACCGAATTTGAGGATGTTCAGGCATATTGCACCGCCCTCCGGCCTCATGGC‑GCAGAgggggggtggtacggccccctcccga                                 >  1:468632/1‑217 (MQ=255)                                                                                                                                                                                                              tGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCGGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGCCACCGGGGGGGGTTCCGTCggctctttatac                                >  2:385547/1‑228 (MQ=255)                                                                                                                                                                                                                              gtggTGCCTAAAAATGGGTTTCTGTTTTTTGCAAAAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTGTTCAGGGTAAAGGGCGGGATGATAAAGCGCATGCCAATATTGTCCGAAGCTATGGGGCGCGCTGGCTGACCGAATGTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGAAGCGGGGCGTGTTACGGCGCACTACGTAAAACGGCGCAGGCGAAg                 <  1:191376/235‑1 (MQ=255)                                                                                                                                                                                                                                   tccGGAAATTGGTTGTGAGTTTTATTCAAAAAGCCGTTGCGGATGGGTTGCT‑TAACGTGTTGATTATTCATGGTAAAGGGCGGGATGGTAAAGCGAATGCCAAGATTGTCCGCAGCTAGGTGGCGCGCTGGCTGACCGAATTTGATTATGTTCAGGCATATTGCCCCGCGCTACCGCAGCATGGCGGCAGCGGGGCGTGTTACGTCGAACTACGTAAAACGGCGCAGGCGAAGCAAGaa           <  1:324865/238‑1 (MQ=255)                                                                                                                                                                                                                                   gCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGGGCCCGGCGAAGCAAGaa           >  2:452057/1‑239 (MQ=255)                                                                                                                                                                                                                                         aaaTGGTTTTCAGGTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGGCGCAGGCGAAGCAAGAAAACTgg     <  2:468632/239‑1 (MQ=255)                                                                                                                                                                                                                                            tGGTGTTCAGTTTTATTCAACAATCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATGGCGCCGGGGTACCGCGTCATGGCGGCAGCGGGGGGTGTTACATCGACCTAAGTACAACTGCGAAGGGGAAGCGAGGAAAACGGGAg  >  1:1155658/1‑239 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             GATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCG‑AGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGGCGCAGGCGAAGCAAGAAAACTGGGAG  >  NC_000913/1406048‑1406518

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