Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,117,903	Δ1 bp	5.9%	coding (502/1281 nt)	wcaK ←	colanic acid biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,117,903	0	G	.	5.9%	223.0 / 7.7	68	coding (502/1281 nt)	wcaK	colanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base . (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

CTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCG  >  NC_000913/2117756‑2118048                                                                                                                                                    |                                                                                                                                                  cttctGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTg                                                                                                                                                   >  4:272775/1‑148 (MQ=255)   tctGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACTTAGTTCGCCAGTTGGt                                                                                                                                                 >  1:225949/1‑148 (MQ=255)   tctGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGt                                                                                                                                                 >  6:221896/1‑148 (MQ=255)    ctGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGt                                                                                                                                                 >  8:233188/1‑147 (MQ=255)       tgtgGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTaa                                                                                                                                              <  3:117019/147‑1 (MQ=255)       tgtgGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTaa                                                                                                                                              >  4:117019/1‑147 (MQ=255)       tgtgGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACAGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTaa                                                                                                                                              <  6:162463/147‑1 (MQ=255)            tgATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCt                                                                                                                                        <  2:225949/148‑1 (MQ=255)               tCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCAt                                                                                                                                     <  8:32007/148‑1 (MQ=255)                 gACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATcc                                                                                                                                   <  2:354538/148‑1 (MQ=255)                  aCCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCt                                                                                                                                  <  6:192019/148‑1 (MQ=255)                   ccAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTg                                                                                                                                          <  3:94130/139‑1 (MQ=255)                   ccAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCACCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTGGGTTAAATGGCTCATCCTg                                                                                                                                 >  4:390122/1‑148 (MQ=255)                   ccAGCAACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTg                                                                                                                                          >  4:94130/1‑139 (MQ=255)                           gcgGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGc                                                                                                                         >  6:3284/1‑148 (MQ=255)                                  cGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACt                                                                                                                  >  2:153145/1‑148 (MQ=255)                                  cGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACt                                                                                                                  >  5:410712/1‑148 (MQ=255)                                          tGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCaa                                                                                                          >  6:329757/1‑148 (MQ=255)                                           gTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGTATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATcc                                                                                                                                   <  5:70259/122‑1 (MQ=255)                                           gTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGTATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATcc                                                                                                                                   >  6:70259/1‑122 (MQ=255)                                             tCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAac                                                                                                                     >  8:209933/1‑134 (MQ=255)                                             tCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAac                                                                                                                     <  7:209933/134‑1 (MQ=255)                                             tCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCa                                                                                                       >  2:218238/1‑148 (MQ=255)                                             tCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATAGCTCATCCTGGAACGGGCCAACACTGAGCCCAATCa                                                                                                       >  2:251557/1‑148 (MQ=255)                                                aCTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATaa                                                                                                    <  7:233188/148‑1 (MQ=255)                                                             tgATATTGCTGCGTTTCCTCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCttttt                                                                                       <  2:74827/148‑1 (MQ=255)                                                             tgATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCttttt                                                                                       <  3:86669/148‑1 (MQ=255)                                                                  ttGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGcc                                                                                   <  5:3284/147‑1 (MQ=255)                                                                  ttGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCa                                                                                  <  6:371965/148‑1 (MQ=255)                                                                     ctgcGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCATCACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATaa                                                                               <  1:251557/148‑1 (MQ=255)                                                                       gcGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATcc                                                                                                                                   <  5:9677/94‑1 (MQ=255)                                                                       gcGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATcc                                                                                                                                   >  6:9677/1‑94 (MQ=255)                                                                        cGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACg                                                                            <  6:410712/148‑1 (MQ=255)                                                                              atcaGAGCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCaa                                                                       <  5:221896/147‑1 (MQ=255)                                                                               tcaGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACg                                                                            <  6:35671/141‑1 (MQ=255)                                                                               tcaGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACg                                                                            >  5:35671/1‑141 (MQ=255)                                                                                 aGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTg                                                                   <  3:272775/148‑1 (MQ=255)                                                                                 aGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTg                                                                   >  1:213274/1‑148 (MQ=255)                                                                                  gATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCGTTTTCGCCATAAACGTGCAAAGTGc                                                                  <  2:393257/148‑1 (MQ=255)                                                                                   aTCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGa                                                                                                                               <  8:102044/86‑1 (MQ=255)                                                                                   aTCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGa                                                                                                                               >  7:102044/1‑86 (MQ=255)                                                                                       aaGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGtt                                                             <  1:153145/148‑1 (MQ=255)                                                                                             aCCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACt                                                       >  1:422188/1‑148 (MQ=255)                                                                                             aCCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACt                                                       >  3:384973/1‑148 (MQ=255)                                                                                               cGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGc                                                     <  3:153482/148‑1 (MQ=255)                                                                                                           gATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCa                                            >  6:99608/1‑145 (MQ=255)                                                                                                            aTCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATaga                                        <  5:329757/148‑1 (MQ=255)                                                                                                                          gTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAAc                          >  1:412091/1‑148 (MQ=255)                                                                                                                            gACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTg                                                                          <  7:190511/98‑1 (MQ=255)                                                                                                                            gACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTg                                                                          >  8:190511/1‑98 (MQ=255)                                                                                                                            gACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACg                         >  4:181804/1‑147 (MQ=255)                                                                                                                             agcAAAAACGTAGTTCGCCAGTTTGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGAt                       <  3:259426/146‑1 (MQ=255)                                                                                                                             aCCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGAt                       <  5:271494/148‑1 (MQ=255)                                                                                                                              ccAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATc                      >  4:392641/1‑148 (MQ=255)                                                                                                                              ccAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATc                      >  5:369761/1‑148 (MQ=255)                                                                                                                                 aaaaCGTAGTTCGCCAGTT‑GTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACaaaaa                            <  2:360109/138‑1 (MQ=255)                                                                                                                                 aaaaCGTAGTTCGCCAGTT‑GTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACaaaaa                            <  2:354106/138‑1 (MQ=255)                                                                                                                                 aaaaCGTAGTTCGCCAGTT‑GTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACaaaaa                            >  1:360109/1‑138 (MQ=255)                                                                                                                                 aaaaCGTAGTTCGCCAGTT‑GTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACaaaaa                            >  1:354106/1‑138 (MQ=255)                                                                                                                                  aaaCGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATccgcc                  >  1:15858/1‑148 (MQ=255)                                                                                                                                   aaCGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATccgccg                 <  4:384973/148‑1 (MQ=255)                                                                                                                                       tAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCt             <  2:213274/148‑1 (MQ=255)                                                                                                                                           tCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGa          >  4:102305/1‑147 (MQ=255)                                                                                                                                           tCGCCACTTCGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGa          >  1:343382/1‑147 (MQ=255)                                                                                                                                             gCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATa        <  5:99608/147‑1 (MQ=255)                                                                                                                                              ccAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGAt                                     <  4:277918/117‑1 (MQ=255)                                                                                                                                              ccAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGAt                                     >  3:277918/1‑117 (MQ=255)                                                                                                                                              ccAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAAt      <  3:390122/148‑1 (MQ=255)                                                                                                                                                  ttGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCAAGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGcg  >  2:326915/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCG  >  NC_000913/2117756‑2118048

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