Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,611	G→A	Y146Y (TAC→TAT)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,611	0	G	A	98.4%	200.2 / NA	64	Y146Y (TAC→TAT)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (40/21);  minor base C (1/0);  total (41/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AGAAACTGGCTGGCGGCAGGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTT‑CAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTC  >  NC_000913/69382‑69854                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       ataaaCGGGCTGGCGGCAAGCCCCCCCACCTTGCGTGCCACAGAGTTTTATGCCCGGCGCTGCCACGTCCTCGTCAAATTTCCTGCGGGCGGGTGTTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCTGGAGGGGGCAGATTGCGCCACGGCGCTGTCCTGGCGATTCACA‑GCGGGTATTTTTCCCAAGAACAATTCGCTGGAAAAATTCCCGCCAATATAGCGGGAGTAGTCAACATTg                                                                                                                                                                                                                                   <  1:536245/249‑1 (MQ=255)     acgggcgggcggaAGGCCGCCCCAGCTTTCGTGCCACAGAAATTTATGCCCGGCGCGGCAACGCCCGCGACGAATATCCGGGGGGGGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCGGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCg                                                                                                                                                                                                                               <  1:665064/244‑1 (MQ=255)       gggctggcGGCAAGCCGCCCATTCTTTCGTGCCACAGAGTTTTTTGCCCGGCGCTTCAACGTCCGCGACAAATATCCGGCGGGCGGGTGGCCCCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAACCCACGATGCGGAAGATTGCGCCACGGCGCTGTCCTGGCGAGCCACATGCAGG‑ATTTTTGCCCAGACCCTTTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGgc                                                                                                                                                                                                                             <  1:588295/250‑1 (MQ=255)        ggctggcGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGgcg                                                                                                                                                                                                                            >  1:754382/1‑251 (MQ=255)               ggcAGGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCTCATGCAGG‑ATTTTTGCCCAGAACCATTCGCTAGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGca                                                                                                                                                                                                                     <  1:405960/251‑1 (MQ=255)                         cccAGCTTTCGTGCCACAGAGATTTAGGCCGGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACGGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAACCCCCGAGGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGACCCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTa                                                                                                                                                                                                           <  1:309263/251‑1 (MQ=255)                           cAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAAt                                                                                                                                                                                                         >  1:464959/1‑251 (MQ=255)                                 ccgggccccaGAGTTTTAGGCCCGGCGCGGCACCGTCCGCGACGAATACCCTGCGGGCGGGTGGTCCCGGAAAGCAGAGCTGGCCCCCAGT‑CGCACAGCTCATCCCACGATGCGGAAGATTGCGCCACGGCGCTGTCCGGGCGAGTCACATGCAGG‑ATTTTTGCCCAGACCCATCCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCtt                                                                                                                                                                                                    <  1:270117/242‑1 (MQ=255)                                              aTTTGTGCCCGGCGCTGCAACGTCCGGTACAAATACCCTGGGGGCGGGGGGTCCCGGAAAGCAGGGCGGGCCCCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGTTTGCGCCCGGGCGCTGTCCTGGCGAGTCACATGCAGG‑TTTTTTCCCCAGAACAATTCGGTGGAATAAATCCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAAcc                                                                                                                                                                                     <  1:385782/251‑1 (MQ=255)                                                ttaggCCCGGCGCGCCAACGCCCGCGACGAATACCCGGCGGGCGGGTGGTACCGAAAAGCAGAGCGGGCCCCCGGT‑CGCACAGCCCAACCCAGGAGGCGGAAGATGGCGCCACGGCGCTGTCCTGGCGAGTCACAGGCAGG‑ATTTTTGCCCAGAACCATCCGCGGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCg                                                                                                                                                                                    <  1:439804/246‑1 (MQ=255)                                                ttATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATATCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGc                                                                                                                                                                                   >  1:213160/1‑251 (MQ=255)                                                     cccGGCGCTGCACCGTCCGCGACGAATTTCCGGCGGGCGGGTGGTACGGAAAAGCAGAGGTGCCCCCCAGC‑CCCCCAGCTCAATCCCCGAGGCGCCAGATGGCGCCACGGCGCGGTCCTGGCGAGTCACATTCAGG‑ATTTTTCCCCAAACCCATCCGCTGAAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAgtgt                                                                                                                                                                              <  1:812523/251‑1 (MQ=255)                                                        ggCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGt                                                                                                                                                                           >  1:411179/1‑251 (MQ=255)                                                        ggCGCTGCAACGCCCGCGACAAATTCCCGGGGGGCGGGTGGTACGGAAAAGAAGAGCTGGCACCATGT‑CGCACAGCACAACCCACGTTGCGGCAGATTGCCCCCGGGCGCTGTCCGGGCGAGCCACAGGCAGG‑TTTTTGGCCCAGAACCATTCGCGGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGt                                                                                                                                                                           <  1:746724/251‑1 (MQ=255)                                                         gcgcTGCAAGGTCCGCGACGAATACCCGGGGGGGGGGTGGTCCCGGAAAGCAGAGCGGGCCCCCAGT‑CGAACAGCTAAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCAAGAACCATTCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTc                                                                                                                                                                          <  1:963112/251‑1 (MQ=255)                                                          cgcTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCt                                                                                                                                                                         >  2:422019/1‑251 (MQ=255)                                                            cTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCACCCGCAGTGTGGTCttt                                                                                                                                                                       >  2:399096/1‑251 (MQ=255)                                                                 cggccggcgacaattccccggggggcggtTGGTGCCGGAAGGCGGAGCTGG‑ACCCCGTGCGCACATCTCAATCCCCGATGCGGAAGATGGCGCCACGGCGCTGTCGGGGCTAATCACTTGCAAG‑ATTTGTGCCCAAACCCTTTCGCGGAAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCT‑GCGCTTCTTTCAACCGCAGTGTGGTCTTTCcaca                                                                                                                                                                  <  2:714051/229‑1 (MQ=255)                                                                 aCGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCcac                                                                                                                                                                   >  1:424484/1‑250 (MQ=255)                                                                 aCGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGAAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCcac                                                                                                                                                                   >  1:296258/1‑250 (MQ=255)                                                                   gTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑CTTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCACCCGCAGTGTGGTCTTTCCACaa                                                                                                                                                                 >  2:398015/1‑250 (MQ=255)                                                                   gTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCAGGCGCGTGGCACAACCGGGAAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCcaaaa                                                                                                                                                                 >  1:784929/1‑247 (MQ=255)                                                                        ggacgaATAGCCTGCGGGCGGGTTGTCCCGAAAAGCAGTGCTGCCACCCCGC‑CGCCCTTCTAAATCCATGATCCGGCAGATTGCGCCACGGCGTTGTCCGGGCAAGTCACATGCAGT‑ATTTTTCCCCATAACCATTCGCTTGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGaa                                                                                                                                                           <  1:679925/250‑1 (MQ=255)                                                                          acgaATACCCGGCGGGGGGGTGATACCGGACAGCAGAGCGGGCCCCCAGG‑CGCACAGCCAAACCCACGATGCGGCAGATTGCGCCCCGGCGCGGTCCTGGCGAGTCACATTGAGG‑ATTTTTGCCGAGACCCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACa                                                                                                                                                         <  1:377864/251‑1 (MQ=255)                                                                            gaATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATc                                                                                                                                                       >  1:975851/1‑251 (MQ=255)                                                                                 accTGCGGGGGGGTGGTACCGGAAAACAGGGCGGGCCCCCCGG‑CGCAAACCACTACCCCCGATGGGGCAGATTGCCCCCCGGCTCTGCCCGGGCGAGTCACATTCAGG‑TTTTTTGCCCAAAACCATTCGCTGGAAAAAATCCCGCCAATTTAGCGGGAGTAGTCAACATTGCCCGGCGCCTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGt                                                                                                                                                   <  1:706884/249‑1 (MQ=255)                                                                                   cTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGCATAAATACCGCCAATATAGCGGGCGTAGGCAACATTGCCCGGCGCATGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGGGTTGTCTTTCCACAATACGATCACCGCGTTcc                                                                                                                                                >  2:235364/1‑250 (MQ=255)                                                                                      cgggcgggTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCATGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCCCTTC‑GCTTCTTTCAACCGCAGAGTGGGCTTTCCACAATACGAACGTCGCGTTCGGGt                                                                                                                                             >  1:565438/1‑251 (MQ=255)                                                                                        ggggggTGGAACCGGAAAGCAGGGGTGGCACCCAGG‑CGCACAGCTCAAGCCACGATGCGGCAGATTGCGCCCCGGCGCTGTCCTGGCGAGTCACAGGCAGG‑ATTTTGGCCCAGAACCATTCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGttt                                                                                                                                           <  1:412543/248‑1 (MQ=255)                                                                                         gcgggTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATCCGAACATCGCGTTCGGGtttt                                                                                                                                          >  2:953606/1‑251 (MQ=255)                                                                                             gtggtACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGc                                                                                                                                      >  1:629577/1‑251 (MQ=255)                                                                                               ggtACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCa                                                                                                                                     >  1:667874/1‑250 (MQ=255)                                                                                                      gAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCgg                                                                                                                             >  2:293988/1‑251 (MQ=255)                                                                                                       aaaGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCggg                                                                                                                            >  1:295173/1‑251 (MQ=255)                                                                                                          gCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGCACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATCCGAACATCGCGTTCGGGTTTTCGGCAAAATCCGTgcgc                                                                                                                         >  2:763796/1‑251 (MQ=255)                                                                                                            agagCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGgcgc                                                                                                                         >  1:993362/1‑249 (MQ=255)                                                                                                                 tGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGcc                                                                                                                   >  1:548974/1‑250 (MQ=255)                                                                                                                      cccAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACg                                                                                                             >  1:690433/1‑251 (MQ=255)                                                                                                                               cacaGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTA‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACACTACGAACATCGCGTTCGGGTTTTCGGCACACTCCGGGCGCAGCGCCAGCACGTTACCGTc                                                                                                     >  1:147477/1‑251 (MQ=255)                                                                                                                                acaGCCCAATCCACGAGGCGGCAGATGGCCCCCCGGCGCTGTCCGGGCGAGTCACATGCAGG‑ATTTTTGCCAAAACCCATCCGCGGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCg                                                                                                    <  1:293988/251‑1 (MQ=255)                                                                                                                                         tCCACGATGCGGCAGATTGAGACACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATACATACCGCCAATATAGCGGGCGTAGTCAACATTGCCCGGCGCGTGGCACAAAAGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGGTTTCGGCAAACTACGGGCGCAGCGACAGCACggtaccggaggcctccatc                                                                                           >  2:897329/1‑233 (MQ=255)                                                                                                                                              gATGCGGCAGAATGAGAAACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATGCGCTGGAATAAATACAGCCGATATAGCGGGAGTAGTCAACATTGCCCGGGGCGAGACACAAACGGTTCATCTATTC‑GCATCTTTCAACAGAAGTGTGCACTTTTCCCAAAACGAAAATCGTGTTAGGTGTTTCGTGAAACTATGGACGCAGCGCCAGAACGTTTACGTCGGGCTCACTCGGTGc                                                                                      >  2:745035/1‑251 (MQ=255)                                                                                                                                                    gCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACAGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCATCACGTTACCGTCGGCATCAATCGGTGCGGGCgg                                                                                >  2:530673/1‑250 (MQ=255)                                                                                                                                                     cAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCTTCAATCGGTGCGGGCG‑Tc                                                                              >  2:93996/1‑251 (MQ=255)                                                                                                                                                            ggccgcGGCGTTGTCGTGGCTGGTCACAGGCAGG‑TTTTTTGCCCAGACCCATCCGCGGAAATAAATCCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCGCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCgg                                                                       <  1:789772/246‑1 (MQ=255)                                                                                                                                                                    cgcTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGAGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATAAATCGGTGCGGGCG‑TCGAGCCGGTACTGTc                                                                >  1:503330/1‑250 (MQ=255)                                                                                                                                                                     gcTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGAGTGGCACAAACGGGTAATCTCTTC‑GCGTCTTTCAACAGCAGGGTGGGATTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCAAGTTACCGTGGGAATCAATCGGGCCGGGCG‑TCTAGCCGGTCCTggca                                                               >  2:209287/1‑247 (MQ=255)                                                                                                                                                                        gTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACAGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACg                                                           <  1:530673/250‑1 (MQ=255)                                                                                                                                                                            tggcggggcacaTGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCGATATAGCGGGAGTAGACAACCTTGCCCGGCGCGTGGAACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGGGTTCGGGTTTTCGGCAAACTCAGGGCGCAGCGCAAGCAGGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGAACTGTCA‑ACGcgca                                                       >  2:41344/9‑247 (MQ=255)                                                                                                                                                                               cGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATccc                                                   >  1:887059/1‑251 (MQ=255)                                                                                                                                                                                      aaTGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATCAATACCGCCAATCTAGCGGGAATAGTCAACATTGCCCGGCGCGAGGCACAAACGGGTCATCTCTTC‑GCATCTTTAAACCGCAGTGTGGTCTTTCCCCAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCCTCAATCGGTGCGTGCG‑TCGAGCCGGAACAGTCA‑ACTCCAATCCCGCCcac                                             >  2:849121/2‑250 (MQ=255)                                                                                                                                                                                        tGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCAGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGc                                          >  1:957688/1‑251 (MQ=255)                                                                                                                                                                                           aGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTgc                                       >  1:787385/1‑251 (MQ=255)                                                                                                                                                                                                 ttttGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCt                                  >  1:647870/1‑251 (MQ=255)                                                                                                                                                                                                     tcccaaaacctgtcgcggaaaaaaataCCGCCAATATAGCGGGAGTAGGCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGTCG‑TCGAGCCGGTACTGGCA‑AAGCGAATCCCGACCACAGCTGCGCGCTGtt                               <  1:803287/228‑1 (MQ=255)                                                                                                                                                                                                           aaCCATTCCCTGGACTAAAAACCGCCAATATAGCGGGAGTAGTCAACATTCACCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTACACAATCCGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGCCGCAGCGCCAGCCCGGTACCGTCGGCATCCACAGGCGAGGGCG‑TAGCGCCGGTACTGTAC‑ACACCAATCCCGCCcaacgatgcgcgcagtacgaagcc                        >  2:458370/1‑229 (MQ=255)                                                                                                                                                                                                             ccttctgttgaaaaaaataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTaa                      <  1:978456/237‑1 (MQ=255)                                                                                                                                                                                                             ccATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTa                       >  2:756455/1‑250 (MQ=255)                                                                                                                                                                                                               aTTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCAAGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGCTCGACGCTAAGc                    >  1:745034/1‑251 (MQ=255)                                                                                                                                                                                                                            aatcccgccaATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACgg       <  1:347679/244‑1 (MQ=255)                                                                                                                                                                                                                            aataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACgg       <  1:668023/251‑1 (MQ=255)                                                                                                                                                                                                                               accgccaATATAGCGGGCGGAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATAGCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCCCAAGACGAACATCGGTTTCGGATTTTCGGAAAACTCGGTGCGCAGCCCAAGCACGGTACAGTCGCCTTCAATCGACGCGG‑CGATCGAGCCGGAACTG‑CACACGCCAGTCCCTACCACATCGGCGCCCGGTTCGCGGCTAAGCTCAGCAAACACGGttt    >  2:587592/5‑251 (MQ=255)                                                                                                                                                                                                                                ccgccaATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGtttt   >  1:691152/4‑251 (MQ=255)                                                                                                                                                                                                                                 cgccaATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTc  >  1:270450/3‑251 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       AGAAACTGGCTGGCGGCAGGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTT‑CAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCG‑TCGAGCCGGTACTGTCA‑ACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTC  >  NC_000913/69382‑69854

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