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.

GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑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)
                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                           <  2:499872/100‑1 (MQ=255)
                                                ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                         >  1:1141105/1‑100 (MQ=255)
                                                    taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑tt                                                   >  1:986737/1‑100 (MQ=255)
                                                    taCAGATTGTTGACAAAGTGCGCTTTGTTTATCACGGATGCAGCGTTAACTCCTTATCCGGCCTACAAAATCGTGCAAATTCATTACATTGCAGAAtt‑‑tt                                                   <  2:85273/100‑1 (MQ=17)
                                                    taCAGATTGTTGAAAAAGTGCGCTTTGTTTATGCCGGACGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑tt                                                   >  2:1871446/1‑100 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑g                                               >  1:2857584/1‑100 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑g                                               <  2:2293156/100‑1 (MQ=255)
                                                       aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑Gt                                              >  1:701618/1‑100 (MQ=255)
                                                        gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCATATTCAATACATTGCAGAATT‑‑TTC‑‑GTa                                             <  1:1843893/100‑1 (MQ=255)
                                                            gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGTTTTATCCTGTTTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGcc                                         <  1:983872/100‑1 (MQ=18)
                                                             ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTATGCCt                                        >  2:1931382/1‑100 (MQ=255)
                                                               gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGa                                      >  2:1712624/1‑100 (MQ=255)
                                                                aCAAAGTGCGCTTTGTTTATGCAGGATGCGGCGTTAAAGCCTTATCAGGCATACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGAt                                     >  2:2140109/1‑100 (MQ=18)
                                                                    aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGc                                 >  1:2827266/1‑100 (MQ=255)
                                                                      tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCgg                               <  1:330353/100‑1 (MQ=255)
                                                                       gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATCAGCGGa                              >  2:1217667/1‑100 (MQ=255)
                                                                        cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                             >  1:587145/1‑100 (MQ=255)
                                                                        cgcTTTGTTGATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                             <  2:2570162/100‑1 (MQ=255)
                                                                         gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                            >  2:373331/1‑100 (MQ=255)
                                                                          cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                            <  1:405163/99‑1 (MQ=255)
                                                                             tttgtaggCCGGATAAGTCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTACAGAGAT‑‑TGC‑‑GTAGGCCTGATAAGCGTAg                             >  2:1462092/8‑95 (MQ=11)
                                                                             tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                        >  1:2836239/8‑100 (MQ=18)
                                                                             tGTTTCTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAAT‑‑TTC‑‑GTACGCCTGATAAGCGGAgcgc                          >  2:1938451/1‑98 (MQ=255)
                                                                             tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAATGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                        >  2:1671444/1‑100 (MQ=255)
                                                                              ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATC‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATc                       >  2:2399942/7‑100 (MQ=18)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCTTGCAAATTCAATACATTGCAGAAAT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                       >  1:1652348/1‑100 (MQ=255)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGGATT‑‑TTC‑‑GTAGGCCTGATCAGCGGAGCGCATc                       <  1:2044564/100‑1 (MQ=255)
                                                                               tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGC‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                      >  1:1815373/6‑100 (MQ=17)
                                                                                gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  2:3020255/5‑100 (MQ=18)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  2:3026071/4‑100 (MQ=21)
                                                                                   ggccagataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATC‑‑TCCGTGTAGGCCTGACAAGCGTAGCGCATCAGGCa                  >  2:2877078/11‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑CCT‑‑GTAGGCCGGATAAGCGTAGCGCATCAGGCa                  >  2:2050239/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                  >  1:509184/2‑100 (MQ=21)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:1106044/2‑100 (MQ=21)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑‑TAT‑‑GTAGACCTGATAAGCATAGCGCATCAGGCa                  >  1:1212855/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGAGAC‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                  >  1:401355/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATA‑‑CGT‑‑ATAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  1:321438/2‑100 (MQ=17)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                 >  1:370551/1‑100 (MQ=255)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                 >  1:1233024/1‑100 (MQ=255)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACACAAACATGCAAATTCAATAGATTGCAGAGAT‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  1:25540/1‑99 (MQ=17)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat                >  1:494082/1‑97 (MQ=21)
                                                                                     ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                 >  1:449408/9‑99 (MQ=21)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  2:945169/9‑98 (MQ=21)
                                                                                     ccGGATCCGGCGTTAACGCCTTATCCAGCCTACAAAATTGTGCAAAATCAATTGATTGCACAGCT‑‑AGC‑‑GTAGGCCTGATGCGCTACGCTCATCAGGc                   <  2:1483766/97‑1 (MQ=18)
                                                                                      cggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                >  2:1357700/8‑97 (MQ=18)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  1:1458984/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCGTCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGGTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct              >  1:2580452/7‑96 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  2:1559778/7‑100 (MQ=14)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct              >  2:756655/7‑96 (MQ=25)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt               >  1:3034549/7‑99 (MQ=17)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat             >  1:1066030/7‑98 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACG‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta              >  2:2269240/7‑96 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCATACAAATCATGCAAATCCAAAATATTGCAGGAGA‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACAttta              >  1:2756742/7‑98 (MQ=18)
                                                                                        gataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt             >  2:2192871/6‑100 (MQ=17)
                                                                                         aTGCGGGGTTAACTCCTTATCCGGCATACAAAATGGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTTATAAGCGGAGTGCATCAGGGATTtttg            <  2:1270139/100‑1 (MQ=18)
                                                                                           gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc          >  1:895262/1‑100 (MQ=255)
                                                                                            cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct         <  1:2913962/100‑1 (MQ=255)
                                                                                            cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct         >  1:2748791/1‑100 (MQ=255)
                                                                                             ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt        >  1:1614734/1‑100 (MQ=255)
                                                                                               cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      >  2:636789/1‑100 (MQ=255)
                                                                                                 ttAACGCCTTATCCGGCCTACAAAATCGTGTAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      <  2:2555109/98‑1 (MQ=255)
                                                                                                 ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc    <  2:475307/100‑1 (MQ=255)
                                                                                                  tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc    >  2:1710077/1‑99 (MQ=255)
                                                                                                   aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt  >  1:2196911/1‑100 (MQ=255)
                                                                                                   |                                                                                                        
GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑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: