Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,610,076	T→G	D202A (GAT→GCT)	uxaB ←	tagaturonate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,610,076	0	T	G	97.7%	149.5 / ‑5.0	43	D202A (GAT→GCT)	uxaB	tagaturonate reductase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (20/22);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.58e-01

ATTGCCACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGG  >  NC_000913/1609830‑1610319                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     aTTGCCACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAg                                                                                                                                                                                                                                                  <  2:124748/251‑1 (MQ=255)    gCCACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGCGCGGATATTTGTCCAGACGCAGTTAAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGagccagcc                                                                                                                                                                                                                                                >  2:411035/1‑250 (MQ=255)     ccACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCAt                                                                                                                                                                                                                                              >  1:168789/1‑251 (MQ=255)          ttGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAAt                                                                                                                                                                                                                                         >  1:313906/1‑251 (MQ=255)              ggtCTTTGTACGGTTTTATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAAt                                                                                                                                                                                                                                     <  1:57384/249‑1 (MQ=255)                tCTTTGTACTGTGTAATTGCGTTAACAATCATCACGTTGAGCGGATATTGGGCCAGACGCAGTTCAGTCGCTAAGGATTTCGTTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCTGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATg                                                                                                                                                                                                                                    <  1:216984/250‑1 (MQ=255)                          gTTTAAGATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCACCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGt                                                                                                                                                                                                                         <  2:345685/251‑1 (MQ=255)                                   cGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCActc                                                                                                                                                                                                                 <  1:122863/250‑1 (MQ=255)                                                    ttGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGTGGATAACCGGTAATGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCAc                                                                                                                                                                                               >  1:184960/1‑251 (MQ=255)                                                            atatTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCAc                                                                                                                                                                                       <  1:226409/251‑1 (MQ=255)                                                                 tGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAAc                                                                                                                                                                                  >  1:138110/1‑251 (MQ=255)                                                                         cGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCTCAGCACCAGTTCACGCAACGCGTC‑Gcc                                                                                                                                                                         >  2:281777/1‑251 (MQ=255)                                                                                               tCGGTCCCTGAATCACAAACATGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCAAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTGCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTcaca                                                                                                                                                    <  2:251510/250‑1 (MQ=255)                                                                                                  gTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGa                                                                                                                                                <  1:281777/251‑1 (MQ=255)                                                                                                       tGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCacac                                                                                                                                                   >  2:127539/1‑243 (MQ=255)                                                                                                       tGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCacac                                                                                                                                                   <  1:127539/243‑1 (MQ=255)                                                                                                                      tAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATc                                                                                                                            >  1:269384/1‑251 (MQ=255)                                                                                                                             gtTCAGCGGTGTCGAGAAACCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCa                                                                                                                     <  2:357220/251‑1 (MQ=255)                                                                                                                                cAGCTGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCg                                                                                                                  <  1:225337/251‑1 (MQ=255)                                                                                                                                          gagaAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGc                                                                                                        <  1:10202/251‑1 (MQ=255)                                                                                                                                              aaaCCAGCGTGATAACCCAGTTCTTCTACCAGTTTAGCCATTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATGGCGCAGCACCAGTGCTCGCAACGCGTC‑GCCATTGTAGTCAATCAACTCTCTCGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTaaa                                                                                                    <  1:200609/251‑1 (MQ=255)                                                                                                                                               aaCCAACGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCATCGTAGAACAGAAGCTGTTAGCCTGAGACAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAg                                                                                                   <  2:189807/251‑1 (MQ=255)                                                                                                                                                    cggTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCAACGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGtt                                                                                               <  2:353221/248‑1 (MQ=255)                                                                                                                                                      gTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGtt                                                                                               >  1:353221/1‑248 (MQ=255)                                                                                                                                                       tGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGa                                                                                            >  2:124270/1‑250 (MQ=255)                                                                                                                                                       tGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTAAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGTAATGGCTAAAGCGTTCGa                                                                                            >  2:318343/1‑250 (MQ=255)                                                                                                                                                        gATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAAc                                                                                                                                                         >  1:99618/1‑188 (MQ=255)                                                                                                                                                         aTAACCCAGTTCTACTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAAca                                                                                         <  2:269384/251‑1 (MQ=255)                                                                                                                                                                    cttcttcCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTc                                                                              >  1:379098/1‑251 (MQ=255)                                                                                                                                                                     ttcttcCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTcc                                                                             >  2:376511/1‑250 (MQ=255)                                                                                                                                                                      tcttcCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAgcg                                                                                                                                                                                                     <  2:293508/131‑1 (MQ=255)                                                                                                                                                                      tcttcCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAgcg                                                                                                                                                                                                     >  1:293508/1‑131 (MQ=255)                                                                                                                                                                                    cgcTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATTATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCATACGGGTCAGTTTTGCCGGATAGc                                                              <  1:393304/249‑1 (MQ=255)                                                                                                                                                                                           acgcgcgGATAACCGTTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTGAGCCTGAGCCAGCCAGTGAATAAATTCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCTTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATTGCTTACCGc                                                       <  1:53428/250‑1 (MQ=255)                                                                                                                                                                                                           tAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAAt                                       >  1:250621/1‑251 (MQ=255)                                                                                                                                                                                                            aaCGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGAATAGCTTACCGCTGGCGCGTCATCGAAtt                                      >  1:151258/1‑251 (MQ=255)                                                                                                                                                                                                                 tACGGTCCAACAGCGTAGAACAGAAGCTGTTAGCCTAAGCCAGACATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAAATCACACGGAATGATGATCCACCCTTTATCCATCGCACCGTAGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCAGGATAGCTTACCACTGACGCGTCAACGACTTTATCg                                 >  2:216082/1‑251 (MQ=255)                                                                                                                                                                                                                  aCGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCg                                 <  1:301542/250‑1 (MQ=255)                                                                                                                                                                                                                   cGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGcc                               >  1:62410/1‑251 (MQ=255)                                                                                                                                                                                                                      tCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGc                            >  1:148915/1‑251 (MQ=255)                                                                                                                                                                                                                                   aaCAGAAGCTGTTAGCCTGAGCCAGCCATTGCATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGAt                >  2:88972/1‑250 (MQ=255)                                                                                                                                                                                                                                             ggtAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGGAGTCAAGCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGGCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTc       <  2:122258/247‑1 (MQ=255)                                                                                                                                                                                                                                                aGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCCCTCTTGCGCATAGCGCAGCCCCAGTTCACGCAACGC‑TCTGCCATTGTAGTCAATCAACTCCCACGGAATGATGAGCCAACCTTTATCCAGCGCCCCGTTGAAATGGCTAAAGCGTTCGAACATCAGACGGGTCATTTTTGCCGGATAGCTTTCCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCGGATACCCGCTTCggtg   <  1:266938/250‑1 (MQ=255)                                                                                                                                                                                                                                                aGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGCTACCCGCTTCggtgg  >  1:124984/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     ATTGCCACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTC‑GCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGG  >  NC_000913/1609830‑1610319

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