Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	520,993	T→G	100%	V193G (GTC→GGC)	ybbP →	putative ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	520,993	0	T	G	100.0%	305.4 / ‑6.7	89	V193G (GTC→GGC)	ybbP	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (42/46);  minor base A (1/0);  total (43/46)

GGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGG  >  NC_000913/520898‑521089                                                                                                |                                                                                                 ggATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAcc                                                                                              >  2:1186044/1‑100 (MQ=255) ggATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAcc                                                                                              <  1:384760/100‑1 (MQ=255) gatctCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAcc                                                                                              <  1:1555315/96‑1 (MQ=255)  gATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCt                                                                                             <  2:991790/100‑1 (MQ=255)  gATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCt                                                                                             <  2:92884/100‑1 (MQ=255)   aTTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTg                                                                                            >  1:1331643/1‑100 (MQ=255)    ttCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTgg                                                                                           >  1:1121243/1‑100 (MQ=255)    ttCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTgg                                                                                           >  1:1809148/1‑100 (MQ=255)    ttCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTgg                                                                                           <  1:1571261/100‑1 (MQ=255)      ccGGTTTTAACCCCTTCCAGATGGCTCGGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         <  2:82199/100‑1 (MQ=255)      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAccagatc                                                                                         >  2:1555304/1‑95 (MQ=255)      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         <  1:91098/100‑1 (MQ=255)      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         >  1:631386/1‑100 (MQ=255)       cGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGt                                                                                        <  1:2266967/100‑1 (MQ=255)         gTTTTAACCCCTTCCAGATGGCTCCGCGTCGGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTta                                                                                      <  1:750985/100‑1 (MQ=255)           tttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTtata                                                                                    <  2:2246440/100‑1 (MQ=255)           tttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTtata                                                                                    <  2:2473802/100‑1 (MQ=255)           tttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTtata                                                                                    <  2:644853/100‑1 (MQ=255)           tttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGACGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGGCGGGGCACCTGGCGTtata                                                                                    >  2:302677/1‑100 (MQ=255)                 cccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcg                                                                              >  2:2259436/1‑100 (MQ=255)                 cccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcg                                                                              <  1:410707/100‑1 (MQ=255)                 cccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcg                                                                              <  2:376777/100‑1 (MQ=255)                   cTTCCAGATGGCTCCGCGTCTGATGATGAATCTGTCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTCATATATTcggc                                                                            <  2:1346647/100‑1 (MQ=255)                   cTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggc                                                                            <  1:1280479/100‑1 (MQ=255)                    ttCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcg                                                                           <  1:2006666/100‑1 (MQ=255)                     tCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcgg                                                                          <  1:1622059/100‑1 (MQ=255)                         gATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAAc                                                                      >  1:286074/1‑100 (MQ=255)                          aTGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACg                                                                     <  1:1055669/100‑1 (MQ=255)                           tGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACga                                                                    >  1:1627798/1‑100 (MQ=255)                             gCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACgaga                                                                  >  1:2365648/1‑100 (MQ=255)                             gCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACgaga                                                                  >  2:897329/1‑100 (MQ=255)                                  gcgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCa                                                              <  2:398293/99‑1 (MQ=255)                                    gTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCt                                                           >  2:2351289/1‑100 (MQ=255)                                     tCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTc                                                          >  1:624643/1‑100 (MQ=255)                                       tgatgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCg                                                         >  2:160414/1‑99 (MQ=255)                                        gatgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGAc                                                       >  2:1487661/1‑100 (MQ=255)                                         atgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACg                                                      <  1:972049/100‑1 (MQ=255)                                          tgatgaATCTGGCGGATGTCGATACAACCGGAGCCGTCCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGCCAACGAGAACAAGCTCGAc                                                       >  2:2031030/1‑98 (MQ=255)                                          tgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACgg                                                     >  1:1920727/1‑100 (MQ=255)                                          tgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACgg                                                     >  1:641422/1‑100 (MQ=255)                                          tgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACgg                                                     >  1:1916595/1‑100 (MQ=255)                                           gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGc                                                    >  1:2043179/1‑100 (MQ=255)                                           gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGc                                                    <  1:2353162/100‑1 (MQ=255)                                               aaTCTGGTGGATGTCGATAAATCCGGAACCTTGCAGCCGGGGAGTCGGGGCACCTGGCATTATAAATTCAGCAGGAACGAAAACCAGCTCGACGGCTATg                                                <  2:733237/100‑1 (MQ=255)                                                aTCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGTTCGACGGCTATga                                               <  1:1741554/100‑1 (MQ=255)                                                 tCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATgag                                              <  2:624636/100‑1 (MQ=255)                                                  cTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAGCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATgaga                                             >  1:1516262/1‑100 (MQ=255)                                                   tGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATgaga                                             >  2:436833/1‑99 (MQ=255)                                                    ggCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGaaa                                           >  2:283136/1‑100 (MQ=255)                                                    ggCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGaaa                                           >  1:652907/1‑100 (MQ=255)                                                      cGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATg                                         >  1:894/1‑100 (MQ=255)                                                       ggATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATgg                                        <  1:1410792/100‑1 (MQ=255)                                                       ggATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATgg                                        <  1:788629/100‑1 (MQ=255)                                                       ggATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATgg                                        <  1:1395845/100‑1 (MQ=255)                                                         aTGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgtt                                      <  2:2469924/100‑1 (MQ=255)                                                          tGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttg                                     >  2:971254/1‑100 (MQ=255)                                                          tGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttg                                     <  1:127170/100‑1 (MQ=255)                                                          tGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgtt                                      <  2:572624/99‑1 (MQ=255)                                                            tCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgtt                                   >  1:1130457/1‑100 (MQ=255)                                                             cGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTa                                  >  1:98361/1‑100 (MQ=255)                                                             cGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTa                                  >  2:2069606/1‑100 (MQ=255)                                                             cGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTa                                  >  1:846495/1‑100 (MQ=255)                                                              gATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAc                                 <  1:998181/100‑1 (MQ=255)                                                               aTAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAcc                                <  1:713672/100‑1 (MQ=255)                                                               aTAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCGCGACGGCTATGAGAAATGGGGGTTAc                                 >  2:1939881/1‑99 (MQ=255)                                                                tAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCCGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCt                               >  1:641802/1‑100 (MQ=255)                                                                 aaaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTc                              <  2:719514/100‑1 (MQ=255)                                                                  aaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCt                               >  2:1556133/1‑98 (MQ=255)                                                                    aCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGc                           >  1:855479/1‑100 (MQ=255)                                                                    aCCGGAGCCGTGCAGCCGGGGAGTCGGGACACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGCGAAATGGTTGTTACCTc                              >  2:1223772/1‑97 (MQ=255)                                                                     ccGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCt                          <  1:530965/100‑1 (MQ=255)                                                                     ccGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCt                          >  2:530958/1‑100 (MQ=255)                                                                      cGGGGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCtt                         <  2:1589091/100‑1 (MQ=255)                                                                          gCCGTGCAGCCGGGGAGACGGGGGACCTGGTGTTATAAATTCAGCGGCAACGAGAACCAGATAGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAc                     <  2:95789/100‑1 (MQ=255)                                                                           ccGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAcc                    <  2:641415/100‑1 (MQ=255)                                                                           ccGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAcc                    <  1:503616/100‑1 (MQ=255)                                                                           ccGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAcc                    <  1:503615/100‑1 (MQ=255)                                                                             gTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCg                  <  2:2468304/100‑1 (MQ=255)                                                                               gCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaa                <  2:2237978/100‑1 (MQ=255)                                                                               gCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaa                >  1:985141/1‑100 (MQ=255)                                                                                cAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaac               <  1:1323221/100‑1 (MQ=255)                                                                                  gCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaacaa             >  1:991461/1‑100 (MQ=255)                                                                                   ccGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaacaac            >  1:1926455/1‑100 (MQ=255)                                                                                    cGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACg           >  1:1474944/1‑100 (MQ=255)                                                                                     ggggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGc          >  2:1853256/1‑100 (MQ=255)                                                                                       ggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTg        <  1:2474977/100‑1 (MQ=255)                                                                                        gAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTgg       <  2:1174697/100‑1 (MQ=255)                                                                                          gTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTa     <  1:2440444/100‑1 (MQ=255)                                                                                            cGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACg   >  1:2425465/1‑100 (MQ=255)                                                                                             ggggCACCTGGCGTTATAAATTCGGCGGTAACGAGGACCAGCTCGACGGCAATGAGAAATGGATGGTAAGTCAGCTTAAACCCGAACAACGCTGGTACg   >  2:245798/1‑99 (MQ=255)                                                                                             ggggCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGCCGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACACCGCTGGtcc    >  2:1212906/1‑96 (MQ=255)                                                                                             ggggCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACgg  >  1:2123019/1‑100 (MQ=255)                                                                                                |                                                                                                 GGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGG  >  NC_000913/520898‑521089

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