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	98.4%	206.8 / ‑4.3	63	G148G (GGT→GGC)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (40/21);  total (41/22)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

AGGATCGGGTCCAGCTCATCAAAGAAACTGGCTGGCGGCAGGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTT‑CAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGG  >  NC_000913/69360‑69849                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      aGGATCGGGTCCAGCTCCTCAAAGAAACTGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCCGA‑CGCCCAGCTCAATCCCCGATGCGGCAGAGTGCCCCACGGCCCTGTCCTGGGGGGTCACATGCAGG‑GTTTTTGCCCAGAACCCTTGCCTGGACTAAATACCCCCAAt                                                                                                                                                                                                                                                  >  2:642246/1‑251 (MQ=255)  ggATCGGGTCCAGCCCATCAAAGAAACTGGCTGGCGGCAAGCCGCCACAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGACGAGTCACATGCAGG‑ATTTTTGCCCAGAACCCTTCGCTGGAATAAATACCGCCAAt                                                                                                                                                                                                                                                  >  1:641292/1‑250 (MQ=255)   gATCGGGTCCAGCTCATCAAAGAAACTGGCTGGCGGCAAGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAAtat                                                                                                                                                                                                                                                >  1:342718/1‑249 (MQ=255)                       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‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCTTCAATCGGTGCGGGCGTc                                                                        >  2:93996/1‑251 (MQ=255)                                                                                                                                                                                  ggccgcGGCGTTGTCGTGGCTGGTCACAGGCAGG‑TTTTTTGCCCAGACCCATCCGCGGAAATAAATCCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCGCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCgg                                                                 <  1:789772/246‑1 (MQ=255)                                                                                                                                                                                          cgcTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGAGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATAAATCGGTGCGGGCGTCGAGCCGGTACTGTc                                                          >  1:503330/1‑250 (MQ=255)                                                                                                                                                                                           gcTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGAGTGGCACAAACGGGTAATCTCTTC‑GCGTCTTTCAACAGCAGGGTGGGATTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCAAGTTACCGTGGGAATCAATCGGGCCGGGCGTCTAGCCGGTCCTggca                                                         >  2:209287/1‑247 (MQ=255)                                                                                                                                                                                              gTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACAGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACg                                                      <  1:530673/250‑1 (MQ=255)                                                                                                                                                                                                  tggcggggcacaTGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCGATATAGCGGGAGTAGACAACCTTGCCCGGCGCGTGGAACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGGGTTCGGGTTTTCGGCAAACTCAGGGCGCAGCGCAAGCAGGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGAACTGTCAACGcgca                                                  >  2:41344/9‑247 (MQ=255)                                                                                                                                                                                                     cGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATccc                                              >  1:887059/1‑251 (MQ=255)                                                                                                                                                                                                            aaTGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATCAATACCGCCAATCTAGCGGGAATAGTCAACATTGCCCGGCGCGAGGCACAAACGGGTCATCTCTTC‑GCATCTTTAAACCGCAGTGTGGTCTTTCCCCAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCCTCAATCGGTGCGTGCGTCGAGCCGGAACAGTCAACTCCAATCCCGCCcac                                        >  2:849121/2‑250 (MQ=255)                                                                                                                                                                                                              tGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCAGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGc                                     >  1:957688/1‑251 (MQ=255)                                                                                                                                                                                                                 aGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTgc                                  >  1:787385/1‑251 (MQ=255)                                                                                                                                                                                                                       ttttGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCt                             >  1:647870/1‑251 (MQ=255)                                                                                                                                                                                                                           tcccaaaacctgtcgcggaaaaaaataCCGCCAATATAGCGGGAGTAGGCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGTCGTCGAGCCGGTACTGGCAAAGCGAATCCCGACCACAGCTGCGCGCTGtt                          <  1:803287/228‑1 (MQ=255)                                                                                                                                                                                                                                 aaCCATTCCCTGGACTAAAAACCGCCAATATAGCGGGAGTAGTCAACATTCACCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTACACAATCCGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGCCGCAGCGCCAGCCCGGTACCGTCGGCATCCACAGGCGAGGGCGTAGCGCCGGTACTGTACACACCAATCCCGCCcaacgatgcgcgcagtacgaagcc                   >  2:458370/1‑229 (MQ=255)                                                                                                                                                                                                                                   ccttctgttgaaaaaaataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTaa                 <  1:978456/237‑1 (MQ=255)                                                                                                                                                                                                                                   ccATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTa                  >  2:756455/1‑250 (MQ=255)                                                                                                                                                                                                                                     aTTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCAAGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGCTCGACGCTAAGc               >  1:745034/1‑251 (MQ=255)                                                                                                                                                                                                                                                  aataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACgg  <  1:668023/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      AGGATCGGGTCCAGCTCATCAAAGAAACTGGCTGGCGGCAGGCCGCCCCAGCTTTCGTGCCACAGAGATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGT‑CGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGG‑ATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTT‑CAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGG  >  NC_000913/69360‑69849

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