Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,978 A→C 25.0% intergenic (‑188/+120) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9780AC25.0% 515.5 / 167.6 243intergenic (‑188/+120)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (69/101);  new base C (47/2);  total (140/103)
Fisher's exact test for biased strand distribution p-value = 2.88e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-03

Only 100 of 247 total aligned reads displayed.

AAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705880‑3706077
                                                                                                  |                                                                                                           
  gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                          <  1:2897167/100‑1 (MQ=255)
  gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                          <  1:280471/100‑1 (MQ=255)
   tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                         <  2:1559957/100‑1 (MQ=255)
    aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                        >  2:1773937/1‑100 (MQ=255)
    aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                        <  1:2523573/100‑1 (MQ=255)
     ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                       <  2:2869453/100‑1 (MQ=255)
      tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                      >  1:1022670/1‑100 (MQ=255)
       gctgcCGGATGACGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                     <  1:1093019/100‑1 (MQ=255)
        ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                    >  2:3712696/1‑100 (MQ=255)
            cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                >  1:803558/1‑100 (MQ=255)
             ggATGGCGCTACACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTCTACGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                               <  1:2461596/100‑1 (MQ=255)
              gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                              <  1:727678/100‑2 (MQ=37)
               aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                             <  2:382467/100‑1 (MQ=255)
                   cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                        <  2:1804473/100‑1 (MQ=255)
                      tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                                     >  1:2976622/1‑100 (MQ=255)
                       gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGTGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                    >  2:1201301/1‑100 (MQ=255)
                        cACGCAGAGGCACAGCGCCACCGGGTATACATATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                   <  2:505950/100‑1 (MQ=255)
                         aCGCAGCGCCACAGCGCCACCGGGTCTACAGATTGTTGACAAAGCGCGCATCCGTTTCACCCGATGCGACGTTAACGCCTTATCCCGCCTCCAACATc                                                                                    >  1:1735371/1‑98 (MQ=16)
                          cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                 <  2:1893263/100‑1 (MQ=255)
                           gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                <  2:2958851/100‑1 (MQ=255)
                            cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGTCCTACAAAATCGTGCa                                                                               >  1:3459258/1‑100 (MQ=255)
                               agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                            <  1:2733075/100‑1 (MQ=255)
                                   acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                        <  2:3577927/100‑1 (MQ=255)
                                    caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                       <  2:2370798/100‑1 (MQ=255)
                                    caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                       <  2:401821/100‑1 (MQ=255)
                                        gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                                   >  1:1187763/1‑100 (MQ=255)
                                         ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                  >  1:112978/1‑100 (MQ=255)
                                           cgatcttataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                <  2:1431508/94‑1 (MQ=255)
                                             cGGGTATACAGATTGTTGACAAAGTGCTCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑Ag                                                           <  1:998186/100‑1 (MQ=255)
                                             cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑Ag                                                           >  2:1520647/1‑100 (MQ=255)
                                             cGGGTATACAGATTGTTGACAAAGTGCGCTTGGTTTATGCCGGATGAGGAGTTAACGCCTTAGCCATCCTACAAAATCGTGCGAATTGAATACATTGc                                                                >  1:3643028/1‑98 (MQ=17)
                                                 tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtt                                                     >  2:3153657/1‑100 (MQ=255)
                                                     cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑g                                               <  2:646015/100‑1 (MQ=255)
                                                     cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑g                                               <  2:1187356/100‑1 (MQ=255)
                                                     cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑g                                               <  1:2832051/100‑1 (MQ=255)
                                                      aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑Gt                                              >  2:3621596/1‑100 (MQ=255)
                                                       gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTa                                             <  1:1220134/100‑1 (MQ=255)
                                                       gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTa                                             <  2:2976643/100‑1 (MQ=255)
                                                        aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAg                                            <  1:2373010/100‑1 (MQ=255)
                                                          tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGc                                          >  2:2109402/1‑100 (MQ=255)
                                                             tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTg                                       >  2:282792/1‑100 (MQ=255)
                                                              gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGa                                      >  1:2021679/1‑100 (MQ=255)
                                                                 aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATaa                                   <  1:848009/100‑1 (MQ=255)
                                                                   aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGc                                 <  1:1661067/100‑1 (MQ=255)
                                                                     tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCgg                               >  1:3081429/1‑100 (MQ=255)
                                                                      gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGCCCTGATAAGCGGa                              <  2:2538387/100‑1 (MQ=255)
                                                                       cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAg                             <  1:3000541/100‑1 (MQ=255)
                                                                        gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                            <  1:1773930/100‑1 (MQ=255)
                                                                         cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                           >  1:372536/1‑100 (MQ=255)
                                                                           ttGTGTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAGAGTGCTAATTAAATACATTGC‑‑TGTATTTTC‑‑GTAGGCCTGATAAGCgg                               <  1:440247/94‑1 (MQ=17)
                                                                             gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                       <  1:1275269/100‑1 (MQ=255)
                                                                              tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                      <  2:3181627/100‑1 (MQ=255)
                                                                              tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                      <  1:1201295/100‑1 (MQ=255)
                                                                              tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                      >  2:500551/1‑100 (MQ=255)
                                                                              tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  1:3605245/6‑100 (MQ=18)
                                                                              tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                      >  1:961340/6‑100 (MQ=17)
                                                                              tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  1:1468470/6‑100 (MQ=21)
                                                                               ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                     <  1:1171033/100‑1 (MQ=255)
                                                                               gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGC‑‑AGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     <  1:2937904/96‑1 (MQ=17)
                                                                               gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAAACCTGCAAATTCAATATATTGC‑‑AGGATACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  1:494504/5‑100 (MQ=18)
                                                                                 aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑ATAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  2:3550985/3‑100 (MQ=17)
                                                                                 aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  1:2568609/3‑100 (MQ=17)
                                                                                 aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  1:3307256/3‑100 (MQ=17)
                                                                                 aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  2:3500583/3‑100 (MQ=18)
                                                                                 aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAGAAAAT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                  >  1:1938436/3‑99 (MQ=21)
                                                                                  ggtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAACGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  1:1978380/11‑100 (MQ=21)
                                                                                  ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:1296614/2‑100 (MQ=17)
                                                                                  ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:2468600/2‑100 (MQ=17)
                                                                                  ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑AAGATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:2696508/2‑100 (MQ=17)
                                                                                  gaccggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTTCCAAATTCAATATATTGC‑‑AGAAATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:3456266/11‑100 (MQ=18)
                                                                                   gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                 <  2:3267846/100‑1 (MQ=255)
                                                                                   gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                 >  2:17402/1‑100 (MQ=255)
                                                                                   gCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTGC‑‑AGAGATTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  2:2427896/1‑99 (MQ=14)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGC‑‑AATATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg                  >  2:2232253/1‑97 (MQ=21)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  2:1675406/1‑99 (MQ=21)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  2:1096657/1‑99 (MQ=21)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag                 >  2:2581391/1‑99 (MQ=18)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGC‑‑AGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  1:2602598/1‑99 (MQ=17)
                                                                                    ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  2:746243/9‑98 (MQ=17)
                                                                                    ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  2:2291139/9‑98 (MQ=18)
                                                                                    ccggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                >  1:2325546/9‑98 (MQ=25)
                                                                                    ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑AAGATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                >  1:318330/9‑98 (MQ=17)
                                                                                    ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                >  2:1962587/1‑100 (MQ=255)
                                                                                    ccGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGC‑‑AGAAATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                >  1:1059579/1‑98 (MQ=17)
                                                                                    ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                >  2:87263/1‑98 (MQ=18)
                                                                                     cggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt               >  2:3379908/8‑97 (MQ=18)
                                                                                      ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt              >  2:1349059/7‑100 (MQ=18)
                                                                                      ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtat              >  1:2271363/7‑96 (MQ=17)
                                                                                      ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt              <  1:1389415/100‑1 (MQ=255)
                                                                                        aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg            <  1:1847113/100‑1 (MQ=255)
                                                                                        aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg            <  2:774892/100‑1 (MQ=255)
                                                                                         tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc           >  1:3397427/1‑100 (MQ=255)
                                                                                           cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct         >  2:3567587/1‑100 (MQ=255)
                                                                                             gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt       <  2:1649900/100‑1 (MQ=255)
                                                                                              cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      <  2:2026018/100‑1 (MQ=255)
                                                                                              cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      <  2:2004902/100‑1 (MQ=255)
                                                                                              cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      <  1:2109391/100‑1 (MQ=255)
                                                                                              cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAGTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      <  1:2169642/100‑1 (MQ=255)
                                                                                                 tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa   <  2:1385288/100‑1 (MQ=255)
                                                                                                  aaCGCCTTATCCGGCCTACAAAATCGTGCAGATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt  >  2:2664900/1‑100 (MQ=255)
                                                                                                  |                                                                                                           
AAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705880‑3706077

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