Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,171	G→A	G389G (GGC→GGT)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,171	0	G	A	97.6%	138.2 / NA	42	G389G (GGC→GGT)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/19);  minor base C (1/0);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTT  >  NC_000913/66932‑67408                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              tGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGACGAAGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGCAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGATTGGGt                                                                                                                                                                                                                                    >  1:320392/1‑251 (MQ=255)                   cgcATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGc                                                                                                                                                                                                                  >  2:184629/1‑251 (MQ=255)                    gcATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCAACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTTGCCCGGTTTGGGTATTGAAGATCAGTCGGGcc                                                                                                                                                                                                                 >  2:282944/1‑250 (MQ=255)                      atatCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAgg                                                                                                                                                                                                               >  1:233169/1‑251 (MQ=255)                       tatCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAgg                                                                                                                                                                                                               >  2:210089/1‑250 (MQ=255)                                        gCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAAt                                                                                                                                                                                             >  1:31717/1‑251 (MQ=255)                                            ggCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTTGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATAcc                                                                                                                                                                                          <  1:184629/250‑1 (MQ=255)                                             gCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCg                                                                                                                                                                                         >  1:48745/1‑250 (MQ=255)                                                   gACGGTATCGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCt                                                                                                                                                                                  <  1:304031/251‑1 (MQ=255)                                                    aCGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTg                                                                                                                                                                                 >  1:12085/1‑251 (MQ=255)                                                          tGGTGCGCGCCACCAGCGCGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGGTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTc                                                                                                                                                                           >  1:58385/1‑251 (MQ=255)                                                                   ccaccaGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGAGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGGGGATCgg                                                                                                                                                                  >  2:220313/1‑251 (MQ=255)                                                                      ccaGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGttt                                                                                                                                                               <  1:242754/251‑1 (MQ=255)                                                                                              gCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGTAATCGACGGGCa                                                                                                                                       <  2:200112/251‑1 (MQ=255)                                                                                                    ggCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACtt                                                                                                                                  <  1:210089/250‑1 (MQ=255)                                                                                                           ccGGTTGCGCTTGCCACAGCGCATTCGCCCCCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCat                                                                                                                           <  2:204025/250‑1 (MQ=255)                                                                                                               ttGCGCTTTCCACAGAGCATTCGCCACCGGCAGTTTCGGCCGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATggg                                                                                                                      >  2:277989/1‑251 (MQ=255)                                                                                                                tGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATggg                                                                                                                      >  1:99167/1‑250 (MQ=255)                                                                                                                            aGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAg                                                                                                         <  1:229452/251‑1 (MQ=255)                                                                                                                               gcATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCCCCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGCTCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTc                                                                                                      >  2:303278/1‑251 (MQ=255)                                                                                                                                    cGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTa                                                                                                  <  1:138804/250‑1 (MQ=255)                                                                                                                                     gCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCCGCATCTGGGAGACGAGCACCAGGTCAttgcc                                                                                                >  2:171681/1‑248 (MQ=255)                                                                                                                                                     gCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGAt                                                                                  >  1:320629/1‑249 (MQ=255)                                                                                                                                                               cGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGtcct                                                                      <  2:95748/251‑1 (MQ=255)                                                                                                                                                                ggTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCGGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGGTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGtcctc                                                                     <  2:84934/251‑1 (MQ=255)                                                                                                                                                                 gtgtTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATAAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGtcctcc                                                                    >  1:291509/1‑251 (MQ=255)                                                                                                                                                                 gggtTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCGCCAGGTCATTTCCTTTCTCGAAGTTACAGGTGTAGtcctcc                                                                    <  2:144387/249‑1 (MQ=255)                                                                                                                                                                   gttgtCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCa                                                                   <  1:261765/246‑1 (MQ=255)                                                                                                                                                                               gcGATACAGTTAACCAGCAGACGGTAAAGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTTATTACCTTTCGCGAAGGGATAGGTGTAGTCCTCCATAAAGGAGGTgcc                                                      <  1:308732/250‑1 (MQ=255)                                                                                                                                                                                          aaCCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGa                                            <  1:54290/250‑1 (MQ=255)                                                                                                                                                                                             cagcaGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCgg                                        <  1:219762/247‑1 (MQ=255)                                                                                                                                                                                               gcaGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGCATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGCGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATACAGGAGGTGCCGCCCTGCAGCCCGGtt                                      >  2:56018/3‑251 (MQ=255)                                                                                                                                                                                                 aGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGa                                    <  1:245600/251‑1 (MQ=255)                                                                                                                                                                                                  gACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCGGCAGACCGGTTGAc                                   <  2:240240/251‑1 (MQ=255)                                                                                                                                                                                                   aCGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATTTGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGGGATAGGAGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGAc                                   >  1:209825/1‑250 (MQ=255)                                                                                                                                                                                                     ggTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATc                                >  1:60906/1‑251 (MQ=255)                                                                                                                                                                                                       tAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCAc                              >  2:42078/1‑251 (MQ=255)                                                                                                                                                                                                          cGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCtt                           >  1:146176/1‑251 (MQ=255)                                                                                                                                                                                                              cGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCat                        <  1:284497/250‑1 (MQ=255)                                                                                                                                                                                                                    gaTCAATCAGGCTGGCGCCAATCGCTGGACCGGCTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGCTGACATCACCTTCATGATGCg                  >  2:188307/1‑250 (MQ=255)                                                                                                                                                                                                                         atcagGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGc              >  1:226275/1‑249 (MQ=255)                                                                                                                                                                                                                                   cGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGtt  <  2:45018/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              TGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTT  >  NC_000913/66932‑67408

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