Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,903,571	A→G	100%	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	94.1%	291.4 / 3.2	95	T164A (ACC→GCC)	manY	mannose‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (48/41);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGA  >  NC_000913/1903474‑1903667                                                                                                  |                                                                                                 cAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGcc                                                                                                <  1:1944314/100‑1 (MQ=255) cAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGcc                                                                                                >  2:592184/1‑100 (MQ=255) cAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGcc                                                                                                <  1:592181/100‑1 (MQ=255)      tGACAGTGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGa                                                                                           <  2:2223762/100‑1 (MQ=255)       gACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGaa                                                                                          <  1:312539/100‑1 (MQ=255)         cAGCGATTTCATTGATCCTCGTGTCTTCTTTGTTCCTGCAAGCAATGTGTGTGGATATTCCGGCCGTCATCTTTGCGCTTTCTGTTGGTGCCAGCGAAGt                                                                                        <  1:644135/100‑1 (MQ=255)          aGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTa                                                                                       <  1:2411062/100‑1 (MQ=255)           gCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTAc                                                                                      >  2:1589512/1‑100 (MQ=255)           gCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTAc                                                                                      >  2:1451695/1‑100 (MQ=255)            tgATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACa                                                                                     <  1:468366/99‑1 (MQ=255)             gATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACaa                                                                                    >  2:468368/1‑99 (MQ=255)              aTTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGa                                                                                   >  1:1899772/1‑100 (MQ=255)                ttCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAAc                                                                                 <  1:424246/100‑1 (MQ=255)                ttCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAAc                                                                                 >  1:1368162/1‑100 (MQ=255)                  ccTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACAt                                                                               >  2:2273345/1‑100 (MQ=255)                  ccTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCAGCCGTCATCATTGCGGTGTCTGTTGGTGCCAGCGAAGTACAGAACAt                                                                               >  1:1750582/1‑100 (MQ=255)                    tGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGc                                                                             >  1:2531592/1‑100 (MQ=255)                     ggATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCt                                                                            <  1:690366/100‑1 (MQ=255)                      gATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTg                                                                           <  1:2041111/100‑1 (MQ=255)                        tCCTCGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGa                                                                          <  1:2343972/99‑1 (MQ=255)                        tCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGaa                                                                         <  1:2363644/100‑1 (MQ=255)                         ccACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGTCCGTCATCGTTGCGCGGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAAt                                                                        >  1:1156428/1‑100 (MQ=255)                         ccACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAAt                                                                        >  2:2565348/1‑100 (MQ=255)                         ccACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAAt                                                                        >  1:269102/1‑100 (MQ=255)                          cACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGTTGAATg                                                                       <  2:2374096/100‑1 (MQ=255)                          cACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATg                                                                       >  1:2508704/1‑100 (MQ=255)                           aCGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGc                                                                      <  1:24800/100‑1 (MQ=255)                           aCGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGc                                                                      <  2:1634190/100‑1 (MQ=255)                           aCGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGc                                                                      <  2:1376603/100‑1 (MQ=255)                             gTTTCTTCTTTCTTCCTGCAAGCCATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTGGGTGCCTGCGAAGTACAGAACATGCTGAATGCGa                                                                    >  1:693607/1‑100 (MQ=255)                             gTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGa                                                                    >  1:1722121/1‑100 (MQ=255)                                 cttctCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTcc                                                                >  2:714267/1‑100 (MQ=255)                                 ccactcgattCCGGCAAGCAATTCGGGTTGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTcc                                                                <  2:476717/92‑1 (MQ=255)                                   tctcTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCgg                                                              >  2:327223/1‑100 (MQ=255)                                    ctctGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGa                                                             <  2:259305/100‑1 (MQ=255)                                    ctctGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGa                                                             >  1:200297/1‑100 (MQ=255)                                    ctctGTTCCTCCAAGCAATGCGTGTGGCTATTCCGTCCGTCATCGTTGCGCTGTTGGTTGGTGCCAGCGAAGTAGAGTACATTCTGAATGCGATTACGGa                                                             <  1:74314/100‑1 (MQ=255)                                      ctGCTCCTGCAAGATATACGTGTCGATATTACCGTCATCATCGTTGCGCTGTCTGCTGGTGCCAGCGAATTACAGAACATGCTGAATGCGAGTCCg                                                               <  1:1507728/96‑1 (MQ=255)                                       tGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgt                                                          >  2:2234463/1‑100 (MQ=255)                                       tGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgt                                                          >  1:2593700/1‑100 (MQ=255)                                       tGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgt                                                          >  1:1209104/1‑100 (MQ=255)                                        gTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtg                                                         >  2:369355/1‑100 (MQ=255)                                         ttCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtgg                                                        <  2:1899782/100‑1 (MQ=255)                                         tctcTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtgg                                                        <  1:1079890/97‑1 (MQ=255)                                          tCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggt                                                       >  2:2142489/1‑100 (MQ=255)                                          tCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggt                                                       >  2:439583/1‑100 (MQ=255)                                            cTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggaga                                                     >  2:1079894/1‑97 (MQ=255)                                            cTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGa                                                     >  1:1828833/1‑100 (MQ=255)                                             tGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAc                                                    <  2:200299/100‑1 (MQ=255)                                             tGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAc                                                    <  2:582336/100‑1 (MQ=255)                                             tGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAc                                                    >  1:2258797/1‑100 (MQ=255)                                              gcaagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCGGGAAGTGGTGAcc                                                   <  1:513474/100‑1 (MQ=255)                                              gcaagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAcc                                                   >  2:722428/1‑100 (MQ=255)                                               caagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCa                                                  >  2:55131/1‑100 (MQ=255)                                                aagcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCaa                                                 >  1:1867989/1‑100 (MQ=255)                                                 agcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAAt                                                <  1:2181086/100‑1 (MQ=255)                                                 agcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAAt                                                >  2:237656/1‑100 (MQ=255)                                                  gcaaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATg                                               >  2:360580/1‑100 (MQ=255)                                                   caaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATgg                                              <  2:1624708/100‑1 (MQ=255)                                                   caaTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATgg                                              >  2:1774594/1‑100 (MQ=255)                                                     aTGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTc                                            >  2:1890427/1‑100 (MQ=255)                                                         gtgtGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGaa                                        <  1:1251318/100‑1 (MQ=255)                                                         gtgtGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGaa                                        <  1:36051/100‑1 (MQ=255)                                                          tgtgGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAat                                       >  1:578362/1‑100 (MQ=255)                                                           gtgGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCAAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAata                                      >  1:355419/1‑100 (MQ=255)                                                             ggCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATc                                    <  2:202883/100‑1 (MQ=255)                                                                tATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGACCATGCTGAATGCGACTCCGGAAGTGGTGACCAATGGTCTGAATATCGcc                                 >  2:592413/1‑100 (MQ=255)                                                                 aTTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCg                                >  1:111565/1‑100 (MQ=255)                                                                  ttCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTATGAATATCGCCgg                               <  1:432415/100‑1 (MQ=255)                                                                   tCCGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGt                              >  2:2305932/1‑100 (MQ=255)                                                                    ccGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTg                             >  2:559045/1‑100 (MQ=255)                                                                    ccGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTg                             >  2:945205/1‑100 (MQ=255)                                                                    ccGGCCGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTg                             >  1:2276817/1‑100 (MQ=255)                                                                     cGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTgg                            <  1:592410/100‑1 (MQ=255)                                                                        ccGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCAt                         <  2:899879/100‑1 (MQ=255)                                                                         cGTCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATg                        <  2:1757147/100‑1 (MQ=255)                                                                           tCATCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGAt                      <  1:762333/100‑1 (MQ=255)                                                                           tCATCGTTGCGCTGTCTGTTGGTGCCAGCAAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGCGACCAATGGTCTGAATATCGCCGGTGGCATGAt                      <  1:1552196/100‑1 (MQ=255)                                                                            naTCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAACAATGGTCTGAATATCGCCGGTGGCATGATc                     >  1:1266490/2‑100 (MQ=255)                                                                             aTCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCg                    <  2:1461980/100‑1 (MQ=255)                                                                              tCGTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgt                   <  2:1114393/100‑1 (MQ=255)                                                                                gTTGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtgg                 >  1:307992/1‑100 (MQ=255)                                                                                 ttGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtggt                <  1:1466869/100‑1 (MQ=255)                                                                                  tGCGCTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCgtggtg               >  2:1425707/1‑100 (MQ=255)                                                                                    cgcTGTCTGTTGCCGCCAGCGAAGTACAGAACATGCTGACTGCTATTACGGAAGTGGTGACCAATGGTATGAATATCGCCGGTGGCATGATc                     >  1:1375702/1‑92 (MQ=255)                                                                                     gcTGTCTGTTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggtt            <  2:2351869/100‑1 (MQ=255)                                                                                       tgtctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttgg          >  2:881344/1‑100 (MQ=255)                                                                                        gtctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttggt         <  2:418697/100‑1 (MQ=255)                                                                                        gtctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGtggttggt         >  2:1992954/1‑100 (MQ=255)                                                                                          ctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTa       >  2:280622/1‑100 (MQ=255)                                                                                          ctgtTGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTa       >  1:147036/1‑100 (MQ=255)                                                                                             ttGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGc    <  2:2178042/100‑1 (MQ=255)                                                                                              tGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATgtg   <  2:1754456/100‑3 (MQ=255)                                                                                              tGGTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCg   <  2:307994/100‑1 (MQ=255)                                                                                               ggTGCCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGa  >  1:48141/1‑100 (MQ=255)                                                                                                  |                                                                                                 CAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGA  >  NC_000913/1903474‑1903667

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: ‑