Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 825,545 T→G V188V (GTA→GTC ybhR ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913825,5450TG98.6% 496.5 / ‑2.0 143V188V (GTA→GTCybhRputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (68/72);  total (68/75)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 143 total aligned reads displayed.

CGATGAAGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATTACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTACCACCAGCTCGC  >  NC_000913/825448‑825639
                                                                                                 |                                                                                              
cGATGAAGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGCCGGAAAGTGAAGTGACGATCATGAc                                                                                              >  1:1310697/1‑100 (MQ=255)
cGATGAAGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGAc                                                                                              <  1:143205/100‑1 (MQ=255)
 gATGAAGATCTGCCAGGTGGTTAGGGGCTAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTAGCGGGCGACGGAAAGTGAAGTGACGATCATGACg                                                                                             >  2:188535/1‑100 (MQ=255)
 gATGAAGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACg                                                                                             >  2:430592/1‑100 (MQ=255)
 gATGAAGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACg                                                                                             >  2:3089412/1‑100 (MQ=255)
      aGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGAt                                                                                        >  1:343189/1‑100 (MQ=255)
       gATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATa                                                                                       >  1:1245494/1‑100 (MQ=255)
       gATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATa                                                                                       >  2:3355395/1‑100 (MQ=255)
         tCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgt                                                                                     <  2:3220158/100‑1 (MQ=255)
         tCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgt                                                                                     >  2:2356111/1‑100 (MQ=255)
           tGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgtgg                                                                                   <  1:438892/100‑1 (MQ=255)
           tGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgtgg                                                                                   <  2:395583/100‑1 (MQ=255)
           tGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgtgg                                                                                   <  2:2260822/100‑1 (MQ=255)
            gCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgtggt                                                                                  >  1:1287062/1‑100 (MQ=255)
            gCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgtggt                                                                                  <  1:1900128/100‑1 (MQ=255)
            gCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgtggg                                                                                  >  1:4082804/1‑99 (MQ=255)
              cAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGa                                                                                <  1:890230/100‑1 (MQ=255)
                 gtggtgAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGTCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCa                                                                             <  1:1846091/100‑1 (MQ=255)
                  tggtgAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCAt                                                                            <  1:2735701/100‑1 (MQ=255)
                  tggtgAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCAt                                                                            <  2:2282002/100‑1 (MQ=255)
                  tggtgAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCAACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCAt                                                                            <  2:355129/100‑1 (MQ=255)
                    gtgAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATgg                                                                          <  1:3527915/100‑1 (MQ=255)
                     tgAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGc                                                                         >  1:2216374/1‑100 (MQ=255)
                     tgAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGc                                                                         >  1:931484/1‑100 (MQ=255)
                      gAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCg                                                                        >  2:1459157/1‑100 (MQ=255)
                      gAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCg                                                                        >  1:3474012/1‑100 (MQ=255)
                        gcggcgAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGAt                                                                      <  2:2902263/100‑1 (MQ=255)
                         cggcgAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATc                                                                     >  1:734577/1‑100 (MQ=255)
                          ggcgAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCa                                                                    >  2:2176930/1‑100 (MQ=255)
                           gcgAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAg                                                                   <  1:2392855/100‑1 (MQ=255)
                            cgAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGt                                                                  <  1:3612757/100‑1 (MQ=255)
                             taaaCCAGTAGCGGATCGAGCGGACGTTGTTGACGCTCGCGGGCGACGGTAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTg                                                                 <  2:3938670/99‑1 (MQ=255)
                             gAAACCAGTAGCTGATCGAGCGTAGCTTGTCGACGTTCGCGGGCGCCGGAAAGTGAAGCGACGATCATGACGCCGACAGTGGTGATCATGGCGATCAGTg                                                                 >  2:2711517/1‑100 (MQ=255)
                                aCCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACg                                                              >  1:1080856/1‑100 (MQ=255)
                                  cAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGc                                                            >  1:1979061/1‑100 (MQ=255)
                                    gTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGcac                                                          <  2:2754928/100‑1 (MQ=255)
                                    gTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGcac                                                          <  1:1477667/100‑1 (MQ=255)
                                     tAGCTGATCGAGCGTACCTTGTTCACGTTCGCGTGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGcacc                                                         <  2:2981730/100‑1 (MQ=255)
                                     tAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGcacc                                                         <  1:4042729/100‑1 (MQ=255)
                                      aGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGcacca                                                        >  1:3635578/1‑100 (MQ=255)
                                          gATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACaaa                                                    <  2:3601420/100‑1 (MQ=255)
                                           aTCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAAc                                                   <  1:2600141/100‑1 (MQ=255)
                                           aTCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAAc                                                   >  2:665753/1‑100 (MQ=255)
                                           aTCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAAc                                                   <  1:922393/100‑1 (MQ=255)
                                            tCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAAcc                                                  <  1:3145951/100‑1 (MQ=255)
                                             cGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCa                                                 >  2:3102660/1‑100 (MQ=255)
                                              gAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCAt                                                <  1:2174475/100‑1 (MQ=255)
                                               aGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCAtt                                               >  2:3054916/1‑100 (MQ=255)
                                                gCGTACCTTGTTCACGTTCGCGGGCGACGCAAAGTGAAGTGACTATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCAttt                                              <  1:3723347/100‑1 (MQ=255)
                                                 cGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACTACAAACCATTTg                                             <  1:3518840/100‑1 (MQ=255)
                                                 cGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTg                                             <  2:2187656/100‑1 (MQ=255)
                                                 cGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTg                                             <  2:2437102/100‑1 (MQ=255)
                                                 cGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTg                                             >  2:3239251/1‑100 (MQ=255)
                                                 cGTACCTTGTTCACGTTCGCGGGCGACAGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTg                                             <  1:2834679/100‑1 (MQ=255)
                                                  gTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGt                                            <  2:1394824/100‑1 (MQ=255)
                                                   tACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTa                                           >  1:2970940/1‑100 (MQ=255)
                                                    aCCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAg                                          <  2:1547640/100‑1 (MQ=255)
                                                     ccTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGt                                         >  1:171659/1‑100 (MQ=255)
                                                     ccTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGt                                         >  1:3110242/1‑100 (MQ=255)
                                                     ccTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGt                                         >  1:1865531/1‑100 (MQ=255)
                                                      cTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTc                                        <  2:583787/100‑1 (MQ=255)
                                                       ttGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGATGTTCATGGCGATCAGTGACGGCACCACAATCCATTTGTAGTCg                                       <  1:3869160/100‑1 (MQ=255)
                                                         gTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCCCAAACCATTTGTAGTCgag                                     >  2:3247435/1‑100 (MQ=255)
                                                          ttCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCgaga                                    >  2:736176/1‑100 (MQ=255)
                                                           tCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGAt                                   <  2:1222933/100‑1 (MQ=255)
                                                           tCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGAt                                   <  2:1222192/100‑1 (MQ=255)
                                                            cACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGAtt                                  <  1:734327/100‑1 (MQ=255)
                                                             aCGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTc                                 <  1:2532927/100‑1 (MQ=255)
                                                             aCGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTc                                 >  1:3945068/1‑100 (MQ=255)
                                                              cGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCg                                >  1:1559082/1‑100 (MQ=255)
                                                                ttCGCGGGCGACGGAAAGTGAAGTGACGATCATTACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCggg                              <  2:2650095/100‑1 (MQ=255)
                                                                 tCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGt                             <  1:835340/100‑1 (MQ=255)
                                                                  cgcgGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGtt                            >  1:1984724/1‑100 (MQ=255)
                                                                    cgGGCGACGGAAAGTGAAGTGACGATCATTACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTtat                          <  2:2424324/100‑1 (MQ=255)
                                                                     gggCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTtata                         <  1:3136606/100‑1 (MQ=255)
                                                                      ggCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATAc                        >  1:1578396/1‑100 (MQ=255)
                                                                       gCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATAcc                       <  1:1717277/100‑1 (MQ=255)
                                                                       gCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCAAAAACCATTTGTAGTCGAGATTCGGGTTATAcc                       >  2:4001493/1‑100 (MQ=255)
                                                                         gACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAg                     >  1:1984008/1‑100 (MQ=255)
                                                                           cGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGtt                   >  2:968352/1‑100 (MQ=255)
                                                                            ggAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTg                  <  2:2977735/100‑1 (MQ=255)
                                                                             gAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGc                 >  1:1194878/1‑100 (MQ=255)
                                                                             gAAAGTGAAGTCACGATCATAACGCCGATAGTCGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGc                 <  2:1548460/100‑1 (MQ=255)
                                                                              taagtgaagtgACGATCATGACGCCGATAGTGGTGATCATGGCGATCACTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCg                <  2:3888535/99‑1 (MQ=255)
                                                                              aaaGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCg                <  2:337083/100‑1 (MQ=255)
                                                                               aagtgaagtgACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGt               >  2:1662179/1‑100 (MQ=255)
                                                                                agtgaagtgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTa              >  2:3182011/1‑100 (MQ=255)
                                                                                agtgaagcgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGTACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTa              <  2:3716514/100‑1 (MQ=255)
                                                                                 gtgaagtgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTac             >  1:1882739/1‑100 (MQ=255)
                                                                                   gaagtgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTacca           <  2:2808773/100‑1 (MQ=255)
                                                                                   gaagtgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTacca           >  2:1000242/1‑100 (MQ=255)
                                                                                    aagtgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTaccac          <  2:2363309/100‑1 (MQ=255)
                                                                                    aagtgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTaccac          >  1:2776326/1‑100 (MQ=255)
                                                                                     agtgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTaccacc         <  2:3447214/100‑1 (MQ=255)
                                                                                       tgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTACCACCAg       >  1:3758551/1‑100 (MQ=255)
                                                                                       tgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTACCACCAg       <  1:3317425/100‑1 (MQ=255)
                                                                                        gaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTACCACCAGc      >  2:2880115/1‑100 (MQ=255)
                                                                                          cGATCATGACGCCGATAGTGGTTATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTCCCACCAGCTc    >  2:3477040/1‑100 (MQ=255)
                                                                                          cGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTCTAGTCGCGATTCGGGTTATACCAGTTGCGTACCACCAGCTc    <  2:191578/100‑1 (MQ=255)
                                                                                            aTCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTACCACCAGCTCGc  <  1:1249714/100‑1 (MQ=255)
                                                                                                 |                                                                                              
CGATGAAGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATTACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTACCACCAGCTCGC  >  NC_000913/825448‑825639

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