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%	116.6 / 6.6	56	V102G (GTG→GGG)	smrA	DNA endonuclease
Reads supporting (aligned to +/- strand):  ref base T (15/26);  new base G (12/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGGCGCAGGCGAAGCAAGAAAACTGGGAG  >  NC_000913/1406048‑1406518                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             gATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGtt                                                                                                                                                                                                                                           <  1:848839/239‑1 (MQ=255)             cATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGTAAGGGGTGCAACAATGGGTGCCGCGAAAG‑CCGCGCCGTGGGAAATATCCGCAACAGGCGGGCCCGAATCTTTTGCCCCCACCGGTAGAAAATGGACCAACAAAGGGGTTGAGTTTTAttga                                                                                                                                                                                                                               >  2:436456/1‑237 (MQ=255)              aTCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGGGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGGGCCGCAAAATGGGGGTACGTTTTATTca                                                                                                                                                                                                                               >  1:595424/1‑238 (MQ=255)               ttccaCGCGTAACCAACGTTCGCCGGAGCGTTTCAACCCGCGGCAGGTTGAAAGTTTCCTCCCCCGCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTcaac                                                                                                                                                                                                                             <  2:442730/235‑1 (MQ=255)                    cgcgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGGCTGAATCTTTTGCGCCAGCCGGGGGAAGAGGGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGc                                                                                                                                                                                                                         >  1:575012/1‑238 (MQ=255)                    actgTAACCACCGTTCCCCGCACGGTTCCGCCCCGCTGCAGTTTGACAATTCCCTCACCACCGGGTTTCTCG‑AAATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGcc                                                                                                                                                                                                                        <  2:155512/236‑1 (MQ=255)                           ccAACGTGCGCCGCAGCGTATCGACCCGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCgg                                                                                                                                                                                                                 <  1:496579/239‑1 (MQ=255)                             aaCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGAt                                                                                                                                                                                                               <  1:848698/239‑1 (MQ=255)                                gTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAAAATGGGGTGCTGGATAAA‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGAgggg                                                                                                                                                                                                            >  1:1154449/1‑238 (MQ=255)                                   cgccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGGGTGCCGCAAAATGGGGTTCAGTTTTATTTAACAAGCCCTGGCGGAGGGtttt                                                                                                                                                                                                         >  1:1098790/1‑238 (MQ=255)                                   cgccgcAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGGAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGGGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTTAACAAGCCCTGGCGGATggcgg                                                                                                                                                                                                          >  1:155512/1‑235 (MQ=255)                                              atctccacGCTGCATGTTGATAATTTCCTCACCGCCGGATTTCGCG‑ACATCCTCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCt                                                                                                                                                                                             <  1:1016735/234‑1 (MQ=255)                                                             cTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGGTTGCG‑TAACGTGCTGATTATTCATGGTaa                                                                                                                                                                              >  1:880039/1‑239 (MQ=255)                                                                   aaTTTCCTCTCCCCCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCg                                                                                                                                                                        <  1:432322/239‑1 (MQ=255)                                                                         cTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGcc                                                                                                                                                                                                                        <  1:151472/186‑1 (MQ=255)                                                                              caccGGATTTGTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCa                                                                                                                                                                                    <  1:725800/216‑1 (MQ=255)                                                                                ccGTATTTCTCGTACA‑CATCCTACTTAGTCAGCGCCTGTGGTTTGGGCGGGAATGGTTATAGAAGCGGGGGGTGGATAAA‑CGAAGAAGGGAAAAACAGCCGCAAGAGACGGTCCTGACTCTTTTGCGCCAGCCGGTGGAAGAAAGCCGCAACATGGGGTTCAGTTTTATTCAACAAGCCATGGCGGATGGTGTGCG‑TAACTTGCTGATTATTCATGGTAAAGGGCGGGATGATaaagct                                                                                                                                                           <  2:1013431/239‑4 (MQ=255)                                                                                ccGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCCAAATGGGGTTCAGGTTTATTCAACAAGCCCTGGCGGATGGTTTGGG‑TAACGTGCTGGTTTTTCATGGTAAAGGGCGGGATGGTAAATc                                                                                                                                                            >  1:896413/1‑238 (MQ=255)                                                                                ccGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGGGGAAGAGAGCCGCAAAATGGGGTTCCGTTTTATTCCACAAGCCCTGGCGGATGGTTTTCG‑TAACGTTCTGATTTTTTATGGTAAAGGGCGGGGTGATAAAtg                                                                                                                                                            >  1:1070824/1‑237 (MQ=255)                                                                                ccGGAGTCCTCG‑AATTCACCCCACAACGGCAGCCTGTGGGGTTTCGGGGGGAAGGGTTGAAACATGGGGTGCTGGATAAG‑CTGCGCAGGGGTAAATATCCGCAACAGGCGAGGCGGAATCTTTTGCGCCAGCAGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCg                                                                                                                                                           <  2:809305/239‑1 (MQ=255)                                                                                acGGATTTCCCG‑GCAACATCCCACTAAGGCAGGCGCTTGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCg                                                                                                                                                           <  1:127949/238‑1 (MQ=255)                                                                                        ctcG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGATGCCACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGGAACAGGCGTGCCCGACTCTTTTGCGCCGGCCGGGGGTAGAGTGCCCCAAACTGGGGTTCAGCTTTATTCTAAAAGCCCTGGCGGCTAGTTTGCGCAAAGGTGCT‑ATTATTCATTGTAAAAAGCGGGATGGTAAAACgacagccaa                                                                                                                                                   >  2:921802/1‑231 (MQ=255)                                                                                               atcatcCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGGGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCCGGCGGATGGGTTTCG‑TAACGTGCTGATTTTTCCTGGTAAAGGGCGGGGTGATAAATCGCATGCCCATATTg                                                                                                                                              >  1:594154/1‑238 (MQ=255)                                                                                                tcatcCCAGGAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTTCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTc                                                                                                                                            <  1:93180/239‑1 (MQ=255)                                                                                                  atcCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTTCGCCAGCCGGTGGAAGAGAGCCGCAAAATGGGGTTCAGTTTTATTCAAAAAGCCCTGCCGGAGGGTTTGTG‑GAACGTGCTGATTATTTCTTGTAAAAGGCGGGATGATAAATGGCATGCCAATAtt                                                                                                                                               >  1:661247/1‑234 (MQ=255)                                                                                                        cTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGGAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGGTTGCG‑GAACGTGCTGGTTATTCATGGTAAAGGTCGGGATGGTAAATGGCATGCCAAAATTTTccccgacta                                                                                                                                    >  1:410749/1‑232 (MQ=255)                                                                                                                 ccGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGGGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGGTTATTCATGGTAAAGGGCGGGATGATAAATc                                                                                                                                                            >  1:984068/1‑206 (MQ=255)                                                                                                                         gtgtCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCgctggctg                                                                                                                    <  1:921802/235‑1 (MQ=255)                                                                                                                           ttCGGGGGGAAGGGTTGCAACAGGGGGGGCTGGATAAG‑TTGGCGAGGGGTAAATAACCCGAACAGGCGGGGCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTTCGCAGCTATGTGGCGCGCTGGCTGAcc                                                                                                                 <  2:880039/239‑1 (MQ=255)                                                                                                                                       gattGCAACATGGGTTTC‑GGGTAAGCCGGCGACGGGGTAAATGCCCGAAACAGCGGAGACGTAACTTTTTGCGCCAGCCGGTGGAGGCGTGCCTCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAAGGTGCTGAGTATTCAGGGTACAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGGTGACCTAATTtgatg                                                                                                       <  2:595424/235‑1 (MQ=255)                                                                                                                                        gTTGCAACATGGGGTGGTGGATAAG‑CTGCGAAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTtgatga                                                                                                      <  1:757068/237‑1 (MQ=255)                                                                                                                                             aaacgggggTTCTGGATTAC‑GTGCGCAGGGGTAAATTTCGCCAACGGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAg                                                                                               <  2:1098790/234‑1 (MQ=255)                                                                                                                                              ccATGGGGGGCTGTATAAG‑CTGGGCAGGGGTAAATATCGCCAACAGGCGAGCCGTAATCTTTGGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAg                                                                                               <  2:594154/237‑1 (MQ=255)                                                                                                                                               cATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCggg                                                                                                                                                                      <  1:88524/166‑1 (MQ=255)                                                                                                                                                         tGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGGGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGGAAAGGGCGGGAGGGTAAATCGCAGGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCGGACGGAATTTGGTGATGTTCAGGGACAATGCAcg                                                                                   >  2:1269498/1‑238 (MQ=255)                                                                                                                                                         tGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCAAATTGCAc                                                                                    >  1:517646/1‑238 (MQ=255)                                                                                                                                                           gATCAG‑CTGCACGGTGTTAAATATCCGCAATCGGCGGGCTTTAAATTTTTGCGCCAGCCGGTTGAAGAGTTCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACcg                                                                                  <  2:1154449/238‑1 (MQ=255)                                                                                                                                                                g‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAata                                                                                                                                                 >  1:270843/1‑170 (MQ=255)                                                                                                                                                                g‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGGGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGTCGGGATGATAAATCGCCTGCCCATATTGTCCGCAGCCATGTGACGCGCTGGCTGACCGAATTTGATGAGGTTCAGGCATATTGCACCGCGCTa                                                                             >  1:1285202/1‑238 (MQ=255)                                                                                                                                                                          tGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTtga                                                                                                         <  1:870947/201‑1 (MQ=255)                                                                                                                                                                          tGGTAAATATCCGAAACAGGCGAGCCTGAATCTGTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTtgatg                                                                                                       <  1:1233550/203‑1 (MQ=255)                                                                                                                                                                            gTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGACTCATggc                                                                 >  1:44649/1‑239 (MQ=255)                                                                                                                                                                                  atCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGc                                                           <  1:913525/239‑1 (MQ=255)                                                                                                                                                                                    ccGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCGGCTATTTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGAACCGCGCTACCGCTCCATGGCGGCAACgg                                                         >  1:841565/1‑239 (MQ=255)                                                                                                                                                                                    ccGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCCCGCTACCGCATAATGGCGGCCGCgg                                                         >  1:1045581/1‑239 (MQ=255)                                                                                                                                                                                      gCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCgggg                                                       >  1:59949/1‑239 (MQ=255)                                                                                                                                                                                         aCAGGCGAGAATGAAACTTTTGAGCCAGCCGGTGGAAGAGTGCCGCAAAATGGAGTTCAATTTTATTCAACAAGCCCTGGCGGATGATTTGCC‑CAACGTGCTGATTATTCATGGTACAGGACGGGCTGATCAATCGCATGCCAATCTTATCCCTACGTCAGTGGCCCGATCGCGGACACAGTTTGATGATGTGCATGCCAATTGCCCGGCGCACCCGCCACAAGCCGGCAGCGGcgcgc                                                    >  2:137571/1‑238 (MQ=255)                                                                                                                                                                                              cgggCCTGAATTTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCGGGGCgtg                                                   <  1:1269498/232‑1 (MQ=255)                                                                                                                                                                                                  ccTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGGGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAACGGGGCGTGTTAAGTc                                            >  1:613985/1‑238 (MQ=255)                                                                                                                                                                                                  ccTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCGGGGCGGGGTACGTc                                            >  1:661920/1‑238 (MQ=255)                                                                                                                                                                                                             tGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCGGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGCCACCGGGGGGGGTTCCGTCggctctttatac                                >  1:385547/1‑228 (MQ=255)                                                                                                                                                                                                                             gtggTGCCTAAAAATGGGTTTCTGTTTTTTGCAAAAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTGTTCAGGGTAAAGGGCGGGATGATAAAGCGCATGCCAATATTGTCCGAAGCTATGGGGCGCGCTGGCTGACCGAATGTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGAAGCGGGGCGTGTTACGGCGCACTACGTAAAACGGCGCAGGCGAAg                 <  2:191376/235‑1 (MQ=255)                                                                                                                                                                                                                                  tccGGAAATTGGTTGTGAGTTTTATTCAAAAAGCCGTTGCGGATGGGTTGCT‑TAACGTGTTGATTATTCATGGTAAAGGGCGGGATGGTAAAGCGAATGCCAAGATTGTCCGCAGCTAGGTGGCGCGCTGGCTGACCGAATTTGATTATGTTCAGGCATATTGCCCCGCGCTACCGCAGCATGGCGGCAGCGGGGCGTGTTACGTCGAACTACGTAAAACGGCGCAGGCGAAGCAAGaa           <  2:324865/238‑1 (MQ=255)                                                                                                                                                                                                                                  gCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGGGCCCGGCGAAGCAAGaa           >  1:452057/1‑239 (MQ=255)                                                                                                                                                                                                                                        aaaTGGTTTTCAGGTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGGCGCAGGCGAAGCAAGAAAACTgg     <  1:468632/239‑1 (MQ=255)                                                                                                                                                                                                                                           tGGTGTTCAGTTTTATTCAACAATCCCTGGCGGATGGTTTGCG‑TAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATGGCGCCGGGGTACCGCGTCATGGCGGCAGCGGGGGGTGTTACATCGACCTAAGTACAACTGCGAAGGGGAAGCGAGGAAAACGGGAg  >  2:1155658/1‑239 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             GATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCG‑ACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAG‑CTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCG‑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 < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑