Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,665	Δ1 bp	coding (384/1701 nt)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,665	0	C	.	90.9%	199.6 / 12.9	55	coding (384/1701 nt)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base . (35/15);  total (37/18)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

GAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCT‑T‑CAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGG  >  NC_000913/69424‑69914                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                              gagaTTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑Gcttct‑ttc                                                                                                                                                                                                                                                     >  1:100664/1‑248 (MQ=255)     tttATGCCCGGCGCTGCACCGCCCGCGACGAATACCCGGCGGGCGGGTGGTCCCGAAAAGCAGAGCTGGCACCCAGTCGCACAGCCCAATCCACGATGCGGCAGATGGCGCCCCGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGAAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAAcc                                                                                                                                                                                                                                                 <  1:64586/250‑1 (MQ=255)       tATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTCCCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGcc                                                                                                                                                                                                                                              >  1:223932/1‑250 (MQ=255)        aggCCCGGCGCTGCACCGCCCGCGACGAATACCCGGGGGGCGGGTGGCACCGAAAAGCAGAGCGGGCCCCCAGGCGCAAAGCCCAACCCAGGATGCGGCAGATGGCGCCCCGGCGGTGTCGGGGCGAGTCACAGGGAGAATTTTTCCCAAAAACCA‑TCCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACAGTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCa                                                                                                                                                                                                                                              <  1:273203/248‑1 (MQ=255)        aggCCCGGCGCTGCAACGCCCGCGATAAATATCCGGGGGGCGGGTGGTCCCGGAAAGCAGAGCTGGCCCCCAGTCGCACAGCTCATCCCACGAGGCGGCAGATTGCGCCCCGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGACCCA‑TCCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCa                                                                                                                                                                                                                                              <  1:422765/248‑1 (MQ=255)        aTGCCGGGCGCTGCAACGTCCGCGACGAATACCCTGGGGGCGGGTGGTACCGGAAAGCAGAGCTGGCCCCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCGCGGCGCGGTCCGGGCGAGTCACATGCAGGATTTTTGCCCAGACCCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAg                                                                                                                                                                                                                                             <  1:671868/251‑1 (MQ=255)         cccccGGCGCTGCAACGCCCCCGACAAATCCCCTGCGGGCGGGTGGTCCCGAAAAGCAGAGCGGGCACCCGGTCCCACAGCCCAACCCACGATGCGGCAGATTGCGCCACGGCGCTGTCTTGGCGAGTCACATGCAGGATTTTTGCCCA‑AACCAATTCGCTGGAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAgt                                                                                                                                                                                                                                            <  1:541647/249‑1 (MQ=255)             cGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCAGGCACCCAGTCGCACAGCTCCATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGGCAACATTGCCCGGCGCGTGGACCAAACGGGGAATCTCTTC‑GCTTCT‑TTCACACGCAGTGTgg                                                                                                                                                                                                                                        >  1:487190/1‑251 (MQ=255)                          gccacgacgaATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCT‑T‑CAACCGAAGTGTGGTCTTTCCACAATa                                                                                                                                                                                                                           <  1:249900/247‑1 (MQ=255)                              cgacgaATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGCGTAGTCAACATTGACCGGCGCGTGGCACAAACGGGGAATCTCTTC‑GCTTCT‑TTCAACCGACGTGGGGTCTTTCCACAATCCGa                                                                                                                                                                                                                        >  2:158367/1‑250 (MQ=255)                                acgaATATCCTGCGGGGGGGTGGTCCCGAAAAGCAGAGCGGGCCCCCAGCCGCACAGCCCAACCCACGAGGCGGCAGATTGCCCCCCGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTCCCCAGACCCA‑TTCGCTGAAATAAATTCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑a                                                                                                                                                                                                                    <  1:347942/251‑1 (MQ=255)                                   aaTATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCACTATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCCACCGCAGTGTGGTCTTTCACCAATACgtaccacg                                                                                                                                                                                                                  >  2:498320/1‑244 (MQ=255)                                      atCCTGCGGGCTGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGAGGCACAACCGGGAAATCTCTTCCGCTTCT‑T‑CAACCGCAGGGAGGTCTTTCCACAATACGAAC‑ATcgcat                                                                                                                                                                                                              >  2:482452/1‑249 (MQ=255)                                         cTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCg                                                                                                                                                                                                           >  1:410978/1‑251 (MQ=255)                                                   gtggtACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGc                                                                                                                                                                                                 >  1:154773/1‑251 (MQ=255)                                                          cGGAAAGGAAAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGAGGCAGATTGCGTCACGGCGCTGTCCTGGCGAGTCACATGCCGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATAACGCAAATATAGCGGGAGTAGTCACCTTTGACAGGCGCGCGGCAAAAACGGGTAAT‑TTATGCGCTTCT‑TTCAACCGCGGTGTGGTATTTCCACAATAC‑AACGCTCGCGTTCGGGTTTTCggacaacacc                                                                                                                                                                                          >  2:231178/1‑243 (MQ=255)                                                          cGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCGGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCCCAAACGGGGAATCTCTTC‑GCTTCT‑TTCAACCGCCGTGTGGTCTTTCCACAATACGAAC‑CTCGCGTTCGGGTTTTCGGCAAAcacc                                                                                                                                                                                          >  2:377896/1‑248 (MQ=255)                                                             aaggcaggcgggcccccattccgcccaGCTCAACCCACGTTGCGACAGTGTGCCCCACGGCGCTGTCCTGGAGAGTCACATGCAGTATTTTTCCCAATACCCA‑TCCGGGGAAATAAAATCCGCCAATATAGCGGGAGTAGTCAACATGGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCggg                                                                                                                                                                                       <  2:471775/226‑1 (MQ=255)                                                              caGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCCCAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCgggg                                                                                                                                                                                      >  2:13156/2‑250 (MQ=255)                                                               cgCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGAGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGgcg                                                                                                                                                                                     >  2:486958/2‑251 (MQ=255)                                                               aGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCCGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGgcg                                                                                                                                                                                     >  2:121711/1‑251 (MQ=255)                                                                   gagCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAgc                                                                                                                                                                                 >  1:185512/1‑251 (MQ=255)                                                                       tGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGcc                                                                                                                                                                              >  1:526489/1‑250 (MQ=255)                                                                            cacAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑ATCGCTGGAATAAATACCGCCAATCTAGCGGGAGTAGTCAACATTGCCCGGCGCGAGGCACAAACGGGCAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGATGTCTTTCCACAATCCCAAC‑ATCGCGCTCGGGTTTTCGGTCAACGCAGGGCGCACCGcc                                                                                                                                                                              >  2:157761/3‑245 (MQ=255)                                                                                   gCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGt                                                                                                                                                                 >  1:120786/1‑251 (MQ=255)                                                                                                   gATGCGGCAGATTGCGCCACGGCGCTGTCCTGTCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCACACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCCTCAATCGGTGc                                                                                                                                                 >  2:40920/1‑251 (MQ=255)                                                                                                       cGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCCTCAATCGGTCCGGGc                                                                                                                                             >  2:334486/1‑251 (MQ=255)                                                                                                          cAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTc                                                                                                                                          >  1:169799/1‑251 (MQ=255)                                                                                                             aTTGCGCCCCGGCGCGGCCCTGGCGAGTCAAATGAAGGTTTTTGGCCCAAAACCA‑TTCGCTGAAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAg                                                                                                                                       <  1:11674/251‑1 (MQ=255)                                                                                                                    cACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCT‑T‑CAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGTCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTAc                                                                                                                                <  1:482452/251‑1 (MQ=255)                                                                                                                      cGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTg                                                                                                                              >  1:30256/1‑251 (MQ=255)                                                                                                                         cgcTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTc                                                                                                                            >  1:82449/1‑250 (MQ=255)                                                                                                                                    cGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGAGCGGGCGTCGAGCCGGTACTGTCAACGCCAATccc                                                                                                                >  1:55103/1‑251 (MQ=255)                                                                                                                                        tCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTc                                                                                                                                                                                                    >  1:707175/1‑163 (MQ=255)                                                                                                                                               cAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGCAGTCAACATTGCCCGGCGCGTGGCCCAAACGGGAAATCTCTTC‑GCTTCT‑TTCAACCGCAGCGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGGTTTCGGCAAACTCCGGgc                                                                                                                                                                                      >  1:407827/1‑170 (MQ=255)                                                                                                                                                  gATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTgcgc                                                                                                 >  2:114230/1‑251 (MQ=255)                                                                                                                                                         gCCCAGGACCA‑TTCGCTGGAATCAATACCGACACTATAGCGGGAGTAGTCAACATTGCACGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCCTCT‑TTCCACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTATCGGCAAACTACGGGCGCAGCGCCAGCAAGTTACCGTCGGTTTCAATCGGTGAGGGCACCGAGCCGGGCCTGTCAACGCCAATCCCG‑ACCCCAGCTGCGCGCTGttt                                                                                          >  2:561342/1‑250 (MQ=255)                                                                                                                                                           ccaaacccA‑TTCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAAc                                                                                                                                                                                             <  1:129280/145‑1 (MQ=255)                                                                                                                                                           ccAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAAc                                                                                                                                                                                             >  2:129280/1‑151 (MQ=255)                                                                                                                                                              gAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGc                                                                                     >  1:591652/1‑251 (MQ=255)                                                                                                                                                                acccttccgcggaaaaaaataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCt                                                                                    <  1:377896/235‑1 (MQ=255)                                                                                                                                                                       cGCTGGAATAAATACCGCCAATATAGCGGGCGTAGTCAACATTGCCCGGCGCGCGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGCAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCCGCGTCAGCACGTTACCGTCGGCATCAATCGGTGCGGGTGTCGAGCAGGTACTGGCAACGCCAATCCCA‑ACCACGGCTGCGCGCTCTTCCACGCTAAActct                                                                             >  2:240288/1‑251 (MQ=255)                                                                                                                                                                             ataatatcccgccaATGTAGCGGGAGTAGTCAACATTGCCCGTCGTGTGGCACAACCGGGTAATCTCTTC‑GCTTCT‑TGCAACCGCAGTGTGGTCTTTCCACAATCCTAAC‑ATCGCGTGCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACTCTAAGCTCTGCAAGc                                                                       <  2:254461/240‑1 (MQ=255)                                                                                                                                                                               taaataCCGCCCATATAGCGGGAGTAGTCAACATTGCCCAGCGCGTGGCCCAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCCCGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTCCTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGCCGCTAAGCTCTGccagacc                                                                     >  2:353139/1‑245 (MQ=255)                                                                                                                                                                                aaataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGt                                                                                                                                                                 >  2:235697/1‑159 (MQ=255)                                                                                                                                                                                 aaaacccccaATATAGCGGGAGTAGGCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCCCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAACCGGGGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACgg                                                                   <  2:185512/244‑1 (MQ=255)                                                                                                                                                                                                 cGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTc                                                   >  1:151853/1‑251 (MQ=255)                                                                                                                                                                                                            aacattGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCa                                         >  1:552361/5‑250 (MQ=255)                                                                                                                                                                                                                   cccGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTc                                 >  1:279913/1‑251 (MQ=255)                                                                                                                                                                                                                          cgTGGCACAAACGGGCAATCTATTC‑GCTGCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTGGGCAAACTCCGGGCGCAGCGACAGCACGTTACCGTCGGCATCAATCGCCGCGGGCGTCGGGCGGGTACTGTCAACGCCAATCCCGCACC‑CAGCTGCGCGCTGTTCGACGCTACGATCTGCAAGCACGGTTTTCAGAGCCGCTTACATTGCCTCATTGTAGTccggcgg                           >  2:176268/1‑244 (MQ=255)                                                                                                                                                                                                                           gTGGCACATACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCACACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGCGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCCCGGTTTTCAGTGCCGCTGCCATTGCCTCAATGTAGTTACGCGgat                         >  2:471955/1‑251 (MQ=255)                                                                                                                                                                                                                               cacaAACGGGTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGatgatg                     <  2:526489/251‑1 (MQ=255)                                                                                                                                                                                                                                      gggTAATCTCT‑GCGCTCCTGT‑CAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACGCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAAc              <  1:629945/251‑1 (MQ=255)                                                                                                                                                                                                                                       ggTAATCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGAGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGTGCGTCGAGCCGGTACTGTCAACGCCAATCTCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGAATTATGACGGAACt             >  1:167012/1‑251 (MQ=255)                                                                                                                                                                                                                                          aaTCTCTTC‑GCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGt          <  1:114230/251‑1 (MQ=255)                                                                                                                                                                                                                                          aaTCTCTTC‑GCTTCT‑TGCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑CTCGCGTTCGGGTTTTCGGCAACCTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGt          <  1:56896/251‑1 (MQ=255)                                                                                                                                                                                                                                               tcttcGCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTAtt      >  1:386380/5‑251 (MQ=255)                                                                                                                                                                                                                                                  tcGCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCg    <  1:587025/249‑1 (MQ=255)                                                                                                                                                                                                                                                   cGCTTCT‑TTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCggg  <  1:498320/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                              GAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCT‑T‑CAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGG  >  NC_000913/69424‑69914

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