Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 801,422 G→A 65.5% A222T (GCA→ACA)  ybhH → putative PrpF family isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913801,4220GA65.5% 135.1 / 125.7 136A222T (GCA→ACA) ybhHputative PrpF family isomerase
Reads supporting (aligned to +/- strand):  ref base G (27/20);  new base A (47/42);  total (74/62)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

TGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTT  >  NC_000913/801324‑801521
                                                                                                  |                                                                                                   
 gACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGAGAAAAACAGCTTATGAATTACAGGCGGAACTGAACGCCGACAAAACa                                                                                                   >  2:3897910/1‑100 (MQ=255)
   ccTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACatt                                                                                                 >  1:1489349/1‑100 (MQ=255)
    cTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACatta                                                                                                <  1:2443601/100‑1 (MQ=255)
     tGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACattat                                                                                               <  1:151076/100‑1 (MQ=255)
       tATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCAttatta                                                                                             <  1:950853/100‑1 (MQ=255)
       tATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTAAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCAttatta                                                                                             >  2:2212737/1‑100 (MQ=255)
        aTCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAg                                                                                            >  2:1413694/1‑100 (MQ=255)
        aTCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAg                                                                                            <  1:2851795/100‑1 (MQ=255)
         tCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGc                                                                                           <  1:866925/100‑1 (MQ=255)
         tCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGc                                                                                           <  1:3504947/100‑1 (MQ=255)
          cGATATGGCGAGGCGAGTCGTAATTAGTCCGGATGAATATCTGGGAAAAAGAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGcc                                                                                          <  2:4044199/100‑1 (MQ=255)
           gATATGGCGATGCCAGTCGTCATTATTACGGCGGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGcc                                                                                          <  2:1617972/99‑1 (MQ=255)
            atatGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCg                                                                                        <  1:1347055/100‑1 (MQ=255)
             tatGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGc                                                                                       >  1:3116511/1‑100 (MQ=255)
               tGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCAt                                                                                     >  1:1346557/1‑100 (MQ=255)
                   gATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGaa                                                                                 <  1:1590800/100‑1 (MQ=255)
                   gATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGaa                                                                                 <  2:917236/100‑1 (MQ=255)
                    aTGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAAt                                                                                <  2:183628/100‑1 (MQ=255)
                    aTGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAAt                                                                                <  1:3728650/100‑1 (MQ=255)
                     tGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGACCTGGATGCCGACCAAACATTATTAGCCCGCATTGACTc                                                                               >  1:3360260/1‑100 (MQ=255)
                     tGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATc                                                                               <  1:3452718/100‑1 (MQ=255)
                      gCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATc                                                                               <  2:2602257/99‑1 (MQ=255)
                       ccAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTa                                                                             >  2:3261620/1‑100 (MQ=255)
                        cAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTAt                                                                            >  2:2688075/1‑100 (MQ=255)
                        cAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTAt                                                                            >  1:2296676/1‑100 (MQ=255)
                         aGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATc                                                                           >  2:1151890/1‑100 (MQ=255)
                            cgtcATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGt                                                                        >  2:2229325/1‑100 (MQ=255)
                            cgtcATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGt                                                                        >  2:2230620/1‑100 (MQ=255)
                            cgtcATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGt                                                                        >  1:3130975/1‑100 (MQ=255)
                                 ttattCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACa                                                                   >  1:30642/1‑100 (MQ=255)
                                 ttattCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACa                                                                   >  1:30643/1‑100 (MQ=255)
                                    ttCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGc                                                                >  2:2129329/1‑100 (MQ=255)
                                    ttCCGGCTGAATATCTGGGAAAAACAGGTTATGAATGACCGGCGGAACGGGATGCCTACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGc                                                                <  1:823884/100‑1 (MQ=255)
                                     tCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCg                                                               <  1:4056376/100‑1 (MQ=255)
                                      tcGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCg                                                               <  2:2041221/98‑1 (MQ=255)
                                       cGGCTGAATAGCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCggg                                                             <  1:1956489/100‑1 (MQ=255)
                                        ggCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGt                                                            >  1:3122517/1‑100 (MQ=255)
                                         gCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCgggca                                                           >  2:2504878/1‑98 (MQ=255)
                                           tGAATATATGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTaaa                                                         >  1:2173903/1‑100 (MQ=255)
                                           tGAATATATGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGGCCGCATTGAATCTATCCGTCTACAAGCGGGTaaa                                                         >  2:2169517/1‑100 (MQ=255)
                                            gAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAg                                                        <  1:1278302/100‑1 (MQ=255)
                                            gAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAg                                                        >  1:1867681/1‑100 (MQ=255)
                                             aaTATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGc                                                       >  1:3748253/1‑100 (MQ=255)
                                             aaTATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGc                                                       <  1:3218783/100‑1 (MQ=255)
                                             aaTATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGc                                                       <  1:1579708/100‑1 (MQ=255)
                                              atatCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGc                                                       <  2:2398933/99‑1 (MQ=255)
                                               tatCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCaa                                                     <  2:2187834/100‑1 (MQ=255)
                                                   tGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATggg                                                 >  1:539945/1‑100 (MQ=255)
                                                   tGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATggg                                                 <  1:153397/100‑1 (MQ=255)
                                                   tGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATggg                                                 >  1:383871/1‑100 (MQ=255)
                                                    gggAAAAAGAGGTTATGCAGTACCGGCGGAACTGGATGGTGACAAAACATTAGTAGAGCGGAATGAGTCTATCCGGCTACAAGCGGGTAAAGCAATgggg                                                <  2:3556965/100‑2 (MQ=255)
                                                    gggAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGc                                                <  2:3464720/100‑1 (MQ=255)
                                                    gggAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGc                                                >  2:2924650/1‑100 (MQ=255)
                                                     ggAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACCAAACATTAGTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCt                                               <  2:505569/100‑1 (MQ=255)
                                                     ggAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCGCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGc                                                >  2:3930427/1‑99 (MQ=255)
                                                      gAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCtt                                              >  1:1415580/1‑100 (MQ=255)
                                                       aaaaaCAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTa                                             <  1:2473148/100‑1 (MQ=255)
                                                        aaaaCAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAg                                            >  2:832574/1‑100 (MQ=255)
                                                         aaaCAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAgg                                           <  2:2091576/100‑1 (MQ=255)
                                                           aCAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTg                                         >  2:4017334/1‑100 (MQ=255)
                                                            cAGGTTATGAATTACCGGCGGAACTGGATGGCGCCAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGCGCTTAGGTGa                                        >  2:2719073/1‑100 (MQ=255)
                                                            cAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGa                                        <  2:4102700/100‑1 (MQ=255)
                                                            cAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGa                                        <  2:4105131/100‑1 (MQ=255)
                                                             aGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGAt                                       <  2:442774/100‑1 (MQ=255)
                                                              ggTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATg                                      <  1:2440943/100‑1 (MQ=255)
                                                              ggTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATg                                      <  2:1420561/100‑1 (MQ=255)
                                                               gTTATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGt                                     >  2:3238446/1‑100 (MQ=255)
                                                                ttATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTc                                    <  2:3056668/100‑1 (MQ=255)
                                                                ttATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTc                                    <  2:3703564/100‑1 (MQ=255)
                                                                ttATGAATTACCGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTc                                    <  1:152413/100‑1 (MQ=255)
                                                                    gAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTa                                <  1:4017360/100‑1 (MQ=255)
                                                                          ccGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATgg                          <  1:2203683/100‑1 (MQ=255)
                                                                          ccGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATgg                          <  1:4086652/100‑1 (MQ=255)
                                                                          ccGGCGGAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATgg                          >  2:1556293/1‑100 (MQ=255)
                                                                           cggcggAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGGCAGTAATATgg                          >  1:2248799/1‑99 (MQ=255)
                                                                            ggcggAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGtt                        >  2:2007961/1‑100 (MQ=255)
                                                                              cggAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTAt                      >  2:1078754/1‑100 (MQ=255)
                                                                               ggAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTGCAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATc                     >  2:3427274/1‑100 (MQ=255)
                                                                                gAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATcc                    >  1:488896/1‑100 (MQ=255)
                                                                                gAACTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATcc                    <  2:2373011/100‑1 (MQ=255)
                                                                                 aaCTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATccc                   <  2:3122498/100‑1 (MQ=255)
                                                                                  aCTGGATGCCGGCAAAACATTCTTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTCGGTGATGTCAGTAATATGGTTATccat                  >  2:2943036/1‑98 (MQ=255)
                                                                                  aCTGGATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCt                  >  1:1379899/1‑100 (MQ=255)
                                                                                     ggATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTaaa               >  1:1167740/1‑100 (MQ=255)
                                                                                      gATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAAc              >  2:245352/1‑100 (MQ=255)
                                                                                      gATGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAAc              >  2:1149686/1‑100 (MQ=255)
                                                                                       aTGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAAcc             >  2:289816/1‑100 (MQ=255)
                                                                                        tGCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCt            >  2:1715118/1‑100 (MQ=255)
                                                                                         gCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCtg           >  1:2687583/1‑100 (MQ=255)
                                                                                         gCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCtg           >  1:2312800/1‑100 (MQ=255)
                                                                                         gCCGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCATAAAcc             <  2:2512020/98‑1 (MQ=255)
                                                                                          ccGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCtgt          >  2:637338/1‑100 (MQ=255)
                                                                                          ccGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCtgt          >  1:1513786/1‑100 (MQ=255)
                                                                                          ccGACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCtgt          <  1:121092/100‑1 (MQ=255)
                                                                                           cGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCtg           <  2:141829/98‑1 (MQ=255)
                                                                                           cGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGCCAGTAATATGGTTATCCCTAAACCtgt          >  2:2991578/1‑99 (MQ=255)
                                                                                            gACAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGc        <  1:2688091/100‑1 (MQ=255)
                                                                                             aCAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGACTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCt       >  2:609416/1‑100 (MQ=255)
                                                                                              cAAAACATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCtt      >  1:1257048/1‑100 (MQ=255)
                                                                                                  gCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTAttt  >  1:3445690/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TGACCTGTATCGATATGGCGATGCCAGTCGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTT  >  NC_000913/801324‑801521

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