Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,215,832 (C)5→6 coding (645/732 nt) yehW ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,215,8271.C97.9% 336.5 / 0.3 95coding (650/732 nt)yehWputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (50/43);  total (52/43)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 102 total aligned reads displayed.

GCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTC  >  NC_000913/2215728‑2215919
                                                                                                    |                                                                                             
gCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACg                                                                                                <  2:136362/100‑1 (MQ=255)
    gCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑ccc‑c                                                                                          <  1:2229329/100‑1 (MQ=255)
       ggTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑ccc‑ccct                                                                                       >  2:2233233/1‑98 (MQ=255)
        gTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑ccc‑ccctg                                                                                      >  1:385979/1‑97 (MQ=255)
          tataCCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑ccc‑ccctgga                                                                                    <  1:996075/100‑6 (MQ=255)
          tataCCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑ccc‑ccctgga                                                                                    <  1:285792/100‑6 (MQ=255)
          tataCCTTTATTTTGCGTGCTGGCTAAGCACCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑ccc‑ccctgga                                                                                    >  2:1777970/1‑95 (MQ=255)
           ataCCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGAt                                                                                    <  2:2097719/100‑1 (MQ=255)
            taCCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATc                                                                                   <  1:1243176/100‑1 (MQ=255)
            taCCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATc                                                                                   <  2:2532911/100‑1 (MQ=255)
            taCCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATc                                                                                   <  1:1869096/100‑1 (MQ=255)
              ccTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATcac                                                                                 <  1:414783/100‑1 (MQ=255)
               cTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATcaca                                                                                >  1:2640706/1‑100 (MQ=255)
               cTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATcaca                                                                                >  1:1893459/1‑100 (MQ=255)
                tttATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACAt                                                                               <  2:2846879/100‑1 (MQ=255)
                    ttttGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATAcgc                                                                           <  2:2042585/100‑1 (MQ=255)
                     tttGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATAcgcg                                                                          >  1:692413/1‑100 (MQ=255)
                     tttGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATAcgcg                                                                          >  2:1594516/1‑100 (MQ=255)
                     tttGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATAcgcg                                                                          >  1:519778/1‑100 (MQ=255)
                      ttGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCgg                                                                         <  1:849175/100‑1 (MQ=255)
                       tGCGTGCTGGCTAAGCGCCTGTACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGt                                                                        >  1:1360566/1‑100 (MQ=255)
                       tGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGt                                                                        >  2:2538696/1‑100 (MQ=255)
                        gCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTa                                                                       <  1:2288742/100‑1 (MQ=255)
                         cGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTAt                                                                      >  1:268051/1‑100 (MQ=255)
                         cGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTAt                                                                      >  1:267983/1‑100 (MQ=255)
                            gCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTaa                                                                   <  1:144530/100‑1 (MQ=255)
                             cTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTaaa                                                                  <  2:2460406/100‑1 (MQ=255)
                              tGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAAt                                                                 <  1:2245639/100‑1 (MQ=255)
                              tGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAAt                                                                 >  1:815766/1‑100 (MQ=255)
                               ggCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACTCCCC‑CCTGGATCACATATGCGGTATTAAATc                                                                <  1:2074341/100‑1 (MQ=255)
                               ggCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATc                                                                >  2:2654259/1‑100 (MQ=255)
                                gCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATcc                                                               <  1:601526/100‑1 (MQ=255)
                                gCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATcc                                                               <  1:2246909/100‑1 (MQ=255)
                                 cTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCg                                                              <  1:416138/100‑1 (MQ=255)
                                  tAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGc                                                             >  2:1299441/1‑100 (MQ=255)
                                   aaGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCt                                                            >  2:920358/1‑100 (MQ=255)
                                    aGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTa                                                           >  1:1206377/1‑100 (MQ=255)
                                    aGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTa                                                           <  1:1594533/100‑1 (MQ=255)
                                     gcgcCTGCACCAGCCTTTCAAACAGGCGTTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTa                                                           <  2:2963248/99‑1 (MQ=255)
                                       gcCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGc                                                        <  2:1625974/100‑1 (MQ=255)
                                        ccTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGccc                                                      >  2:505996/1‑100 (MQ=255)
                                          tGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCg                                                     >  2:459725/1‑100 (MQ=255)
                                          tGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCg                                                     >  2:375897/1‑100 (MQ=255)
                                          tGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTCGATCACATACGCGGTATTAAATCCGCTAAGccc                                                      <  2:1850060/99‑1 (MQ=255)
                                            cACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGcccatc                                                   <  1:2299585/100‑4 (MQ=255)
                                            cACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCgat                                                   >  1:1462307/1‑100 (MQ=255)
                                             aCCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCgatg                                                  <  2:2264087/100‑1 (MQ=255)
                                             aCCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCgatg                                                  <  2:1883580/100‑1 (MQ=255)
                                             aCCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCgatg                                                  >  2:1775613/1‑100 (MQ=255)
                                              ccAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCgatga                                                 >  1:2672438/1‑100 (MQ=255)
                                                 gCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCgatgatga                                              <  1:1558973/100‑1 (MQ=255)
                                                  ccTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCgatgatgat                                             <  2:1439601/100‑1 (MQ=255)
                                                    tttCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTACGCCCGATGATGATgg                                           >  1:965245/1‑100 (MQ=255)
                                                    tttCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATgg                                           >  1:2870863/1‑100 (MQ=255)
                                                     ttCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATggg                                          >  1:988136/1‑100 (MQ=255)
                                                     ntCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATggg                                          >  1:1693178/2‑100 (MQ=255)
                                                      tCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGc                                         <  2:819336/100‑1 (MQ=255)
                                                       cAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCg                                        >  2:1086450/1‑100 (MQ=255)
                                                       cAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCg                                        >  2:2430189/1‑100 (MQ=255)
                                                       cAAAAAGGCGGTCAGCGATAATCACCGCCAGGGCGAACAGTAACG‑CCCTCCTGGATCAAATACGAGGTATTAAAGCCGCGAAGCCCGatgatgatg                                            >  2:1339606/1‑96 (MQ=255)
                                                           cAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTAcc                                    >  2:2768275/1‑100 (MQ=255)
                                                           cAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTAcc                                    >  2:2897731/1‑100 (MQ=255)
                                                            aGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTAccc                                   <  1:2659530/100‑1 (MQ=255)
                                                                ggTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCg                               >  2:910326/1‑100 (MQ=255)
                                                                     gCGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTg                          >  1:429636/1‑100 (MQ=255)
                                                                      cGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTgg                         <  2:1960890/100‑1 (MQ=255)
                                                                      cGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTgg                         <  2:2864260/100‑1 (MQ=255)
                                                                      cGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTgg                         >  1:591987/1‑100 (MQ=255)
                                                                      cGATGATCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTgg                         >  1:1637002/1‑100 (MQ=255)
                                                                       gatgatCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGc                        >  2:2781181/1‑100 (MQ=255)
                                                                        atgatCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGcc                       <  1:2538719/100‑1 (MQ=255)
                                                                        atgatCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGcc                       >  1:1516872/1‑100 (MQ=255)
                                                                           atCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCt                    >  1:655816/1‑100 (MQ=255)
                                                                           atCGCCGCCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCt                    >  1:959727/1‑100 (MQ=255)
                                                                             cgccgcCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTAc                  >  2:2133042/1‑100 (MQ=255)
                                                                             cgccgcCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTAc                  >  2:499554/1‑100 (MQ=255)
                                                                             cgccgcCAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTAc                  <  2:1543688/100‑1 (MQ=255)
                                                                              gccgccAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTAcc                 <  2:2640682/100‑1 (MQ=255)
                                                                              gccgccAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTAcc                 >  1:376111/1‑100 (MQ=255)
                                                                              gccgccAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTAcc                 >  1:1917696/1‑100 (MQ=255)
                                                                              gccgccAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTAcc                 <  1:2806686/100‑1 (MQ=255)
                                                                               ccgccAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCg                <  2:385977/100‑1 (MQ=255)
                                                                                cgccAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGt               >  1:2279495/1‑100 (MQ=255)
                                                                                  ccAGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTg             <  1:1480370/100‑1 (MQ=255)
                                                                                    aGTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAg           <  1:2781208/100‑1 (MQ=255)
                                                                                     gTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAgn          <  2:2427986/100‑2 (MQ=255)
                                                                                     gTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAgg          <  2:2498000/100‑1 (MQ=255)
                                                                                     gTGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAgg          >  1:725899/1‑100 (MQ=255)
                                                                                      tGCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGc         <  2:388882/100‑1 (MQ=255)
                                                                                       gCCACCAGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCg        >  1:1326982/1‑100 (MQ=255)
                                                                                        ccaccaGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGa       <  2:1206367/100‑1 (MQ=255)
                                                                                         caccaGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAt      >  2:1053918/1‑100 (MQ=255)
                                                                                         caccaGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAt      <  1:2041748/100‑1 (MQ=255)
                                                                                          accaGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtc     <  1:1299452/100‑1 (MQ=255)
                                                                                          accaGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtc     >  1:2019010/1‑100 (MQ=255)
                                                                                            caGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGTCCCTACCGTTGAGGCGAtcgt   >  2:2283643/1‑100 (MQ=255)
                                                                                            caGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtcgt   >  1:1508047/1‑100 (MQ=255)
                                                                                            caGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtcgt   >  2:1888914/1‑100 (MQ=255)
                                                                                             aGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtcgtc  >  2:1703751/1‑100 (MQ=255)
                                                                                             aGTAACGCCCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtcgtc  <  1:1132008/100‑1 (MQ=255)
                                                                                                    |                                                                                             
GCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCC‑CCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTC  >  NC_000913/2215728‑2215919

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