Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,224,895	C→T	L755L (CTG→TTG)	ebgA →	evolved beta‑D‑galactosidase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,224,895	0	C	T	97.8%	147.2 / NA	45	L755L (CTG→TTG)	ebgA	evolved beta‑D‑galactosidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/19);  minor base A (1/0);  total (26/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 = 9.84e-01

ACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAGGTTAAC  >  NC_000913/3224658‑3225137                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   aCGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCGGTTGACTCGCGAGcc                                                                                                                                                                                                                                       <  1:236639/251‑1 (MQ=255)              aTTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAActtc                                                                                                                                                                                                                          <  1:9013/251‑1 (MQ=255)                ccccGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAActtctt                                                                                                                                                                                                                        <  2:111346/250‑1 (MQ=255)                                            acCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGc                                                                                                                                                                                              >  1:165394/1‑249 (MQ=255)                                                tATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACACCCACAAGCAGGAGt                                                                                                                                                                                        >  2:187923/1‑251 (MQ=255)                                                          tATCAGTCCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGAGCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAggg                                                                                                                                                                                <  2:41855/249‑1 (MQ=255)                                                           aTCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAg                                                                                                                                                                                                                                         <  2:98764/191‑1 (MQ=255)                                                           aTCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAg                                                                                                                                                                                                                                         >  1:98764/1‑191 (MQ=255)                                                             cAGTTCCCGCTGAAGGAAAACACCGCGCATCCAGTGCCTTTCGCACCAAATAATGCGCGCCCGCTGCCGCGGGAAGACGATCGTTTGAGCAGCACCGTTCGCGGCTACAACTTCGCGATCCCCTTCCCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGTCTCCCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCATGAGTACGAAGGGCt                                                                                                                                                                              <  2:235853/248‑1 (MQ=255)                                                               gTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACGCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGc                                                                                                                                                                         <  2:59876/251‑1 (MQ=255)                                                                     gCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCg                                                                                                                                                                    >  2:167921/1‑250 (MQ=255)                                                                       tGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAAt                                                                                                                                                                 >  1:154266/1‑251 (MQ=255)                                                                        gAAGGAAAGCACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATc                                                                                                                                                                >  1:164661/1‑251 (MQ=255)                                                                             aaaaCACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTg                                                                                                                                                           >  1:213669/1‑251 (MQ=255)                                                                              aaaCACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGCGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCACGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGGGGCAACCGAATCATTTGc                                                                                                                                                          >  1:205148/1‑251 (MQ=255)                                                                              aaaCACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCAGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGTCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGc                                                                                                                                                          >  2:67040/1‑251 (MQ=255)                                                                               agcacCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCa                                                                                                                                                         <  2:164661/249‑1 (MQ=255)                                                                               aaCACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGc                                                                                                                                                          >  2:55040/1‑250 (MQ=255)                                                                                               tGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACAt                                                                                                                                         >  1:300131/1‑251 (MQ=255)                                                                                                   tttCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACACCTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGCGAATGGCGAAACGCTGTTGACTCGCGAGCCAAAGATCAACGTCTTCAAGCCGATGATCGCCAACCACAAGAAGGAGTACGCAGTGCTGGGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTgc                                                                                                                                     >  2:250519/1‑251 (MQ=255)                                                                                                    ttCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTgcg                                                                                                                                    >  2:202924/1‑251 (MQ=255)                                                                                                        cACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGAc                                                                                                                                >  1:24160/1‑251 (MQ=255)                                                                                                              ataatGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGcc                                                                                                                          <  2:246773/251‑1 (MQ=255)                                                                                                                aatGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGt                                                                                                                        >  1:208099/1‑251 (MQ=255)                                                                                                                         ccGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACag                                                                                                                 >  1:58163/1‑249 (MQ=255)                                                                                                                          cGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGc                                                                                                              >  1:244770/1‑251 (MQ=255)                                                                                                                                     gAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTc                                                                                                                                                                                                                                              <  1:220450/112‑1 (MQ=255)                                                                                                                                     gAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTc                                                                                                                                                                                                                                              >  2:220450/1‑112 (MQ=255)                                                                                                                                           gATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGat                                                                                             <  1:69678/251‑1 (MQ=255)                                                                                                                                           gATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGat                                                                                             >  2:23675/1‑251 (MQ=255)                                                                                                                                                 ttGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGAtcatca                                                                                        <  2:73462/250‑1 (MQ=255)                                                                                                                                                        gCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCCAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGcac                                                                                 >  1:16357/1‑250 (MQ=255)                                                                                                                                                          aCCGTTCGCGGCTACAACTTAGCGATCACCTTCACCAAAATAAGTGGCAAACCGACAACCTGGCAGGCGAATGGCGAATCGATGATGCCTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGAGGCAACCGAATCATTTGCAGATCATGCAGGAACATATGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGt                                                                              >  1:23284/1‑251 (MQ=255)                                                                                                                                                            cGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTa                                                                            <  1:219979/251‑1 (MQ=255)                                                                                                                                                                   ggCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGcccc                                                                     <  2:165394/251‑1 (MQ=255)                                                                                                                                                                   ggCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGcccc                                                                     <  2:184890/251‑1 (MQ=255)                                                                                                                                                                     cTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTAGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCccgc                                                                   >  2:301950/1‑251 (MQ=255)                                                                                                                                                                      tACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCccgcc                                                                  >  1:226099/1‑251 (MQ=255)                                                                                                                                                                                    ccccttcaaaaaaaTGGGGGGAAAACCGACACCCTGGCGGGTAAAGGGCGAACCGCTGTTGACCCGCGAGCAAAAGATCACCTCCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGGGAAGTGCTGATCTTCACCCGCAAAGTTATTCCCCCGCCGGTGTTGGACTTCg                                                    <  2:16357/242‑1 (MQ=255)                                                                                                                                                                                                               gacacccgggCAGGTGAAGGGCGAATCGCTGTTGACTCGCGAGCAAAAGATCAACTTCTTCAAGCCGATGATCGACACCCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGACCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGACCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCa                         <  2:244770/243‑1 (MQ=255)                                                                                                                                                                                                                    ccTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCt                    >  1:117897/1‑251 (MQ=255)                                                                                                                                                                                                                       ggaaGGTGAATGGCCAACCGCTGTTGCCTCGGGAGCAAAAGTTAACCTTCTCCAAGCCGAGGATCGACACCCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGcc                 <  2:300131/248‑1 (MQ=255)                                                                                                                                                                                                                               aaTGGCAAATCGCTGTGGACTCGCAAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCa         <  1:187923/251‑1 (MQ=255)                                                                                                                                                                                                                                aTGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGcc          >  2:127133/1‑249 (MQ=255)                                                                                                                                                                                                                                      aaTCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAGGTTAAc  >  1:20613/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   ACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAGGTTAAC  >  NC_000913/3224658‑3225137

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