Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,605	A→G	G148G (GGT→GGC)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,605	0	A	G	94.1%	149.7 / ‑2.1	51	G148G (GGT→GGC)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (35/13);  total (36/15)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AGAAACTGGCTGGCGGCAGGCCGCCCCAGCTTTCGTGCCACAG‑AGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTT‑CAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCAC  >  NC_000913/69382‑69847                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     aGAAACTGGCGGGCGGCAAGCCCCCCCAGCTTTCG‑GTCCCAGTAGATTGTGGCGCGGCCCGGCAACGTCCGCGACGAATATCCGGGGGGCGGGTGGTACCGGAAAGCAGGGCTGGCCCCCATCCGCACAGCCCAATCCACGAGGCGGCATATTGCGCCAGGGCGCTGTCCGGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGTTGGAAAAAAAACCGCCAATATAGCGGGAGTAGTCAACATTg                                                                                                                                                                                                                           <  2:590467/251‑1 (MQ=255)   aaaCTGGCTAGCGGCAAGCCGCCCCAGCTTTCGTGCCTCAG‑AGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCCCGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATATATCCCGCCAATCTAGCGGGAGGAGTTAACATTGcc                                                                                                                                                                                                                         >  2:304659/1‑251 (MQ=255)                    cccgccccgcctttgttcccACAG‑AGTTTTAGCCCCGGCGCTGCACCGCCCGGGACGATTCCCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCGGGCACCCGGCCCCCCAGCTCAATCCCCGATGCGCCATATTGGCCCGCGGCGCTGCCCTGGCGAGCCACATGCAGGATTTTTGCCCAGACCCA‑TTCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAAc                                                                                                                                                                                                        <  1:267705/233‑1 (MQ=255)                            aGCTTTCGTGCCACAG‑AGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAAtc                                                                                                                                                                                                >  1:206303/1‑251 (MQ=255)                                ttCGTGCCACAG‑CGATTTATGCCCGGCGCTGCAACGTCCGCGACGACTATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCCGGATTTTTGCCCAGAACCA‑TTCGCTGGAATACATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAACCGGGACATCTCtt                                                                                                                                                                                            >  2:273203/1‑251 (MQ=255)                                ttCGTCCCAAAG‑AGATTTACGCCCGGCGCCGCAACGCCCGCGACAAATACCCTGGGGGCGGGGGGTCCCGAAAAGCAGAGCTGCCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCCCGGCGCTTTCCTGGGGAGTCACATGCAGGTTTTTTGCCCAGAACCA‑TTCGCTGGAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAAtctc                                                                                                                                                                                              <  1:614378/249‑1 (MQ=255)                                           g‑agaTTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑Gcttctttc                                                                                                                                                                                 >  1:100664/1‑248 (MQ=255)                                                tttATGCCCGGCGCTGCACCGCCCGCGACGAATACCCGGCGGGCGGGTGGTCCCGAAAAGCAGAGCTGGCACCCAGTCGCACAGCCCAATCCACGATGCGGCAGATGGCGCCCCGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGAAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAAcc                                                                                                                                                                             <  1:64586/250‑1 (MQ=255)                                                  tATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTCCCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGcc                                                                                                                                                                          >  1:223932/1‑250 (MQ=255)                                                   aggCCCGGCGCTGCACCGCCCGCGACGAATACCCGGGGGGCGGGTGGCACCGAAAAGCAGAGCGGGCCCCCAGGCGCAAAGCCCAACCCAGGATGCGGCAGATGGCGCCCCGGCGGTGTCGGGGCGAGTCACAGGGAGAATTTTTCCCAAAAACCA‑TCCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACAGTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCa                                                                                                                                                                          <  1:273203/248‑1 (MQ=255)                                                   aggCCCGGCGCTGCAACGCCCGCGATAAATATCCGGGGGGCGGGTGGTCCCGGAAAGCAGAGCTGGCCCCCAGTCGCACAGCTCATCCCACGAGGCGGCAGATTGCGCCCCGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGACCCA‑TCCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCa                                                                                                                                                                          <  1:422765/248‑1 (MQ=255)                                                   aTGCCGGGCGCTGCAACGTCCGCGACGAATACCCTGGGGGCGGGTGGTACCGGAAAGCAGAGCTGGCCCCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCGCGGCGCGGTCCGGGCGAGTCACATGCAGGATTTTTGCCCAGACCCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAg                                                                                                                                                                         <  1:671868/251‑1 (MQ=255)                                                    cccccGGCGCTGCAACGCCCCCGACAAATCCCCTGCGGGCGGGTGGTCCCGAAAAGCAGAGCGGGCACCCGGTCCCACAGCCCAACCCACGATGCGGCAGATTGCGCCACGGCGCTGTCTTGGCGAGTCACATGCAGGATTTTTGCCCA‑AACCAATTCGCTGGAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAgt                                                                                                                                                                        <  1:541647/249‑1 (MQ=255)                                                        cGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCAGGCACCCAGTCGCACAGCTCCATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGGCAACATTGCCCGGCGCGTGGACCAAACGGGGAATCTCTTC‑GCTTCTTTCACACGCAGTGTgg                                                                                                                                                                    >  1:487190/1‑251 (MQ=255)                                                                     gccacgacgaATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTT‑CAACCGAAGTGTGGTCTTTCCACAATa                                                                                                                                                       <  1:249900/247‑1 (MQ=255)                                                                         cgacgaATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGCGTAGTCAACATTGACCGGCGCGTGGCACAAACGGGGAATCTCTTC‑GCTTCTTTCAACCGACGTGGGGTCTTTCCACAATCCGa                                                                                                                                                    >  2:158367/1‑250 (MQ=255)                                                                           acgaATATCCTGCGGGGGGGTGGTCCCGAAAAGCAGAGCGGGCCCCCAGCCGCACAGCCCAACCCACGAGGCGGCAGATTGCCCCCCGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTCCCCAGACCCA‑TTCGCTGAAATAAATTCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑a                                                                                                                                                <  1:347942/251‑1 (MQ=255)                                                                              aaTATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCACTATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCCACCGCAGTGTGGTCTTTCACCAATACgtaccacg                                                                                                                                              >  2:498320/1‑244 (MQ=255)                                                                                 atCCTGCGGGCTGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGAGGCACAACCGGGAAATCTCTTCCGCTTCTT‑CAACCGCAGGGAGGTCTTTCCACAATACGAAC‑ATcgcat                                                                                                                                          >  2:482452/1‑249 (MQ=255)                                                                                    cTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCg                                                                                                                                       >  1:410978/1‑251 (MQ=255)                                                                                              gtggtACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGc                                                                                                                             >  1:154773/1‑251 (MQ=255)                                                                                                     cGGAAAGGAAAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGAGGCAGATTGCGTCACGGCGCTGTCCTGGCGAGTCACATGCCGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATAACGCAAATATAGCGGGAGTAGTCACCTTTGACAGGCGCGCGGCAAAAACGGGTAAT‑TTATGCGCTTCTTTCAACCGCGGTGTGGTATTTCCACAATAC‑AACGCTCGCGTTCGGGTTTTCggacaacacc                                                                                                                      >  2:231178/1‑243 (MQ=255)                                                                                                     cGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCGGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCCCAAACGGGGAATCTCTTC‑GCTTCTTTCAACCGCCGTGTGGTCTTTCCACAATACGAAC‑CTCGCGTTCGGGTTTTCGGCAAAcacc                                                                                                                      >  2:377896/1‑248 (MQ=255)                                                                                                     cGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGgcg                                                                                                                                                                                                                    >  1:557901/1‑159 (MQ=255)                                                                                                        aaggcaggcgggcccccattccgcccaGCTCAACCCACGTTGCGACAGTGTGCCCCACGGCGCTGTCCTGGAGAGTCACATGCAGTATTTTTCCCAATACCCA‑TCCGGGGAAATAAAATCCGCCAATATAGCGGGAGTAGTCAACATGGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCggg                                                                                                                   <  2:471775/226‑1 (MQ=255)                                                                                                         caGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCCCAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCgggg                                                                                                                  >  2:13156/2‑250 (MQ=255)                                                                                                          cgCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGAGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGgcg                                                                                                                 >  2:486958/2‑251 (MQ=255)                                                                                                          aGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCCGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGgcg                                                                                                                 >  2:121711/1‑251 (MQ=255)                                                                                                              gagCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAgc                                                                                                             >  1:185512/1‑251 (MQ=255)                                                                                                                  tGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGcc                                                                                                          >  1:526489/1‑250 (MQ=255)                                                                                                                       cacAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑ATCGCTGGAATAAATACCGCCAATCTAGCGGGAGTAGTCAACATTGCCCGGCGCGAGGCACAAACGGGCAATCTCTTC‑GCTTCTTTCAACCGCAGTGATGTCTTTCCACAATCCCAAC‑ATCGCGCTCGGGTTTTCGGTCAACGCAGGGCGCACCGcc                                                                                                          >  2:157761/3‑245 (MQ=255)                                                                                                                              gCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGt                                                                                             >  1:120786/1‑251 (MQ=255)                                                                                                                                              gATGCGGCAGATTGCGCCACGGCGCTGTCCTGTCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCACACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCCTCAATCGGTGc                                                                             >  2:40920/1‑251 (MQ=255)                                                                                                                                                  cGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCCTCAATCGGTCCGGGc                                                                         >  2:334486/1‑251 (MQ=255)                                                                                                                                                     cAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTc                                                                      >  1:169799/1‑251 (MQ=255)                                                                                                                                                        aTTGCGCCCCGGCGCGGCCCTGGCGAGTCAAATGAAGGTTTTTGGCCCAAAACCA‑TTCGCTGAAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAg                                                                   <  1:11674/251‑1 (MQ=255)                                                                                                                                                               cACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTT‑CAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGTCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTAc                                                            <  1:482452/251‑1 (MQ=255)                                                                                                                                                                 cGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTg                                                          >  1:30256/1‑251 (MQ=255)                                                                                                                                                                    cgcTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTc                                                        >  1:82449/1‑250 (MQ=255)                                                                                                                                                                               cGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGAGCGGGCGTCGAGCCGGTACTGTCAACGCCAATccc                                            >  1:55103/1‑251 (MQ=255)                                                                                                                                                                                   tCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTc                                                                                                                                >  1:707175/1‑163 (MQ=255)                                                                                                                                                                                          cAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGCAGTCAACATTGCCCGGCGCGTGGCCCAAACGGGAAATCTCTTC‑GCTTCTTTCAACCGCAGCGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGGTTTCGGCAAACTCCGGgc                                                                                                                  >  1:407827/1‑170 (MQ=255)                                                                                                                                                                                             gATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTgcgc                              >  2:114230/1‑251 (MQ=255)                                                                                                                                                                                                    gCCCAGGACCA‑TTCGCTGGAATCAATACCGACACTATAGCGGGAGTAGTCAACATTGCACGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCCTCTTTCCACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTATCGGCAAACTACGGGCGCAGCGCCAGCAAGTTACCGTCGGTTTCAATCGGTGAGGGCACCGAGCCGGGCCTGTCAACGCCAATCCCGACCCCAGCTGCGCGCTGttt                       >  2:561342/1‑250 (MQ=255)                                                                                                                                                                                                      ccaaacccA‑TTCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAAc                                                                                                                         <  1:129280/145‑1 (MQ=255)                                                                                                                                                                                                      ccAGAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAAc                                                                                                                         >  2:129280/1‑151 (MQ=255)                                                                                                                                                                                                         gAACCA‑TTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGc                  >  1:591652/1‑251 (MQ=255)                                                                                                                                                                                                           acccttccgcggaaaaaaataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCt                 <  1:377896/235‑1 (MQ=255)                                                                                                                                                                                                                  cGCTGGAATAAATACCGCCAATATAGCGGGCGTAGTCAACATTGCCCGGCGCGCGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGCAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCCGCGTCAGCACGTTACCGTCGGCATCAATCGGTGCGGGTGTCGAGCAGGTACTGGCAACGCCAATCCCAACCACGGCTGCGCGCTCTTCCACGCTAAActct          >  2:240288/1‑251 (MQ=255)                                                                                                                                                                                                                          taaataCCGCCCATATAGCGGGAGTAGTCAACATTGCCCAGCGCGTGGCCCAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCCCGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTCCTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGCCGCTAAGCTCTGccagacc  >  2:353139/1‑245 (MQ=255)                                                                                                                                                                                                                           aaataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGt                                                                                             >  2:235697/1‑159 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     AGAAACTGGCTGGCGGCAGGCCGCCCCAGCTTTCGTGCCACAG‑AGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCA‑TTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTT‑CAACCGCAGTGTGGTCTTTCCACAATACGAAC‑ATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCAC  >  NC_000913/69382‑69847

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