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_000913520,9930TG98.4% 394.5 / ‑3.2 118V193G (GTC→GGC) ybbPputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (71/45);  total (73/45)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 118 total aligned reads displayed.

CCGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACG  >  NC_000913/520896‑521088
                                                                                                 |                                                                                               
  ggATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTTTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAcc                                                                                             <  1:1769962/100‑1 (MQ=255)
  ggATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAcc                                                                                             >  1:1047985/1‑100 (MQ=255)
   gATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCt                                                                                            >  2:2281491/1‑100 (MQ=255)
   gATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGAGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCt                                                                                            >  2:827636/1‑100 (MQ=255)
      tCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGTCTGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGc                                                                                         <  1:77845/100‑1 (MQ=255)
       ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                        >  2:1842544/1‑100 (MQ=255)
       ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGACGGGGCACCTGGCg                                                                                        >  1:2376688/1‑100 (MQ=255)
        cGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGt                                                                                       >  1:2097299/1‑100 (MQ=255)
        cGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGt                                                                                       <  2:274388/100‑1 (MQ=255)
          gTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTta                                                                                     <  2:447158/100‑1 (MQ=255)
              tAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATaaa                                                                                 >  1:316827/1‑100 (MQ=255)
               aaCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAAt                                                                                >  1:2426549/1‑100 (MQ=255)
               aaCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAAt                                                                                <  2:241163/100‑1 (MQ=255)
                aCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGCCGGGGCACCTGGCGTTATAAAtt                                                                               >  2:1598166/1‑100 (MQ=255)
                  cccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACGTGGCGTTATAAATGcg                                                                             >  2:3025687/1‑100 (MQ=255)
                  cccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcg                                                                             <  2:1389745/100‑1 (MQ=255)
                     ttCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcg                                                                          <  2:2857314/100‑1 (MQ=255)
                      tCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcgg                                                                         <  1:2460310/100‑1 (MQ=255)
                        cAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGTGTTATAAATTCGGCGGCa                                                                       >  2:1070084/1‑100 (MQ=255)
                        cAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCa                                                                       <  1:274793/100‑1 (MQ=255)
                        cAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAACTGGAGTTATAAATTCGGCGGCa                                                                       >  1:1982832/1‑100 (MQ=255)
                          naTGGCTCCGCGTCTGATGATGAATCTGGCGGACGTCGATAAAACCGGCACCGTGCAGCCGGGGAGACGGGGCACCTGGCGTTATAAATTCGGCGGCaa                                                                      >  1:1183404/2‑99 (MQ=255)
                           aTGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACg                                                                    >  2:2339735/1‑100 (MQ=255)
                             ggCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCATCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACgag                                                                  <  2:1235658/100‑1 (MQ=255)
                              gCTCCGCGTCTGATGATGAATCTGGCTGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACgaga                                                                 <  1:115993/100‑1 (MQ=255)
                                 ccGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAAcc                                                              >  1:984171/1‑100 (MQ=255)
                                  cgcgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCa                                                             <  1:2811713/100‑1 (MQ=255)
                                    cgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGc                                                           >  1:2655337/1‑100 (MQ=255)
                                      tCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTc                                                         >  1:623753/1‑100 (MQ=255)
                                      tCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGCGAGTCTTGGCTCCTGGTGTTATAAATTCGGCGGCAACGAGAAcc                                                              <  2:1863211/95‑1 (MQ=255)
                                       cTGATGATGAATCTGTCGGATGTCGATAAAACCGGAGCCGAGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGc                                                           >  2:1138894/1‑97 (MQ=255)
                                       cTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCg                                                        >  1:644362/1‑100 (MQ=255)
                                        tgatgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGa                                                       >  1:119257/1‑100 (MQ=255)
                                           tgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACg                                                     >  1:402360/1‑99 (MQ=255)
                                            gatgaATCTGGCGGATGTCGATAAAACCGGAGCTGTGCAGCCGGGGTGTTGGGGTAGTTGGTGTTATAAATTTGGTGGCAACGAGAACCAGCTCGACgg                                                    >  1:2530183/1‑99 (MQ=255)
                                            gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGTCGGGGAGTGGGGGGATTTGGGGTTATAAATTAGGCGGCAACGAGAACCAGCTCGACGGc                                                   >  1:2529610/1‑100 (MQ=255)
                                            gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGc                                                   <  1:1456023/100‑1 (MQ=255)
                                            gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGc                                                   >  1:2407105/1‑100 (MQ=255)
                                              tgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTa                                                 <  2:1711337/100‑1 (MQ=255)
                                              tgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTa                                                 >  2:2617352/1‑100 (MQ=255)
                                              tgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTa                                                 >  1:1226340/1‑100 (MQ=255)
                                                aaTCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATg                                               >  1:469339/1‑100 (MQ=255)
                                                aaTCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGACTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATg                                               <  2:1946410/100‑1 (MQ=255)
                                                   cTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATgaga                                            <  1:500664/100‑1 (MQ=255)
                                                     ggCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGTAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATga                                              <  2:412853/96‑1 (MQ=255)
                                                      gCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAAt                                         <  1:1202024/100‑1 (MQ=255)
                                                      gCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAAt                                         <  1:941814/100‑1 (MQ=255)
                                                       cGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATg                                        <  1:2934297/100‑1 (MQ=255)
                                                        ggATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATgg                                       >  2:1884949/1‑100 (MQ=255)
                                                         gATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgt                                      <  2:2058084/100‑1 (MQ=255)
                                                           tGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttg                                    <  1:2157317/100‑1 (MQ=255)
                                                           tGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttg                                    >  1:2202758/1‑100 (MQ=255)
                                                           tGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttg                                    <  2:2372499/100‑1 (MQ=255)
                                                           tGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttg                                    >  1:55291/1‑100 (MQ=255)
                                                            ggcaaTAAAACCGGAGCCGTGTAACCGGGGAATCGGGGAACCTGGCGTTATAAATTCGGCGGCAACGAAAACCAGCTCGATGGCTATGAGAAATGgtt                                     <  2:1183174/94‑1 (MQ=255)
                                                            gTCGATAAAACCGGAGCCGGGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgt                                   <  1:884552/100‑1 (MQ=255)
                                                             tCGATAAAACCGGATCCGTGCAGCCGGGGAGTCGGGGTACCTGGCGTTATAATTTCGGCAGCAACGAGAACCAGCTCGACGGCTTTGAGAAATGgttg                                    <  2:1875422/98‑1 (MQ=255)
                                                             tCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgtt                                  >  2:2858770/1‑100 (MQ=255)
                                                              cGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGACATGGTTGTTa                                 >  2:751814/1‑100 (MQ=255)
                                                               gATAAAACCGAAGCCGTGCAGCCGGGGAGTCTGGGCACCGGACGTTATAAATTCGGCGGCGAAGAGAACCAGCGCGACGTGTATGAGAAAttg                                       >  2:1766622/1‑91 (MQ=255)
                                                                tctaaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTa                                 <  2:1222387/95‑1 (MQ=255)
                                                                aTAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAcc                               >  1:2950597/1‑100 (MQ=255)
                                                                 tAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCt                              >  1:2477321/1‑100 (MQ=255)
                                                                 tAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCt                              >  1:8965/1‑100 (MQ=255)
                                                                  aaaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTc                             >  1:2079591/1‑100 (MQ=255)
                                                                   aaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTaaga                              >  1:1222398/1‑95 (MQ=255)
                                                                   aaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTc                             >  2:1577375/1‑99 (MQ=255)
                                                                    aaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCACAg                           >  2:1846595/1‑100 (MQ=255)
                                                                    aaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGCTGTTACCTCa                            >  2:1844199/1‑99 (MQ=255)
                                                                     aCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAg                           >  2:1700681/1‑99 (MQ=255)
                                                                     aCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGc                          >  1:472091/1‑100 (MQ=255)
                                                                      ccGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCt                         >  1:1090560/1‑100 (MQ=255)
                                                                      ccGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCt                         <  1:2238943/100‑1 (MQ=255)
                                                                       cGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCtt                        >  2:2088330/1‑100 (MQ=255)
                                                                       cGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGGTGTTACCTCAGCtt                        >  2:2064716/1‑100 (MQ=255)
                                                                          aGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGTCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTaaa                     <  2:1968059/100‑1 (MQ=255)
                                                                           gCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAc                    >  1:1015069/1‑100 (MQ=255)
                                                                           gCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAc                    >  1:1989236/1‑100 (MQ=255)
                                                                           gCCGGGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCTGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAc                    <  2:1296897/100‑1 (MQ=255)
                                                                            ctgTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAcc                   <  1:730934/98‑1 (MQ=255)
                                                                            ccGTGCATCCGGGGAGTCGTGGCCCCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAcc                   <  2:522432/100‑1 (MQ=255)
                                                                            ccGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAcc                   <  2:87995/100‑1 (MQ=255)
                                                                            ccGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAcc                   <  2:2672440/100‑1 (MQ=255)
                                                                            ccGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCCGCAACGCGAACCAGCTCGACGGCTATGAGAAATAGTTGTTACCTCAGCTTAAAcc                   >  1:1670888/1‑100 (MQ=255)
                                                                             ggTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAccc                  >  1:716457/2‑100 (MQ=255)
                                                                             cGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTCCCTCAGCTTAAAccc                  >  1:3006014/1‑100 (MQ=255)
                                                                             cGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAccc                  >  2:226448/1‑100 (MQ=255)
                                                                              gTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAccag                 >  2:730930/1‑98 (MQ=255)
                                                                               tGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGa                <  1:117122/100‑1 (MQ=255)
                                                                                gCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaa               <  1:1319547/100‑1 (MQ=255)
                                                                                   gCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTTGGCTTAAACCCGa                >  2:92420/1‑96 (MQ=255)
                                                                                   gCCGGGGAGTCGGAGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCAATGAGAAATGGGTGTTACCTCAGCTTAAACCCGaaca             >  2:2426430/1‑99 (MQ=255)
                                                                                    ccGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaacaac           >  1:1888662/1‑100 (MQ=255)
                                                                                      ggggggTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTGTGAGAAATGGTTGTTGACTCGGCGTAAGCCCGAACAAGGc         >  2:2890261/1‑100 (MQ=255)
                                                                                       gggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCt        <  1:165225/100‑1 (MQ=255)
                                                                                        ggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTg       <  2:316824/100‑1 (MQ=255)
                                                                                          aGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCACAACTTAAACCCGAAAAACGCTgg      >  1:2353741/1‑99 (MQ=255)
                                                                                           gTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTa    >  1:1425024/1‑100 (MQ=255)
                                                                                           gTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGATGGTa    >  1:1336415/1‑100 (MQ=255)
                                                                                             cGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGCACg  >  2:1111192/1‑100 (MQ=255)
                                                                                                 |                                                                                               
CCGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACG  >  NC_000913/520896‑521088

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