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	98.7%	254.9 / NA	75	D202A (GAT→GCT)	uxaB	tagaturonate reductase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base G (37/36);  minor base A (0/1);  total (38/37)
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.94e-01

GATTGCCACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGG  >  NC_000913/1609829‑1610316                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 gATTGCCACTTTGCGTTCGTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAGAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCa                                                                                                                                                                                                                                               <  2:578866/251‑1 (MQ=255)       cACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCAtt                                                                                                                                                                                                                                         >  1:453279/1‑251 (MQ=255)            ttcgttcgttGTACGGTTTAATATCGTCAACAAGCAGCACGTTGAGCGGATATTTGTCCAGACGTAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAAt                                                                                                                                                                                                                                     <  1:333472/248‑1 (MQ=255)                ttCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAAAAGGTAAAAGCGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGACGCTGTTAGCCTGAGCCAGCCATTGAATAAATg                                                                                                                                                                                                                                >  1:553383/1‑251 (MQ=255)                    ttGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGGCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGc                                                                                                                                                                                                                               <  2:354255/248‑1 (MQ=255)                     tGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTc                                                                                                                                                                                                                            <  1:639857/250‑1 (MQ=255)                     tGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCt                                                                                                                                                                                                                           >  2:244304/1‑251 (MQ=255)                          ggTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGc                                                                                                                                                                                                                        >  1:861351/1‑249 (MQ=255)                                atatCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGccc                                                                                                                                                                                                                 >  2:571379/1‑250 (MQ=255)                                 tatCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGGCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGACGTGCCCAc                                                                                                                                                                                                               >  2:557522/1‑251 (MQ=255)                                  atCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCAc                                                                                                                                                                                                               <  1:627978/250‑1 (MQ=255)                                  atCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCAAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGCTACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTCAGCCAGCCATTGAATAAATGCTTCTGGCAGTGccccct                                                                                                                                                                                                              >  1:379083/1‑248 (MQ=255)                                   tCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCActc                                                                                                                                                                                                             >  2:9105/1‑251 (MQ=255)                                    cGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCActct                                                                                                                                                                                                            <  2:329022/251‑1 (MQ=255)                                     gTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCtt                                                                                                                                                                                                           >  1:172585/1‑251 (MQ=255)                                      tCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCttc                                                                                                                                                                                                          >  1:693191/1‑250 (MQ=255)                                      tCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCACTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGCACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGATTCTGGCAGTGCCCACTCTTg                                                                                                                                                                                                          >  1:294644/1‑251 (MQ=255)                                       caacaaTCAGCACGTTGAACGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAAACGGTAACGATACGGTCCACCAGCGTAGGACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCAACTCTTgc                                                                                                                                                                                                         >  1:206159/1‑251 (MQ=255)                                        aacaaTCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCTGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCATGAGCCAGCCATTGAATAAAAGCTTCTGGCAGTGCCCACTCTGgcg                                                                                                                                                                                                        >  1:778976/1‑251 (MQ=255)                                                gCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTACACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAgcgc                                                                                                                                                                                                >  2:857560/1‑251 (MQ=255)                                                  ccGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGTTCCCTGAATCACAATCAGGGAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTGCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCgcag                                                                                                                                                                                              <  2:290576/250‑1 (MQ=255)                                                          cGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAACCAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGGATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGtt                                                                                                                                                                                      >  1:837574/1‑251 (MQ=255)                                                             atatTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCAc                                                                                                                                                                                   >  1:178003/1‑251 (MQ=255)                                                              tatTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGACACCTGAATCACTCCCAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTACCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCAACTCTTGCGCCTAGCGCCACACTAGTTCACg                                                                                                                                                                                  >  2:406880/1‑251 (MQ=255)                                                                 ttGTACAGACGCAGTTCAGTCGCTAAGGATGTCGGTCCCTGAATCACAAACAGGTAAAAGGGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTGCTTCCAGTTTAGCCACTTCAACGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCa                                                                                                                                                                                <  2:100417/250‑1 (MQ=255)                                                                     ccAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATAGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCCTAGCGCAGCACCAGTTCACGCAAcgcg                                                                                                                                                                           >  2:585419/1‑251 (MQ=255)                                                                        gACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCg                                                                                                                                                                        >  1:322278/1‑251 (MQ=255)                                                                           gCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGGCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCTCCAGTTCACGCAACGCGTCGCCa                                                                                                                                                                     <  1:797548/251‑1 (MQ=255)                                                                             aGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCAtt                                                                                                                                                                   <  2:306208/251‑1 (MQ=255)                                                                                     cGCTAAGGATTTCGGTCCCTTAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGtcaa                                                                                                                                                           <  2:564094/251‑1 (MQ=255)                                                                                     cGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGGTCAGCGGTGTCGATAAAACAATCGTGATGACCAAGGTCTTCTTCCAGTGTAGCCCCTACATCGCGCGTATAACCTGTAACAATACGCTCCACCAGCGTAGAACAGAAGCTGTTAACCTGAGCCAGCCATTGAATAAATGCTACTGGCAGTGCCGACTCTTGCGCCTAGCGCAGGACCAGATAACGCAACTCGTCGCCATTGTAGtaa                                                                                                                                                            >  1:33983/1‑248 (MQ=255)                                                                                       cTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCAtt                                                                                                                                                                   <  1:390939/241‑1 (MQ=255)                                                                                       cTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGCGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATCCGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCAtt                                                                                                                                                                   >  2:390939/1‑241 (MQ=255)                                                                                        tAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTgcg                                                                                                                                                                                                        <  1:872631/203‑1 (MQ=255)                                                                                        tAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTgcg                                                                                                                                                                                                        >  2:872631/1‑203 (MQ=255)                                                                                         aaGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGtcaatca                                                                                                                                                        >  2:53779/1‑250 (MQ=255)                                                                                         aaGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGtcaatc                                                                                                                                                         <  1:303268/249‑1 (MQ=255)                                                                                           ggATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACt                                                                                                                                                     >  1:215546/1‑251 (MQ=255)                                                                                           ggATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCGTCTGGCAGTGCCCACTCTTGCGCATAGAGCAGCACCAGTTCACGCAACGCGACGCCATTGTAGTCAATCAAcc                                                                                                                                                     >  2:410784/1‑250 (MQ=255)                                                                                                 cGGTCCCGGAAACACAAACTGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCAGCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCacac                                                                                                                                                <  2:299951/250‑1 (MQ=255)                                                                                                         gAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGGTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAatga                                                                                                                                         <  2:322278/249‑1 (MQ=255)                                                                                                             cacaAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATc                                                                                                                                    <  1:166314/250‑1 (MQ=255)                                                                                                                aaaCAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTAAATCAACTCACCCGGAATGATGATCCAAc                                                                                                                                >  2:534938/1‑251 (MQ=255)                                                                                                                      gTAAAAGTGTTCAGCGGTGTCGAGAAAACCATAGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTAt                                                                                                                          <  1:585419/251‑1 (MQ=255)                                                                                                                                  aggggTGTCGAGAAAACCATCGTGAAACCCCAGTTCTTCTTCCATGTTAGCCACTTCACCGCGCTGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCga                                                                                                              <  2:210277/248‑2 (MQ=255)                                                                                                                                  aGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGt                                                                                                              <  2:519191/251‑1 (MQ=255)                                                                                                                                            agaAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCATCGCACCGTTGAAATGGCt                                                                                                    >  2:150245/1‑251 (MQ=255)                                                                                                                                              aaaaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCAAACCTTTATCCAGCGCACCGTTGAAATGGCTa                                                                                                   >  2:183588/1‑250 (MQ=255)                                                                                                                                                      cGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCg                                                                                          <  1:560166/251‑1 (MQ=255)                                                                                                                                                      cGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTAGCCAGCGCACCGTTGAAATGGCTAAAGCGTTCg                                                                                          <  1:266705/251‑1 (MQ=255)                                                                                                                                                          aTAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAAca                                                                                      <  1:857560/251‑1 (MQ=255)                                                                                                                                                                cAGTTCTTCTTCAAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGAc                                                                                >  1:279357/1‑251 (MQ=255)                                                                                                                                                                       tcttcCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCa                                                                          <  1:819514/250‑1 (MQ=255)                                                                                                                                                                       tcttcCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGACATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCa                                                                          <  1:378234/250‑1 (MQ=255)                                                                                                                                                                          tcCAGTTGAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTTACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGtt                                                                       <  1:365784/250‑1 (MQ=255)                                                                                                                                                                             aGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGcc                                                                                                                                                                      >  1:398922/1‑152 (MQ=255)                                                                                                                                                                             aGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGcc                                                                                                                                                                      <  2:398922/152‑1 (MQ=255)                                                                                                                                                                                          cATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTAc                                                       <  1:557522/250‑1 (MQ=255)                                                                                                                                                                                           aTCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCg                                                     >  1:236163/1‑251 (MQ=255)                                                                                                                                                                                               cgcgGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTgg                                                 <  2:644510/251‑1 (MQ=255)                                                                                                                                                                                                gcgGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGgc                                                >  1:10636/1‑251 (MQ=255)                                                                                                                                                                                                  ggATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTaaa                                                                                                 <  1:152666/200‑1 (MQ=255)                                                                                                                                                                                                  ggATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTaaa                                                                                                 >  2:152666/1‑200 (MQ=255)                                                                                                                                                                                                        ccGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTaaa                                                                                                 >  1:488216/1‑194 (MQ=255)                                                                                                                                                                                                         cGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCg                                       <  2:10636/251‑1 (MQ=255)                                                                                                                                                                                                              aCGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAAttt                                  >  1:839110/1‑251 (MQ=255)                                                                                                                                                                                                                gATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTAt                                >  1:832897/1‑251 (MQ=255)                                                                                                                                                                                                                     tgTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGAAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGAACGCATACCGCTGGCGCGTCATCGAATTTACCGccc                           >  2:386634/2‑251 (MQ=255)                                                                                                                                                                                                                          accaGCGTAGAACAGAAGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATg                      <  1:355527/251‑1 (MQ=255)                                                                                                                                                                                                                                  aGAACTGATGCTGTTCGCCTCAGCCAGCCATTTAATAAATGCTTATGGAAGTGACCACTCTTGCGCATAGCGCAGCACCAGTGCACGCAACGCGGCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCACCCTTTATCCAGCGCTCCGTTGAAATGTCTAAAGTTGTCGAACAGCTGACGGTTCAGTTTGTCCGGATAGCTTACCGCTGGCGAGTCATCGATTTTATCGCCCTCATGGTAGCTGa              <  2:35404/251‑1 (MQ=255)                                                                                                                                                                                                                                  aGAACAGAAGCTGGTAGCCTGAGCCAGCCATTGAATAAGTACTTCTGGCGGTGCCCACGCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCCTTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCAACAGCACCGTTGTAATGGCTAAAGCGTTCGCGCAGCAGACGGGACAGCTGTGCCGGAAAGCTAACACCTGGCGCGGCATAGGATTTTTCCCCCGCATggaagttg               >  2:277072/1‑244 (MQ=255)                                                                                                                                                                                                                                          aGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCTCATGGTAGCTGATACCCGCt      <  2:713273/251‑1 (MQ=255)                                                                                                                                                                                                                                          aGCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCt      <  1:571379/251‑1 (MQ=255)                                                                                                                                                                                                                                           gCTGTTAGCCTGAGCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTATCTGATACCCgatt     <  1:246207/251‑4 (MQ=255)                                                                                                                                                                                                                                               atAGCCTGGACCAGCTATTGAATAAATGCTTCGGGCAGTTCCCCCCCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGGAGTCAATCAGCTCCCACGGAATGATGATCCAACCTTCATCCAGCGCCCCGTTTAAATGGCTAAAGCGTTCTAACAGCAGACGGGTCAGTTGTGCCGGATAGCGTACCGCTGGCGCGTCATCGAAGTTAACTCCCGCCTGGTAGCTGATACCCTCttgtg  <  2:837574/249‑4 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 GATTGCCACTTTGCGTTCTTTGTACGGTTTAATATCGTCAACAATCAGCACGTTGAGCGGATATTTGTCCAGACGCAGTTCAGTCGCTAAGGATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGCAGTGCCCACTCTTGCGCATAGCGCAGCACCAGTTCACGCAACGCGTCGCCATTGTAGTCAATCAACTCACACGGAATGATGATCCAACCTTTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGG  >  NC_000913/1609829‑1610316

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