Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 520,993 T→G 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_000913520,9930TG99.2% 441.5 / NA 129V193G (GTC→GGC) ybbPputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base G (69/58);  minor base A (1/0);  total (71/58)
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

Only 100 of 133 total aligned reads displayed.

CCGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGG  >  NC_000913/520896‑521089
                                                                                                 |                                                                                                
  ggATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAcc                                                                                              >  2:1379747/1‑100 (MQ=255)
     ttCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTgg                                                                                           >  2:373060/1‑100 (MQ=255)
       ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         >  2:2870307/1‑100 (MQ=255)
       ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         <  1:2966469/100‑1 (MQ=255)
        cGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGt                                                                                        >  2:1209704/1‑100 (MQ=255)
         ggTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGtt                                                                                       >  2:1264507/1‑100 (MQ=255)
         ggTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGtt                                                                                       <  1:1378666/100‑1 (MQ=255)
         ggTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGGCGGGGCACCTGGCGtt                                                                                       >  1:590716/1‑100 (MQ=255)
               aaCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAAt                                                                                 <  2:1245828/100‑1 (MQ=255)
               aaCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAAt                                                                                 >  1:1656196/1‑100 (MQ=255)
               aaCCCCTTCCAGATGGATCCGCGTCTGATGATGAATCCCGCGGATGTCGATAAAACCGGCGCCGTGCAGCCGGGGAGCAGGGGCACCTCGCGGTATAAAt                                                                                 >  1:1884368/1‑100 (MQ=255)
                aCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAAtt                                                                                >  1:383687/1‑100 (MQ=255)
                aCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAAtt                                                                                >  2:1690552/1‑100 (MQ=255)
                aCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCAGGCGTTATAAAt                                                                                 >  1:2503419/1‑99 (MQ=255)
                 ccccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGCGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTc                                                                               >  1:2442080/1‑100 (MQ=255)
                 ccccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTc                                                                               <  1:1500599/100‑1 (MQ=255)
                    cTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGACAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggc                                                                            >  1:2004890/1‑100 (MQ=255)
                     ttCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcg                                                                           <  2:2208374/100‑1 (MQ=255)
                     ttCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcg                                                                           <  2:424775/100‑1 (MQ=255)
                     ttCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggc                                                                            <  1:1537674/99‑1 (MQ=255)
                      tCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcgn                                                                          <  2:1306298/100‑2 (MQ=255)
                         aGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTCATAAATTCggcggcca                                                                       >  1:841766/1‑98 (MQ=255)
                          gATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAAc                                                                      <  1:3215296/100‑1 (MQ=255)
                          gATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAAc                                                                      >  1:223799/1‑100 (MQ=255)
                           aTGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACg                                                                     <  2:223802/100‑1 (MQ=255)
                            tGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACga                                                                    <  1:1975042/100‑1 (MQ=255)
                             ggCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACgag                                                                   <  2:1966841/100‑1 (MQ=255)
                               cTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGaa                                                                 >  2:1012670/1‑100 (MQ=255)
                               cTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGaa                                                                 >  2:85288/1‑100 (MQ=255)
                               cTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGaa                                                                 >  1:1237655/1‑100 (MQ=255)
                                  cgcgTCTGATGATGAATCTGGCGGATGTCGATAAACCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCa                                                              <  2:3474842/100‑1 (MQ=255)
                                  cgcgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCa                                                              >  1:1062193/1‑100 (MQ=255)
                                   gcgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAg                                                             >  1:1555948/1‑100 (MQ=255)
                                    cgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGc                                                            <  1:1553193/100‑1 (MQ=255)
                                    cgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGc                                                            <  1:1690529/100‑1 (MQ=255)
                                       cTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTc                                                          <  1:1340236/99‑1 (MQ=255)
                                       cTGATGAAGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGTAGGCGTGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTc                                                          <  1:2054715/99‑1 (MQ=255)
                                         gatgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGAc                                                       >  1:3220081/1‑100 (MQ=255)
                                          atgatgaATCTGGCGGATGTCGATAAAGCCGGAGCCGTGCAGACGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACg                                                      <  2:2951285/100‑1 (MQ=255)
                                           tgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACgg                                                     >  2:1341222/1‑100 (MQ=255)
                                           tgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACgg                                                     <  1:3322473/100‑1 (MQ=255)
                                           tgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACgg                                                     >  2:1005436/1‑100 (MQ=255)
                                           tgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACgg                                                     >  2:710700/1‑100 (MQ=255)
                                            gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGc                                                    >  2:2829641/1‑100 (MQ=255)
                                             aggaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCt                                                   >  2:3225308/3‑100 (MQ=255)
                                               gaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGGGAACCAGCTCGACGGCTAt                                                 >  1:354010/1‑100 (MQ=255)
                                                aaTCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATg                                                >  2:1438118/1‑100 (MQ=255)
                                                 aTCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATga                                               >  2:1736258/1‑100 (MQ=255)
                                                   ntGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATgaga                                             >  2:1216874/2‑100 (MQ=255)
                                                   cTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATgaga                                             >  1:1761544/1‑100 (MQ=255)
                                                     ngCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATga                                               >  1:3448031/2‑96 (MQ=255)
                                                      gCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAAt                                          >  2:2546957/1‑100 (MQ=255)
                                                       cGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAAt                                          <  2:3390584/99‑1 (MQ=255)
                                                       cGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATg                                         <  2:2052920/100‑1 (MQ=255)
                                                       cGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATg                                         >  2:2019902/1‑100 (MQ=255)
                                                         gATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgt                                       >  2:897919/1‑100 (MQ=255)
                                                             tCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgtt                                   <  1:3095491/100‑1 (MQ=255)
                                                             tCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgtt                                   <  2:2624465/100‑1 (MQ=255)
                                                              cGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTa                                  >  2:518188/1‑100 (MQ=255)
                                                              cGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTa                                  >  1:2624920/1‑100 (MQ=255)
                                                              cGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTAATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTa                                  >  1:1933799/1‑100 (MQ=255)
                                                               gATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAc                                 <  2:3441634/100‑1 (MQ=255)
                                                               gATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAc                                 >  1:349894/1‑100 (MQ=255)
                                                                aTAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAcc                                <  1:2199853/100‑1 (MQ=255)
                                                                 tAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCt                               >  1:2622724/1‑100 (MQ=255)
                                                                  aaaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTc                              <  1:1822669/100‑1 (MQ=255)
                                                                   aaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCa                             <  1:2829597/100‑1 (MQ=255)
                                                                   aaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCa                             >  1:1544356/1‑100 (MQ=255)
                                                                      ccGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCt                          <  1:32101/100‑1 (MQ=255)
                                                                      ccGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCt                          <  1:1496179/100‑1 (MQ=255)
                                                                      ccGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCt                          >  1:1880157/1‑100 (MQ=255)
                                                                            ccGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAGATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAcc                    <  2:2442113/100‑1 (MQ=255)
                                                                             cGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAc                     <  2:1403558/98‑1 (MQ=255)
                                                                              gTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCg                  <  1:603426/100‑1 (MQ=255)
                                                                              gTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAAGGGTTGTTACCTCAGCTTAAACCCg                  >  1:1582130/1‑100 (MQ=255)
                                                                               tGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGa                 <  2:2158813/100‑1 (MQ=255)
                                                                               tGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGa                 <  2:505995/100‑1 (MQ=255)
                                                                                gCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaa                <  2:2706095/100‑1 (MQ=255)
                                                                                 cAGTCGGGGAGTCGGGGCACCTGGCGTTATATATTCGGCGGCAACGAGAGCAAGCTCGACATCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaac               <  1:1601903/100‑1 (MQ=255)
                                                                                  aGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGGCGGCTATGAGGAATGGTTGTTACCTCAGCTTAAACCAGaaca              >  1:152500/1‑100 (MQ=255)
                                                                                  aGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaaca              >  2:90483/1‑100 (MQ=255)
                                                                                  aGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTGGTTACCGCAGCTTAAACCCGaaca              >  1:1079926/1‑100 (MQ=255)
                                                                                   gCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaacaa             >  1:3014548/1‑100 (MQ=255)
                                                                                    ccGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaacaac            <  2:211787/100‑1 (MQ=255)
                                                                                     cGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACg           >  2:2008582/1‑100 (MQ=255)
                                                                                     cGGGGAGTCGGGGCACCTGGCGTTATAAATTCGACGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAGCg           >  1:255380/1‑100 (MQ=255)
                                                                                        ggAGTCGGGGCGCCTGGCGTTATAAATTCGGCGGCAACGAGAACCTGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTg        <  1:848335/100‑1 (MQ=255)
                                                                                        ggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTg        <  1:323138/100‑1 (MQ=255)
                                                                                        ggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTGATTACCTCAGCTTAAACCCGAACAACGCTg        >  1:1483619/1‑100 (MQ=255)
                                                                                         gAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTgg       <  2:659921/100‑1 (MQ=255)
                                                                                         gAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTg        >  1:3127065/1‑99 (MQ=255)
                                                                                          aGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGt      >  1:3254380/1‑100 (MQ=255)
                                                                                           gTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTa     >  2:602344/1‑100 (MQ=255)
                                                                                           gTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGGTGTTACCTCAGCTTAAACCCGAACAACGCTGGTa     >  1:6976/1‑100 (MQ=255)
                                                                                            tCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTAc    <  2:2064918/100‑1 (MQ=255)
                                                                                             cGGGGCACCTGGGGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTa     <  1:1078283/98‑1 (MQ=255)
                                                                                              ggggCACCTGGCGTTATAAATTCGGCGGCACCGAGAACCCGCTCGACGGCGATGAGAAATGGACGTTACCACAGCTTACGCCCGAGAAACGCTGGTACgg  >  1:3012876/1‑100 (MQ=255)
                                                                                              ggggCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACgg  <  1:935466/100‑1 (MQ=255)
                                                                                              ggggCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACgg  >  2:1698593/1‑100 (MQ=255)
                                                                                              ggggCACCTGGCGATATAAATTCGGCGGCACCGAGAACCTGCTCGACGGCTATGAGAAATGGACAAGACCACAGCTTAAAGCCGAGCAACGCgggaa     >  1:3176181/1‑95 (MQ=255)
                                                                                                 |                                                                                                
CCGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGG  >  NC_000913/520896‑521089

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