Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,271,629	T→C	T48A (ACG→GCG)	garK ←	glycerate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,271,629	0	T	C	99.0%	378.0 / ‑4.4	100	T48A (ACG→GCG)	garK	glycerate kinase I
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (51/48);  total (51/49)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

CCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTA  >  NC_000913/3271482‑3271775                                                                                                                                                    |                                                                                                                                                   ccAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTcc                                                                                                                                                 >  1:494256/1‑151 (MQ=255) ccAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTcc                                                                                                                                                 <  2:388058/151‑1 (MQ=255)   aGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCtt                                                                                                                                               >  2:877085/1‑151 (MQ=255)         cTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCa                                                                                                                                          <  2:150928/150‑1 (MQ=255)          tGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCAcc                                                                                                                                        >  1:704336/1‑151 (MQ=255)              cgccgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCg                                                                                                                                    <  1:1010722/151‑1 (MQ=255)                 cgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGc                                                                                                                                  >  1:1378984/1‑150 (MQ=255)                   ccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAAc                                                                                                                               >  2:825656/1‑151 (MQ=255)                     aTTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCg                                                                                                                             >  1:843605/1‑151 (MQ=255)                     aTTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCg                                                                                                                             >  2:460110/1‑151 (MQ=255)                       ttCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCgg                                                                                                                            <  2:483637/150‑1 (MQ=255)                           aTAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAAc                                                                                                                         >  2:496203/1‑149 (MQ=255)                            tAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAAc                                                                                                                         >  1:563089/1‑148 (MQ=255)                             aaaCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGa                                                                                                                      <  2:1116369/150‑1 (MQ=255)                              aaCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGaa                                                                                                                     >  1:349453/1‑150 (MQ=255)                               aCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAAc                                                                                                                   <  2:851614/151‑1 (MQ=255)                               aCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAAc                                                                                                                   <  2:967392/151‑1 (MQ=255)                                cgcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACg                                                                                                                  <  2:573830/151‑1 (MQ=255)                                cgcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACg                                                                                                                  <  2:685237/151‑1 (MQ=255)                                 gcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGt                                                                                                                 <  2:217793/151‑1 (MQ=255)                                    gTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACt                                                                                                              >  1:846528/1‑151 (MQ=255)                                      tttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGGAACGTACt                                                                                                              <  1:865712/149‑1 (MQ=255)                                       ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACtgtg                                                                                                           >  1:336148/1‑151 (MQ=255)                                        tGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGc                                                                                                          >  2:1230920/1‑151 (MQ=255)                                         gCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGc                                                                                                          >  2:465010/1‑150 (MQ=255)                                         gCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCa                                                                                                         <  1:411462/151‑1 (MQ=255)                                          ccATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCAt                                                                                                        <  2:913737/151‑1 (MQ=255)                                           cATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATc                                                                                                       >  1:1072110/1‑151 (MQ=255)                                           cATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACACAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATc                                                                                                       >  1:721907/1‑151 (MQ=255)                                             tCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATc                                                                                                       >  1:4976/1‑149 (MQ=255)                                              cGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAg                                                                                                     <  2:1329003/150‑1 (MQ=255)                                                ccGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGaa                                                                                                  <  1:775738/151‑1 (MQ=255)                                                 cGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGaaa                                                                                                 <  1:932135/151‑1 (MQ=255)                                                 cGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGaaa                                                                                                 <  1:358841/151‑1 (MQ=255)                                                 cGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGaa                                                                                                  >  1:1276675/1‑150 (MQ=255)                                                       tCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTc                                                                                           >  2:120500/1‑151 (MQ=255)                                                         ccccAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTccc                                                                                         >  2:109501/1‑151 (MQ=255)                                                           ccAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCg                                                                                        >  2:1194395/1‑150 (MQ=255)                                                            cAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGa                                                                                       >  2:1272260/1‑150 (MQ=255)                                                              aCTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAAt                                                                                    >  2:1213477/1‑151 (MQ=255)                                                                tGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATcc                                                                                  >  2:922909/1‑151 (MQ=255)                                                                 ggCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCt                                                                                 <  2:96360/151‑1 (MQ=255)                                                                   cATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCttt                                                                               >  2:1357983/1‑151 (MQ=255)                                                                   cATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCttt                                                                               <  2:1297129/151‑1 (MQ=255)                                                                   cATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCttt                                                                               <  2:1035524/151‑1 (MQ=255)                                                                    atccgcTTTCTCGCACCACGGCCCTGTAGCCCAGGCAGGACATTCAGCCCGCTGGGTGGCTACAATCATCGCTTCCTCCGCTCCTTCGCCACCGTCGGCACCCGGAACAGAAACGTAATGTGGATCACGAAAAATTTCCCGAAATCCtt                                                                                <  1:1328941/144‑1 (MQ=255)                                                                     ttCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCttttt                                                                             >  1:1142345/1‑151 (MQ=255)                                                                         cTTTCTCGCCCAGCGGCCCTGGAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCCCTCCCTCGCCACCGTCGGCAACCGGAACAGTAACGTACTGTGCATCAGGAAAGATTTCCCGAAATCATTTTTCTAt                                                                         >  1:743926/1‑151 (MQ=255)                                                                           ttCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATc                                                                        <  1:980793/150‑1 (MQ=255)                                                                            tctcGCCCAGCGGCCCTGTAACCGAGGCGTGACGTTCAGCCCCTTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGCGCATCAGGAAAAAGTTCCCGCAATCCTTTTTCTa                                                                          >  1:141714/1‑147 (MQ=255)                                                                               cGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCt                                                                    <  2:821505/150‑1 (MQ=255)                                                                                  ccAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTgcgc                                                                <  2:868657/151‑1 (MQ=255)                                                                                   cAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTgcgc                                                                >  2:344055/1‑150 (MQ=255)                                                                                   cAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCa                                                               <  2:514402/151‑1 (MQ=255)                                                                                    aGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCaa                                                              <  2:38218/151‑1 (MQ=255)                                                                                     gCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAAc                                                             <  1:1337114/151‑1 (MQ=255)                                                                                        gCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTc                                                          >  2:75262/1‑151 (MQ=255)                                                                                        gCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTc                                                          >  2:300549/1‑151 (MQ=255)                                                                                          ccTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGc                                                        >  1:1307379/1‑151 (MQ=255)                                                                                          ccTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGc                                                        >  2:180756/1‑151 (MQ=255)                                                                                            tGTAACCCAGGCGTGACGTTCAGCCACCTGGGAGGCTGAAATCATCGCATGCAACGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGGACTGTGCATGAGGATAAATTTCCCGAAATCCTTTTTCTAGCGCCTGCGCCACATCg                                                         >  1:1172050/1‑148 (MQ=255)                                                                                             gTACCCCCGGCGTGACGTTCAGCCCCCTGGGTGACTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGTAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCATCCTCGCt                                                       <  1:1319495/149‑1 (MQ=255)                                                                                             gTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTg                                                      <  1:749399/150‑1 (MQ=255)                                                                                               aaCCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGc                                                    >  1:778503/1‑150 (MQ=255)                                                                                               aaCCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGc                                                    <  2:267904/150‑1 (MQ=255)                                                                                                aCCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAg                                                  <  2:10164/151‑1 (MQ=255)                                                                                                 cccAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGc                                                    >  1:623476/1‑148 (MQ=255)                                                                                                  ccAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCACTCATCGCTTCCACCGCTCCTTCGCCACCTTCGGCAGCCGGAACAGAAACTTACTGTCCATAAGGAAAAATTTCTCGAAATCCTGTTGCTATCGCCTGCGCAACCACGCTGGCAGAt                                                <  1:361503/151‑1 (MQ=255)                                                                                                  ccAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGc                                                    >  2:1348192/1‑147 (MQ=255)                                                                                                    aGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATa                                               <  1:1292525/150‑1 (MQ=255)                                                                                                    aGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATa                                               >  2:267038/1‑150 (MQ=255)                                                                                                        gTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAActt                                          >  1:1278707/1‑151 (MQ=255)                                                                                                           aCGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAActttct                                       <  2:128796/151‑1 (MQ=255)                                                                                                           aCGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAActttct                                       <  1:1309279/151‑1 (MQ=255)                                                                                                            cGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAActttctt                                      >  1:758648/1‑151 (MQ=255)                                                                                                            cGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAActttctt                                      >  2:1342432/1‑151 (MQ=255)                                                                                                            cGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAActttctt                                      >  1:985362/1‑151 (MQ=255)                                                                                                                cAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTtata                                  <  1:1172136/151‑1 (MQ=255)                                                                                                                cAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTtat                                   >  1:558815/1‑150 (MQ=255)                                                                                                                 aGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATaa                                 <  1:744420/151‑1 (MQ=255)                                                                                                                    ccccTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCGGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAagag                              <  1:273204/151‑1 (MQ=255)                                                                                                                    ccccTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAagag                              <  2:273131/151‑1 (MQ=255)                                                                                                                      ccTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTc                            >  2:756442/1‑151 (MQ=255)                                                                                                                       cTGGGTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTc                            <  1:347889/150‑1 (MQ=255)                                                                                                                         gggTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGATCACTGCCGTCCTTTTCCTCTTTCCTTCTTTCCCTCTTTCCTTTTTCTCTCGCCTTCTCCACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCt                           <  2:974536/149‑1 (MQ=255)                                                                                                                         gggTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCCTCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCt                           >  1:1305867/1‑149 (MQ=255)                                                                                                                         gggTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCt                           <  1:825620/149‑1 (MQ=255)                                                                                                                         gggTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTg                          <  2:420055/150‑1 (MQ=255)                                                                                                                           gTGGCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTgggg                       >  2:188595/1‑151 (MQ=255)                                                                                                                              gCTGCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCgt                    >  1:337050/1‑150 (MQ=255)                                                                                                                                 gCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGAtt                  >  1:1127804/1‑150 (MQ=255)                                                                                                                                 gCAATCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTa                 <  2:1083719/151‑1 (MQ=255)                                                                                                                                   aaTCATCGCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACg               <  1:894726/151‑1 (MQ=255)                                                                                                                                          gCTTCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTc         >  2:38072/1‑150 (MQ=255)                                                                                                                                             tCCACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCAt       <  1:580813/149‑1 (MQ=255)                                                                                                                                               cACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCt   <  2:485808/151‑1 (MQ=255)                                                                                                                                               cACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCt   <  2:809083/151‑1 (MQ=255)                                                                                                                                               cACCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCt   >  1:1085866/1‑151 (MQ=255)                                                                                                                                                aCCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCt   >  1:1354321/1‑150 (MQ=255)                                                                                                                                                aCCGCTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTa  <  2:259328/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCTGGCAGATAAACTTTCTTTATAAGAGTCTGGGGCGATTACGATTTTCATACCTA  >  NC_000913/3271482‑3271775

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