Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,978 A→C 23.2% 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,9780AC23.2% 323.2 / 93.5 154intergenic (‑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 (52/60);  new base C (29/2);  total (92/62)
Fisher's exact test for biased strand distribution p-value = 9.33e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-02

Only 100 of 158 total aligned reads displayed.

GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705879‑3706077
                                                                                                   |                                                                                                        
gAAGTATTGCTTCCGGATGGCGCTGCACGCAGAGGCCCAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGAATGCGGCGTTa                                                                                                          >  2:1786103/1‑100 (MQ=255)
 aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                         >  2:773845/1‑100 (MQ=255)
     aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                     <  2:250664/100‑1 (MQ=255)
      ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                    <  1:1533050/100‑1 (MQ=255)
      ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                    <  1:2312421/100‑1 (MQ=255)
      ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                    >  1:2293182/1‑100 (MQ=255)
      ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                    <  1:2611395/100‑1 (MQ=255)
      ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                    <  1:1707229/100‑1 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                  >  1:2994121/1‑100 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                  <  2:697540/100‑1 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                               <  2:1827955/100‑1 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                               >  1:1381441/1‑100 (MQ=255)
              ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCAGGc                                                                                            >  1:670220/1‑100 (MQ=255)
              ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGACGTTAACGCCTTATcc                                                                                               >  1:3015248/1‑97 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                           <  1:2701208/100‑1 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                           <  2:19662/100‑1 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGCTAACGCCTTATCCGGcc                                                                                           >  2:75192/1‑100 (MQ=255)
                aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                          <  1:797027/100‑3 (MQ=37)
                aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                          <  1:807423/100‑3 (MQ=37)
                 tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                         <  2:139529/100‑1 (MQ=255)
                   gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGAATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                       >  1:1012399/1‑100 (MQ=255)
                    cgcTGCACGCAGAGGCACAGCGCCACCGGTTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                      <  1:257249/100‑1 (MQ=255)
                    cgcTGCACGCAGAGGCACAGCACCACCGGGTATACAGATTGTTGACAAACTGCGCTTTGTTTATGCCGGATGCCGCGTTAACGCCTTATCCGTCCTACaa                                                                                      >  2:2615250/1‑100 (MQ=25)
                        gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaag                                                                                  <  2:2842761/100‑11 (MQ=21)
                         cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAAGCCTTATCCGGCCTACAAAATCg                                                                                 >  1:468003/1‑100 (MQ=255)
                          aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                >  1:1186042/1‑100 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                >  1:1415153/1‑99 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                               >  2:1958097/1‑100 (MQ=255)
                               gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                           >  2:57883/1‑100 (MQ=255)
                                agGCACAGCGCCACCGGGTATATAGTTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                          <  1:2117489/100‑1 (MQ=255)
                                  gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                        >  2:1190205/1‑100 (MQ=255)
                                     caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                     <  2:1050109/100‑1 (MQ=255)
                                      aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                     <  2:326091/99‑1 (MQ=255)
                                       gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                   <  2:45238/100‑1 (MQ=255)
                                         gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCAGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                                 >  1:1493294/1‑100 (MQ=255)
                                             ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                             <  2:1364976/100‑1 (MQ=255)
                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑g                                                           <  2:499872/100‑1 (MQ=255)
                                                ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTACAAAATCGTGCAAATTCAATACATTGCA‑Gaa                                                         >  1:1141105/1‑100 (MQ=255)
                                                    taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAtt‑‑tt                                                   >  1:986737/1‑100 (MQ=255)
                                                    taCAGATTGTTGACAAAGTGCGCTTTGTTTATCACGGATGCAGCGTTAACTCCTTATCCGGCCTACAAAATCGTGCAAATTCATTACATTGCA‑GAAtt‑‑tt                                                   <  2:85273/100‑1 (MQ=17)
                                                    taCAGATTGTTGAAAAAGTGCGCTTTGTTTATGCCGGACGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAtt‑‑tt                                                   >  2:1871446/1‑100 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑g                                               >  1:2857584/1‑100 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑g                                               <  2:2293156/100‑1 (MQ=255)
                                                       aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑Gt                                              >  1:701618/1‑100 (MQ=255)
                                                        gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCATATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTa                                             <  1:1843893/100‑1 (MQ=255)
                                                            gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGTTTTATCCTGTTTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGcc                                         <  1:983872/100‑1 (MQ=18)
                                                             ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTATGCCt                                        >  2:1931382/1‑100 (MQ=255)
                                                               gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGa                                      >  2:1712624/1‑100 (MQ=255)
                                                                aCAAAGTGCGCTTTGTTTATGCAGGATGCGGCGTTAAAGCCTTATCAGGCATACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGAt                                     >  2:2140109/1‑100 (MQ=18)
                                                                    aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGc                                 >  1:2827266/1‑100 (MQ=255)
                                                                      tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCgg                               <  1:330353/100‑1 (MQ=255)
                                                                       gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATCAGCGGa                              >  2:1217667/1‑100 (MQ=255)
                                                                        cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                             >  1:587145/1‑100 (MQ=255)
                                                                        cgcTTTGTTGATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                             <  2:2570162/100‑1 (MQ=255)
                                                                         gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                            >  2:373331/1‑100 (MQ=255)
                                                                          cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                            <  1:405163/99‑1 (MQ=255)
                                                                             tttgtaggCCGGATAAGTCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTACA‑GAGAT‑‑TGC‑‑GTAGGCCTGATAAGCGTAg                             >  2:1462092/8‑95 (MQ=11)
                                                                             tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCA‑ATTCT‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                        >  1:2836239/8‑100 (MQ=18)
                                                                             tGTTTCTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAAT‑‑TTC‑‑GTACGCCTGATAAGCGGAgcgc                          >  2:1938451/1‑98 (MQ=255)
                                                                             tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAATGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                        >  2:1671444/1‑100 (MQ=255)
                                                                              ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AATC‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATc                       >  2:2399942/7‑100 (MQ=18)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCTTGCAAATTCAATACATTGCA‑GAAAT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                       >  1:1652348/1‑100 (MQ=255)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GGATT‑‑TTC‑‑GTAGGCCTGATCAGCGGAGCGCATc                       <  1:2044564/100‑1 (MQ=255)
                                                                               tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCA‑TGGGC‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                      >  1:1815373/6‑100 (MQ=17)
                                                                                gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑GGAGC‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  2:3020255/5‑100 (MQ=18)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGAT‑CGATC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  2:3026071/4‑100 (MQ=21)
                                                                                   ggccagataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AATC‑‑TCCGTGTAGGCCTGACAAGCGTAGCGCATCAGGCa                  >  2:2877078/11‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑GAAAC‑‑CCT‑‑GTAGGCCGGATAAGCGTAGCGCATCAGGCa                  >  2:2050239/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                  >  1:509184/2‑100 (MQ=21)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGAT‑CGATC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:1106044/2‑100 (MQ=21)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCA‑GGGAC‑‑TAT‑‑GTAGACCTGATAAGCATAGCGCATCAGGCa                  >  1:1212855/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCA‑GAGAC‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                  >  1:401355/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCA‑GGATA‑‑CGT‑‑ATAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  1:321438/2‑100 (MQ=17)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                 >  1:370551/1‑100 (MQ=255)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                 >  1:1233024/1‑100 (MQ=255)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACACAAACATGCAAATTCAATAGATTGCA‑GAGAT‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  1:25540/1‑99 (MQ=17)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCA‑‑ATAT‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat                >  1:494082/1‑97 (MQ=21)
                                                                                     ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCA‑‑ATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                 >  1:449408/9‑99 (MQ=21)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTG‑GGCTC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  2:945169/9‑98 (MQ=21)
                                                                                     ccGGATCCGGCGTTAACGCCTTATCCAGCCTACAAAATTGTGCAAAATCAATTGATTGCA‑CAGCT‑‑AGC‑‑GTAGGCCTGATGCGCTACGCTCATCAGGc                   <  2:1483766/97‑1 (MQ=18)
                                                                                      cggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCA‑CCATC‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                >  2:1357700/8‑97 (MQ=18)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑GAAAC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  1:1458984/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCGTCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GGTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct              >  1:2580452/7‑96 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCA‑TGAGC‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  2:1559778/7‑100 (MQ=14)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct              >  2:756655/7‑96 (MQ=25)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑GGAGC‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt               >  1:3034549/7‑99 (MQ=17)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAGT‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat             >  1:1066030/7‑98 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACG‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta              >  2:2269240/7‑96 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCATAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAGA‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACAttta              >  1:2756742/7‑98 (MQ=18)
                                                                                        gataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑GAGTT‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt             >  2:2192871/6‑100 (MQ=17)
                                                                                         aTGCGGGGTTAACTCCTTATCCGGCATACAAAATGGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTTATAAGCGGAGTGCATCAGGGATTtttg            <  2:1270139/100‑1 (MQ=18)
                                                                                           gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc          >  1:895262/1‑100 (MQ=255)
                                                                                            cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct         <  1:2913962/100‑1 (MQ=255)
                                                                                            cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct         >  1:2748791/1‑100 (MQ=255)
                                                                                             ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt        >  1:1614734/1‑100 (MQ=255)
                                                                                               cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      >  2:636789/1‑100 (MQ=255)
                                                                                                 ttAACGCCTTATCCGGCCTACAAAATCGTGTAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      <  2:2555109/98‑1 (MQ=255)
                                                                                                 ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc    <  2:475307/100‑1 (MQ=255)
                                                                                                  tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc    >  2:1710077/1‑99 (MQ=255)
                                                                                                   aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt  >  1:2196911/1‑100 (MQ=255)
                                                                                                   |                                                                                                        
GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705879‑3706077

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