Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,903,571	A→G	65.4%	T164A (ACC→GCC)	manY →	mannose‑specific enzyme IIC component of PTS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,903,571	0	A	G	65.4%	106.6 / 86.2	96	T164A (ACC→GCC)	manY	mannose‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base A (17/16);  new base G (29/33);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

GGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGG  >  NC_000913/1903472‑1903670                                                                                                    |                                                                                                    ggCAACCTGACAGTGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTa                                                                                                     >  2:1422468/1‑100 (MQ=255) ggCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTa                                                                                                     <  1:684854/100‑1 (MQ=255)  gCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTAc                                                                                                    >  1:2316786/1‑100 (MQ=255)      ccTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGc                                                                                                >  2:2138411/1‑100 (MQ=255)      ccTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGc                                                                                                >  2:2405853/1‑100 (MQ=255)      ccTGACAGCGATTTCCTGGATCCACGTTGCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGc                                                                                                >  2:1945797/1‑100 (MQ=255)       cTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCg                                                                                               >  1:717848/1‑100 (MQ=255)        tGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGa                                                                                              <  1:675839/100‑1 (MQ=255)        tGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGa                                                                                              >  1:64824/1‑100 (MQ=255)          aCAGCGATTTCCTGGATCCACGTTTCTTCTTTGTTCCTGCAAGCATTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAg                                                                                            <  1:2079589/100‑1 (MQ=255)          aCAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAg                                                                                            <  1:404998/100‑1 (MQ=255)           cAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGt                                                                                           >  2:951249/1‑100 (MQ=255)           cAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGt                                                                                           <  1:788447/100‑1 (MQ=255)             gCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTAc                                                                                         >  2:2052180/1‑100 (MQ=255)               gATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAg                                                                                       >  2:2563687/1‑100 (MQ=255)               gATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAg                                                                                       >  1:334058/1‑100 (MQ=255)                 tttCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGaa                                                                                     <  2:2239109/100‑1 (MQ=255)                  ttCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAAc                                                                                    <  1:208462/100‑1 (MQ=255)                    ccTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACAt                                                                                  >  1:2273227/1‑100 (MQ=255)                    ccTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACAt                                                                                  >  2:1358659/1‑100 (MQ=255)                      tGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGc                                                                                <  2:1866840/100‑1 (MQ=255)                        gATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTg                                                                              >  1:1203440/1‑100 (MQ=255)                           ccACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAAt                                                                           <  1:677413/100‑1 (MQ=255)                           ccACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAAt                                                                           >  1:2604743/1‑100 (MQ=255)                              cGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCg                                                                        >  1:1391001/1‑100 (MQ=255)                                tttCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGAt                                                                      <  2:2320434/100‑1 (MQ=255)                                 ttcttcTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGAt                                                                      >  2:946510/1‑99 (MQ=255)                                   cttctCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTcc                                                                   >  2:931751/1‑100 (MQ=255)                                   cttctCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTcc                                                                   <  1:2481035/100‑1 (MQ=255)                                    ttctCTTTTCCTGCAATCTTTGCTTGTGGCTATTACGGTTGTCATCGTTGCACTGTCTGTTGGTGCCAGAGAGGTACAGAACATGCTGAATGCGATTCCg                                                                  <  1:970941/100‑1 (MQ=255)                                    ttctCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCg                                                                  <  2:2011119/100‑1 (MQ=255)                                     tctcTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCgg                                                                 <  2:2595400/100‑1 (MQ=255)                                       tctGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGACCATGCTGAATGCGATTCCGGaa                                                               >  2:1226119/1‑100 (MQ=255)                                       tctGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGaa                                                               >  1:419292/1‑100 (MQ=255)                                       tctGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGaa                                                               >  1:2648994/1‑100 (MQ=255)                                        ctGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAg                                                              <  2:2648962/100‑1 (MQ=255)                                            tCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggt                                                          >  1:1343269/1‑100 (MQ=255)                                            tCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggt                                                          >  2:2604105/1‑100 (MQ=255)                                             ccTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggtg                                                         >  2:1496993/1‑100 (MQ=255)                                             ccTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggtg                                                         >  2:2140327/1‑100 (MQ=255)                                             ccTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggtg                                                         <  2:723/100‑1 (MQ=255)                                              cTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGa                                                        >  2:2161079/1‑100 (MQ=255)                                               tGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAc                                                       <  2:2493222/100‑1 (MQ=255)                                               tGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACATAACATGCTGAATGCGATTCCGGAAGTGGTGa                                                        <  2:2332705/99‑1 (MQ=255)                                               tGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAc                                                       <  1:1730380/100‑1 (MQ=255)                                                gcaagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAcc                                                      >  2:767339/1‑100 (MQ=255)                                                gcaagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAcc                                                      >  1:1926671/1‑100 (MQ=255)                                                gcaagcaaTGAGTGTGGCTATTCCGGCTGTAATCGTTGCGCTGTCTGATGGTGCCAGCGAAGTACACCACATGCTTAATGCGATTCTGGAAGTGGTGatc                                                      >  2:47168/1‑98 (MQ=255)                                                 caggcaaTGCGTGTGGCTATTCCGGACGTCATCGTTGCGATGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGTTTCTGGAAGTGGTTACCa                                                     >  2:2381730/4‑100 (MQ=255)                                                 caagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCa                                                     >  1:288600/1‑100 (MQ=255)                                                    gcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTCTCTGTTAGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATg                                                  >  1:2284169/1‑100 (MQ=255)                                                      aaTGCGTGTTGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATTCGATTCCGTAAGTGGTGACCAATGGt                                                <  2:2638394/100‑1 (MQ=255)                                                       aTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTc                                               <  1:1398948/100‑1 (MQ=255)                                                       aTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTc                                               <  1:2224154/100‑1 (MQ=255)                                                       aTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTc                                               <  1:2138437/100‑1 (MQ=255)                                                        tGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTc                                               <  2:1034299/99‑1 (MQ=255)                                                        tGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGtt                                               <  2:2549325/99‑2 (MQ=255)                                                         gCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTg                                             <  2:1993718/100‑1 (MQ=255)                                                          cGTGTGGCTATTCCGGCCGTTATTGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGa                                            <  1:730097/100‑1 (MQ=255)                                                          cGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGa                                            <  2:64824/100‑1 (MQ=255)                                                           gtgtGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGaa                                           <  1:1072231/100‑1 (MQ=255)                                                              tgGCTATTCCGGCAGTCATCGTTGCGCTGTCTGTTGGTACAAGCGAAGTACAGAACATGCTGAATGCGATTACGGAAGAGGTGACCAATGGTCTGAata                                         >  2:718998/1‑99 (MQ=255)                                                               ggCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATc                                       >  2:2577733/1‑100 (MQ=255)                                                                gCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCg                                      <  2:674051/100‑1 (MQ=255)                                                                   aTTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCg                                   >  2:1101352/1‑100 (MQ=255)                                                                     tCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGt                                 >  1:1756266/1‑100 (MQ=255)                                                                       cGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTgg                               >  1:369892/1‑100 (MQ=255)                                                                        ggCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGc                              <  2:368420/100‑1 (MQ=255)                                                                         gCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCa                             <  2:2054011/100‑1 (MQ=255)                                                                          ccGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCAt                            <  1:1487678/100‑1 (MQ=255)                                                                          ccGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCAt                            >  2:1073876/1‑100 (MQ=255)                                                                          ccGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCAt                            <  1:1922423/100‑1 (MQ=255)                                                                           cGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCCTg                           <  2:369890/100‑1 (MQ=255)                                                                           cGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCacg                           >  1:1305698/1‑98 (MQ=255)                                                                               aTCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCg                       >  2:772708/1‑100 (MQ=255)                                                                                tCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgt                      >  1:2633469/1‑100 (MQ=255)                                                                                 cGTTGCGCTGTCTGTTGGTTTCTTCGAAGTAAAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtg                     <  1:2065511/100‑1 (MQ=255)                                                                                   ttGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtggt                   <  1:931760/100‑1 (MQ=255)                                                                                   ttGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtggt                   <  2:445677/100‑1 (MQ=255)                                                                                    tGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtggtg                  <  1:1359857/100‑1 (MQ=255)                                                                                      cgcTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTggtggt                >  1:542012/1‑100 (MQ=255)                                                                                       gcTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggtt               <  1:1536659/100‑1 (MQ=255)                                                                                       gcTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATAGCCGGTAGCATGAACGTGGtggtt               >  1:576851/1‑100 (MQ=255)                                                                                         tgtctgtTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttgg             <  2:2273590/100‑1 (MQ=255)                                                                                          gtctgtTGGTGCCAGCGAAGTACAGAAGATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttggt            <  2:2284141/100‑1 (MQ=255)                                                                                          gtctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGTTGGCATGATCGTGGtggttggt            <  2:1760953/100‑1 (MQ=255)                                                                                          gtctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttggt            >  1:298376/1‑100 (MQ=255)                                                                                          gtctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttggt            <  2:2273199/100‑1 (MQ=255)                                                                                          gtctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTATGAATATCGCCGGTGGCATGATCGTGGtggttggt            >  1:1064465/1‑100 (MQ=255)                                                                                          gtctgtTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttggt            <  1:1422482/100‑1 (MQ=255)                                                                                          gtctgtTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttggt            >  1:2329890/1‑100 (MQ=255)                                                                                           tctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTggttggtt           <  2:717843/100‑1 (MQ=255)                                                                                           tctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATACGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTggttggtt           <  1:1196772/100‑1 (MQ=255)                                                                                           tctgtTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCTGTAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttggt            <  1:719003/99‑1 (MQ=255)                                                                                            ctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTa          <  2:281732/100‑1 (MQ=255)                                                                                            ctgtTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTa          <  1:442720/100‑1 (MQ=255)                                                                                              gtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATg        <  1:909293/100‑1 (MQ=255)                                                                                              gtTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCAAGATCGTGGTGGTTGGTTATg        >  1:453035/1‑100 (MQ=255)                                                                                                    aCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATgg  <  1:1073886/100‑1 (MQ=255)                                                                                                    |                                                                                                    GGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGG  >  NC_000913/1903472‑1903670

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