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	92.5%	205.0 / 7.9	67	D202A (GAT→GCT)	uxaB	tagaturonate reductase
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base G (32/30);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-01

CGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTG  >  NC_000913/1609842‑1610321                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      cGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATaa                                                                                                                                                                                                                                        >  1:517258/1‑250 (MQ=255)          tACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTgg                                                                                                                                                                                                                              >  2:239508/1‑251 (MQ=255)                       cGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCActct                                                                                                                                                                                                                 <  2:411176/251‑1 (MQ=255)                        gTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCtt                                                                                                                                                                                                                <  1:536637/251‑1 (MQ=255)                          caacaaTCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTgc                                                                                                                                                                                                              >  2:379245/1‑251 (MQ=255)                           aacaaTCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACCGAAGCTGTTAGCCTGAGCCCGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTgcg                                                                                                                                                                                                             >  2:219032/1‑251 (MQ=255)                               aTCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCAt                                                                                                                                                                                                          >  2:391282/1‑250 (MQ=255)                                 cAGCAAGTTGAGCGGATCTTTGTCCAAACGCAGTTCAGACGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAACGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGAAACGGTCCCCCAGAGTAGAACAGAAGCTGTTAGCCTGAGCAAGCCATTGACTACATGCTTCTGGCAGTGCCCACTCTTGCGCATatc                                                                                                                                                                                                       >  2:269446/1‑249 (MQ=255)                                  aGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAgcg                                                                                                                                                                                                      <  2:517258/251‑1 (MQ=255)                                    cACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGACTAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAgcgc                                                                                                                                                                                                     >  2:279725/1‑250 (MQ=255)                                    cACGTTGAGCGGATATTTGGCCAGACGCAGTTCAGTCGCGAAGGATTTCGGCACCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCAGCGTGAGAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCgca                                                                                                                                                                                                    <  2:127450/251‑1 (MQ=255)                                        ttGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCAc                                                                                                                                                                                                >  1:473980/1‑251 (MQ=255)                                               gATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCa                                                                                                                                                                                         <  1:264576/251‑1 (MQ=255)                                                 aatTTGGCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGGTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACg                                                                                                                                                                                       <  2:106690/250‑1 (MQ=255)                                                        ccAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAAcgc                                                                                                                                                                                 <  2:298117/250‑1 (MQ=255)                                                        ccAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAAcgc                                                                                                                                                                                 <  1:53089/250‑1 (MQ=255)                                                          aGACGCAGTTCAGTCGCTAAGGATTTCGGTCCATGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTc                                                                                                                                                                              <  1:150695/251‑1 (MQ=255)                                                            aCGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGc                                                                                                                                                                            >  1:460827/1‑251 (MQ=255)                                                                      gTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGtc                                                                                                                                                                  <  1:76967/251‑1 (MQ=255)                                                                        cGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGtca                                                                                                                                                                 <  1:473456/250‑1 (MQ=255)                                                                         gCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGtcaat                                                                                                                                                               <  1:125340/251‑1 (MQ=255)                                                                                 tttCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTcac                                                                                                                                                       >  1:9371/1‑251 (MQ=255)                                                                                  ttCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTcac                                                                                                                                                       <  1:12279/250‑1 (MQ=255)                                                                                      gTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACgg                                                                                                                                                   >  1:137335/1‑250 (MQ=255)                                                                                      gTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGAGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGACTAAATGCTTCTGGCAGGGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCCACGCGTCGCCATTGTAGTCAATCAACTcaaccgg                                                                                                                                                   >  2:43539/1‑245 (MQ=255)                                                                                        cccTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCgca                                                                                                                                                                                                    <  2:33014/199‑1 (MQ=255)                                                                                        cccTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCgca                                                                                                                                                                                                    >  1:33014/1‑199 (MQ=255)                                                                                            gAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATACCCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAatgat                                                                                                                                             <  2:194506/250‑1 (MQ=255)                                                                                            gAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAatgatg                                                                                                                                            <  1:178504/251‑1 (MQ=255)                                                                                                   aaaCAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAAc                                                                                                                                     >  1:302120/1‑251 (MQ=255)                                                                                                           aaaaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATcc                                                                                                                             >  1:445795/1‑251 (MQ=255)                                                                                                                 gtTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAgcgc                                                                                                                        >  1:530151/1‑250 (MQ=255)                                                                                                                     aGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGCATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCACCCTTTATCCAGCGCACCGt                                                                                                                   >  2:496197/1‑251 (MQ=255)                                                                                                                          tgtCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGaa                                                                                                               >  1:409977/1‑250 (MQ=255)                                                                                                                           gtCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCCGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGCAGCCCAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTAGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGCAAt                                                                                                             >  2:480999/1‑251 (MQ=255)                                                                                                                               agaAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACAAGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAAtagc                                                                                                          >  2:539418/1‑247 (MQ=255)                                                                                                                                               aaCCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAAcag                                                                                          <  2:409977/250‑1 (MQ=255)                                                                                                                                                       tcttcttcCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCTCTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCGGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACggg                                                                                  <  1:219032/250‑1 (MQ=255)                                                                                                                                                             tcCATTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTGTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGtt                                                                            <  2:525364/250‑1 (MQ=255)                                                                                                                                                             tcCAGTTTAGCCACGTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGt                                                                             <  1:74411/249‑1 (MQ=255)                                                                                                                                                                aGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGc                                                                        >  1:527130/1‑251 (MQ=255)                                                                                                                                                                 gTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGcc                                                                       <  1:279725/251‑1 (MQ=255)                                                                                                                                                                 gTTTAGCCACTTAATCGCGCGGATAACCGGTAACGATACGGACCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGACGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGCACAGCAGACGGGTCAGTTTTGcc                                                                       >  1:125785/1‑251 (MQ=255)                                                                                                                                                                        cACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTTCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGc                                                                >  1:247630/1‑251 (MQ=255)                                                                                                                                                                           ttCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGCTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCCGTTTTGCCGGATAGCtt                                                              >  1:37868/1‑250 (MQ=255)                                                                                                                                                                            tCATCGCGCTGATAAACGGTAACGTTACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTTAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAAGCAACGCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGTCTAAAGTGTTCGAACAGCAGACGGGGCAGTTTTGCCGGATAgttt                                                              <  2:32739/249‑4 (MQ=255)                                                                                                                                                                                 gcgcgGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTg                                                       <  1:548669/251‑1 (MQ=255)                                                                                                                                                                                  cgcgGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTgg                                                      >  1:174550/1‑251 (MQ=255)                                                                                                                                                                                     ggATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGCAGCTGTTCGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCAAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGCCGGGTCAGTTTTGCCGGCTAGCTTAACGCTTgcgc                                                   >  1:302111/1‑251 (MQ=255)                                                                                                                                                                                         aaCCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCa                                               <  2:191043/251‑1 (MQ=255)                                                                                                                                                                                              gcaaCGATACGGGCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGaa                                          <  2:392137/249‑1 (MQ=255)                                                                                                                                                                                                aaCGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGCAtt                                        >  1:68654/1‑251 (MQ=255)                                                                                                                                                                                                aaCGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAAtt                                        >  1:420020/1‑251 (MQ=255)                                                                                                                                                                                                 aCGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAAttt                                       <  2:467253/251‑1 (MQ=255)                                                                                                                                                                                                    aTACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATc                                    >  2:117110/1‑251 (MQ=255)                                                                                                                                                                                                        ggTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCTATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGccc                                >  1:194012/1‑251 (MQ=255)                                                                                                                                                                                                         gTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCg                               <  1:28499/251‑1 (MQ=255)                                                                                                                                                                                                            caccaGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTAACGCTGGCGCGTCATCGAATTTATCGCCCGCAt                            >  1:243420/1‑251 (MQ=255)                                                                                                                                                                                                             accaGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATg                           >  2:294536/1‑251 (MQ=255)                                                                                                                                                                                                                 gCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCTGATATCTTTCCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAg                       <  1:362914/251‑1 (MQ=255)                                                                                                                                                                                                                          aGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATAccc              <  2:527130/251‑1 (MQ=255)                                                                                                                                                                                                                               cTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTc         >  1:336/1‑251 (MQ=255)                                                                                                                                                                                                                                tGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCg        >  2:195962/1‑251 (MQ=255)                                                                                                                                                                                                                                 tttAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTc         <  2:445795/248‑1 (MQ=255)                                                                                                                                                                                                                                     gCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCggtggt   <  1:294536/251‑1 (MQ=255)                                                                                                                                                                                                                                      ccTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCCACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGgtggtg  >  2:371909/1‑251 (MQ=255)                                                                                                                                                                                                                                      ccTGAGCCAGCCATTGAATAAACGCCGCTGGCAGTGCCCACGCTTGCGCATAGCGCAGCACTAGGTCACGCAACGCGTCGCCATTGTAGTCAATCAACCCCCACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGgcggtg  <  2:471814/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      CGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTG  >  NC_000913/1609842‑1610321

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