Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,672	+T	coding (377/1701 nt)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,672	1	.	T	92.9%	180.1 / 4.7	58	coding (377/1701 nt)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base T (37/15);  total (38/20)
Fisher's exact test for biased strand distribution p-value = 9.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCT‑T‑C‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATC  >  NC_000913/69424‑69919                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                           gagaTTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑Gcttct‑ttc                                                                                                                                                                                                                                                           >  1:100664/1‑248 (MQ=255)     tttATGCCCGGCGCTGCACCGCCCGCGACGAATACCCGGCGGGCGGGTGGTCCCGAAAAGCAGAGCTGGCACCCAGTCGCACAGCCCAATCCACGATGCGGCAGATGGCGCCCCGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGAAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AAcc                                                                                                                                                                                                                                                      <  1:64586/250‑1 (MQ=255)       tATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTCCCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGcc                                                                                                                                                                                                                                                   >  1:223932/1‑250 (MQ=255)        aggCCCGGCGCTGCACCGCCCGCGACGAATACCCGGGGGGCGGGTGGCACCGAAAAGCAGAGCGGGCCCCCAGGCGCAAAGCCCAACCCAGGATGCGGCAGATGGCGCCCCGGCGGTGTCGGGGCGAGTCACAGGGAGAATTTTTCCCAAAAACCA‑TCCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACAGTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCa                                                                                                                                                                                                                                                   <  1:273203/248‑1 (MQ=255)        aggCCCGGCGCTGCAACGCCCGCGATAAATATCCGGGGGGCGGGTGGTCCCGGAAAGCAGAGCTGGCCCCCAGTCGCACAGCTCATCCCACGAGGCGGCAGATTGCGCCCCGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGACCCA‑TCCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCa                                                                                                                                                                                                                                                   <  1:422765/248‑1 (MQ=255)        aTGCCGGGCGCTGCAACGTCCGCGACGAATACCCTGGGGGCGGGTGGTACCGGAAAGCAGAGCTGGCCCCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCGCGGCGCGGTCCGGGCGAGTCACATGCAGGATTTTTGCCCAGACCCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAg                                                                                                                                                                                                                                                  <  1:671868/251‑1 (MQ=255)         cccccGGCGCTGCAACGCCCCCGACAAATCCCCTGCGGGCGGGTGGTCCCGAAAAGCAGAGCGGGCACCCGGTCCCACAGCCCAACCCACGATGCGGCAGATTGCGCCACGGCGCTGTCTTGGCGAGTCACATGCAGGATTTTTGCCCA‑AACCAATTCGCTGGAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAgt                                                                                                                                                                                                                                                 <  1:541647/249‑1 (MQ=255)             cGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCAGGCACCCAGTCGCACAGCTCCATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGGCAACATTGCCCGGCGCGTGGACCAAACGGGGAATCTCTTC‑GCTTCT‑TTC‑ACACGCAGTGTgg                                                                                                                                                                                                                                             >  1:487190/1‑251 (MQ=255)                          gccacgacgaATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCT‑T‑C‑AACCGAAGTGTGGTCTTTCCACAATa                                                                                                                                                                                                                                <  1:249900/247‑1 (MQ=255)                              cgacgaATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGCGTAGTCAACATTGACCGGCGCGTGGCACAAACGGGGAATCTCTTC‑GCTTCT‑TTC‑AACCGACGTGGGGTCTTTCCACAATCCGa                                                                                                                                                                                                                             >  2:158367/1‑250 (MQ=255)                                acgaATATCCTGCGGGGGGGTGGTCCCGAAAAGCAGAGCGGGCCCCCAGCCGCACAGCCCAACCCACGAGGCGGCAGATTGCCCCCCGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTCCCCAGACCCA‑TTCGCTGAAATAAATTCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑a                                                                                                                                                                                                                         <  1:347942/251‑1 (MQ=255)                                   aaTATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCACTATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑CACCGCAGTGTGGTCTTTCACCAATACgtaccacg                                                                                                                                                                                                                       >  2:498320/1‑244 (MQ=255)                                      atCCTGCGGGCTGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGAGGCACAACCGGGAAATCTCTTCCGCTTCT‑T‑C‑AACCGCAGGGAGGTCTTTCCACAATACGAAC‑ATcgcat                                                                                                                                                                                                                   >  2:482452/1‑249 (MQ=255)                                         cTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCg                                                                                                                                                                                                                >  1:410978/1‑251 (MQ=255)                                                   gtggtACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGc                                                                                                                                                                                                      >  1:154773/1‑251 (MQ=255)                                                          cGGAAAGGAAAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGAGGCAGATTGCGTCACGGCGCTGTCCTGGCGAGTCACATGCCGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATAACGCAAATATAGCGGGAGTAGTCACCTTTGACAGGCGCGCGGCAAAAACGGGTAAT‑TTATGCGCTTCT‑TTC‑AACCGCGGTGTGGTATTTCCACAATAC‑AACGCTCGCGTTCGGGTTTTCggacaacacc                                                                                                                                                                                               >  2:231178/1‑243 (MQ=255)                                                          cGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCGGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCCCAAACGGGGAATCTCTTC‑GCTTCT‑TTC‑AACCGCCGTGTGGTCTTTCCACAATACGAAC‑CTCGCGTTCGGGTTTTCGGCAAAcacc                                                                                                                                                                                               >  2:377896/1‑248 (MQ=255)                                                             aaggcaggcgggcccccattccgcccaGCTCAACCCACGTTGCGACAGTGTGCCCCACGGCGCTGTCCTGGAGAGTCACATGCAGTATTTTTCCCAATACCCA‑TCCGGGGAAATAAAATCCGCCAATATAGCGGGAGTAGTCAACATGGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCggg                                                                                                                                                                                            <  2:471775/226‑1 (MQ=255)                                                              caGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCCCAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCgggg                                                                                                                                                                                           >  2:13156/2‑250 (MQ=255)                                                               cgCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGAGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGgcg                                                                                                                                                                                          >  2:486958/2‑251 (MQ=255)                                                               aGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCCGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGgcg                                                                                                                                                                                          >  2:121711/1‑251 (MQ=255)                                                                   gagCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAgc                                                                                                                                                                                      >  1:185512/1‑251 (MQ=255)                                                                       tGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGcc                                                                                                                                                                                   >  1:526489/1‑250 (MQ=255)                                                                            cacAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑ATCGCTGGAATAAATACCGCCAATCTAGCGGGAGTAGTCAACATTGCCCGGCGCGAGGCACAAACGGGCAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGATGTCTTTCCACAATCCCAAC‑ATCGCGCTCGGGTTTTCGGTCAACGCAGGGCGCACCGcc                                                                                                                                                                                   >  2:157761/3‑245 (MQ=255)                                                                                   gCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGt                                                                                                                                                                      >  1:120786/1‑251 (MQ=255)                                                                                                   gATGCGGCAGATTGCGCCACGGCGCTGTCCTGTCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCACACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCCTCAATCGGTGc                                                                                                                                                      >  2:40920/1‑251 (MQ=255)                                                                                                       cGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCCTCAATCGGTCCGGGc                                                                                                                                                  >  2:334486/1‑251 (MQ=255)                                                                                                          cAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTc                                                                                                                                               >  1:169799/1‑251 (MQ=255)                                                                                                             aTTGCGCCCCGGCGCGGCCCTGGCGAGTCAAATGAAGGTTTTTGGCCCAAAACCA‑TTCGCTGAAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAg                                                                                                                                            <  1:11674/251‑1 (MQ=255)                                                                                                                    cACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCT‑T‑C‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGTCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTAc                                                                                                                                     <  1:482452/251‑1 (MQ=255)                                                                                                                      cGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTg                                                                                                                                   >  1:30256/1‑251 (MQ=255)                                                                                                                         cgcTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTc                                                                                                                                 >  1:82449/1‑250 (MQ=255)                                                                                                                                    cGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGAGCGGGCGTCGAGCCGGTACTGTCAACGCCAATccc                                                                                                                     >  1:55103/1‑251 (MQ=255)                                                                                                                                        tCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTc                                                                                                                                                                                                         >  1:707175/1‑163 (MQ=255)                                                                                                                                               cAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGCAGTCAACATTGCCCGGCGCGTGGCCCAAACGGGAAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGCGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGGTTTCGGCAAACTCCGGgc                                                                                                                                                                                           >  1:407827/1‑170 (MQ=255)                                                                                                                                                  gATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTgcgc                                                                                                      >  2:114230/1‑251 (MQ=255)                                                                                                                                                         gCCCAGGACCA‑TTCGCTGGAATCAATACCGACACTATAGCGGGAGTAGTCAACATTGCACGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCCTCT‑TTC‑CACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTATCGGCAAACTACGGGCGCAGCGCCAGCAAGTTACCGTCGGTTTCAATCGGTGAGGGCACCGAGCCGGGCCTGTCAACGCCAATCCCG‑ACCCCAGCTGCGCGCTGttt                                                                                               >  2:561342/1‑250 (MQ=255)                                                                                                                                                           ccaaacccA‑TTCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAAc                                                                                                                                                                                                  <  1:129280/145‑1 (MQ=255)                                                                                                                                                           ccAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAAc                                                                                                                                                                                                  >  2:129280/1‑151 (MQ=255)                                                                                                                                                              gAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGc                                                                                          >  1:591652/1‑251 (MQ=255)                                                                                                                                                                acccttccgcggaaaaaaataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCt                                                                                         <  1:377896/235‑1 (MQ=255)                                                                                                                                                                       cGCTGGAATAAATACCGCCAATATAGCGGGCGTAGTCAACATTGCCCGGCGCGCGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGCAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCCGCGTCAGCACGTTACCGTCGGCATCAATCGGTGCGGGTGTCGAGCAGGTACTGGCAACGCCAATCCCA‑ACCACGGCTGCGCGCTCTTCCACGCTAAActct                                                                                  >  2:240288/1‑251 (MQ=255)                                                                                                                                                                             ataatatcccgccaATGTAGCGGGAGTAGTCAACATTGCCCGTCGTGTGGCACAACCGGGTAATCTCTTC‑GCTTCT‑TGC‑AACCGCAGTGTGGTCTTTCCACAATCCTAAC‑ATCGCGTGCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACTCTAAGCTCTGCAAGc                                                                            <  2:254461/240‑1 (MQ=255)                                                                                                                                                                               taaataCCGCCCATATAGCGGGAGTAGTCAACATTGCCCAGCGCGTGGCCCAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCCCGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTCCTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGCCGCTAAGCTCTGccagacc                                                                          >  2:353139/1‑245 (MQ=255)                                                                                                                                                                                aaataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGt                                                                                                                                                                      >  2:235697/1‑159 (MQ=255)                                                                                                                                                                                 aaaacccccaATATAGCGGGAGTAGGCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCCCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAACCGGGGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACgg                                                                        <  2:185512/244‑1 (MQ=255)                                                                                                                                                                                                 cGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTc                                                        >  1:151853/1‑251 (MQ=255)                                                                                                                                                                                                            aacattGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCa                                              >  1:552361/5‑250 (MQ=255)                                                                                                                                                                                                                   cccGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTc                                      >  1:279913/1‑251 (MQ=255)                                                                                                                                                                                                                          cgTGGCACAAACGGGCAATCTATTC‑GCTGCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTGGGCAAACTCCGGGCGCAGCGACAGCACGTTACCGTCGGCATCAATCGCCGCGGGCGTCGGGCGGGTACTGTCAACGCCAATCCCGCACC‑CAGCTGCGCGCTGTTCGACGCTACGATCTGCAAGCACGGTTTTCAGAGCCGCTTACATTGCCTCATTGTAGTccggcgg                                >  2:176268/1‑244 (MQ=255)                                                                                                                                                                                                                           gTGGCACATACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCACACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGCGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCCCGGTTTTCAGTGCCGCTGCCATTGCCTCAATGTAGTTACGCGgat                              >  2:471955/1‑251 (MQ=255)                                                                                                                                                                                                                               cacaAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGatgatg                          <  2:526489/251‑1 (MQ=255)                                                                                                                                                                                                                                      gggTAATCTCT‑GCGCTCCTGT‑C‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACGCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAAc                   <  1:629945/251‑1 (MQ=255)                                                                                                                                                                                                                                       ggTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGAGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGTGCGTCGAGCCGGTACTGTCAACGCCAATCTCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGAATTATGACGGAACt                  >  1:167012/1‑251 (MQ=255)                                                                                                                                                                                                                                          aaTCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGt               <  1:114230/251‑1 (MQ=255)                                                                                                                                                                                                                                          aaTCTCTTC‑GCTTCT‑TGC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑CTCGCGTTCGGGTTTTCGGCAACCTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGt               <  1:56896/251‑1 (MQ=255)                                                                                                                                                                                                                                               tcttcGCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTAtt           >  1:386380/5‑251 (MQ=255)                                                                                                                                                                                                                                                cttcccttct‑t‑CCAACCGCAGTTTGGTCTTTCCACAATACGACC‑ACCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCCCCAGCCCGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGtt              <  2:410978/242‑1 (MQ=255)                                                                                                                                                                                                                                                  tcGCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCg         <  1:587025/249‑1 (MQ=255)                                                                                                                                                                                                                                                   cGCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCggg       <  1:498320/251‑1 (MQ=255)                                                                                                                                                                                                                                                        tt‑t‑c‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATc  >  1:400173/2‑250 (MQ=255)                                                                                                                                                                                                                                                        tt‑t‑c‑AACCGCAGGGTGGTCTTTCCAGAATACGAAC‑ATCGCGTTCGGGTTTTCGACATACTCCGGGCGCAGCGCCAGCACTTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTAGTCGGGGCATc  <  1:482143/249‑1 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                           GAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCT‑T‑C‑AACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATC  >  NC_000913/69424‑69919

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