Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 31.4% intergenic (‑195/+113) 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,9850TC31.4% 218.6 / 131.6 153intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (55/50);  new base C (40/8);  total (95/58)
Fisher's exact test for biased strand distribution p-value = 2.78e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-03

Only 100 of 153 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGT  >  NC_000913/3705887‑3706081
                                                                                                  |                                                                                                    
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                     >  1:2994121/1‑100 (MQ=255)
   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                  >  1:1381441/1‑100 (MQ=255)
   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                  <  2:1827955/100‑1 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCAGGc                                                                                               >  1:670220/1‑100 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                               <  2:1736085/100‑1 (MQ=37)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                              <  2:19662/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                              <  1:2701208/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGCTAACGCCTTATCCGGcc                                                                                              >  2:75192/1‑100 (MQ=255)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                             <  1:807423/100‑3 (MQ=37)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                             <  1:300444/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)
            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                         >  1:1569909/1‑100 (MQ=255)
            cgcTGCACGCAGAGGCACAGCACCACCGGGTATACAGATTGTTGACAAACTGCGCTTTGTTTATGCCGGATGCCGCGTTAACGCCTTATCCGTCCTACaa                                                                                         >  2:2615250/1‑100 (MQ=25)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAAGCCTTATCCGGCCTACAAAATCg                                                                                    >  1:468003/1‑100 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                  <  2:306123/100‑1 (MQ=255)
                       gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                              >  2:57883/1‑100 (MQ=255)
                        agGCACAGCGCCACCGGGTATATAGTTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                             <  1:2117489/100‑1 (MQ=255)
                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                             <  2:2172688/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                                                                      <  1:46985/100‑1 (MQ=255)
                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                      <  2:45238/100‑1 (MQ=255)
                                 gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAGAACGTGCAAATTCaacaaat                                                                    >  2:1919836/1‑95 (MQ=17)
                                     ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                                <  2:1364976/100‑1 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑g                                                              <  2:499872/100‑1 (MQ=255)
                                        ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTACAAAATCGTGCAAATTCAATACATTGCA‑Gaa                                                            >  1:1141105/1‑100 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCTGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTTCAAATTCAATACATTGCA‑GAAt‑ttt                                                       >  2:117613/1‑100 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATCACGGATGCAGCGTTAACTCCTTATCCGGCCTACAAAATCGTGCAAATTCATTACATTGCA‑GAAt‑ttt                                                       <  2:85273/100‑1 (MQ=17)
                                            taCAGATTGTTGAAAAAGTGCGCTTTGTTTATGCCGGACGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAt‑ttt                                                       >  2:1871446/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑g                                                   >  1:2857584/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑g                                                   >  1:1702126/1‑100 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCATATTCAATACATTGCA‑GAAT‑TTTC‑‑GTa                                                 <  1:1843893/100‑1 (MQ=255)
                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCTGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAgg                                               <  1:1201207/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGTTTTATCCTGTTTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGcc                                             <  1:983872/100‑1 (MQ=18)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTg                                           >  2:759821/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGa                                          >  2:1712624/1‑100 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGAt                                         <  2:2811126/100‑1 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCAGGATGCGGCGTTAAAGCCTTATCAGGCATACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGAt                                         >  2:2140109/1‑100 (MQ=18)
                                                           aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAg                                      <  1:2655613/100‑1 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGc                                     >  1:2827266/1‑100 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCgg                                   <  1:330353/100‑1 (MQ=255)
                                                                cgcTTTGTTGATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                                 <  2:2570162/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                >  2:373331/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAGAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                               <  2:1944875/100‑1 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                <  1:405163/99‑1 (MQ=255)
                                                                   tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcgc                              >  1:746255/1‑100 (MQ=255)
                                                                     tttgtaggCCGGATAAGTCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTACA‑GAGA‑TTGC‑‑GTAGGCCTGATAAGCGTAg                                 >  2:1462092/8‑95 (MQ=11)
                                                                     tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCA‑ATTC‑TCTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                            >  1:2836239/8‑100 (MQ=18)
                                                                     tGTTTCTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAA‑TTTC‑‑GTACGCCTGATAAGCGGAgcgc                              >  2:1938451/1‑98 (MQ=255)
                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAATGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                            >  2:1671444/1‑100 (MQ=255)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATc                           >  2:2399942/7‑100 (MQ=18)
                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                           <  2:2326732/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGAATCa                          >  1:577411/1‑100 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGACTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TATC‑‑GTAGGCCTGATAAACGGAGCGCATCa                          >  2:2268461/1‑100 (MQ=255)
                                                                       tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCA‑ATTC‑TCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                          >  1:1798739/6‑100 (MQ=18)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑GGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                         >  2:3020255/5‑100 (MQ=18)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAg                         <  1:2399970/96‑1 (MQ=18)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGAT‑CGAT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                        >  2:3026071/4‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCA‑TGGG‑CCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                       >  1:412280/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCA‑AGAT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                       >  1:396157/3‑100 (MQ=17)
                                                                           ggccagataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑CTCCGTGTAGGCCTGACAAGCGTAGCGCATCAGGCa                      >  2:2877078/11‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                      >  1:85464/2‑100 (MQ=14)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                      >  1:509184/2‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCA‑GGGA‑CTAT‑‑GTAGACCTGATAAGCATAGCGCATCAGGCa                      >  1:1212855/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCA‑GGAT‑ACGT‑‑ATAGGCCTGATAAGCGTAGCGCATCAGGCa                      >  1:321438/2‑100 (MQ=17)
                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                     >  1:370551/1‑100 (MQ=255)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                     >  1:2262768/1‑99 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTG‑GGCT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                    >  2:945169/9‑98 (MQ=21)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCA‑TGAG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                    >  1:146806/9‑98 (MQ=17)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                    >  1:2028258/1‑100 (MQ=255)
                                                                             ccGGATCCGGCGTTAACGCCTTATCCAGCCTACAAAATTGTGCAAAATCAATTGATTGCA‑CAGC‑TAGC‑‑GTAGGCCTGATGCGCTACGCTCATCAGGc                       <  2:1483766/97‑1 (MQ=18)
                                                                              cggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCA‑CCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                    >  2:1357700/8‑97 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                 >  1:2576818/7‑98 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                  >  1:1458984/7‑100 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                  >  2:1008095/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCA‑TGAG‑CCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                  >  2:1559778/7‑100 (MQ=14)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCA‑TGAG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                  >  1:131475/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                  >  2:1708710/7‑96 (MQ=25)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑GGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                   >  1:3034549/7‑99 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                  >  2:1239319/7‑96 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                 >  1:1066030/7‑98 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAAC‑GAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                  >  2:2269240/7‑96 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCATAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑ATAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACAttta                  >  1:2756742/7‑98 (MQ=18)
                                                                                 aTGGGGCGTTATCGCTTTATAGGGCGAAAAAAATCGTGGAAATTCACTACATTGCA‑GAAT‑TTTC‑‑GGAGGGCTAATAAGTGGAGCGCATCAGGCATTtttg                <  2:1974961/100‑1 (MQ=16)
                                                                                 aTGCGGGGTTAACTCCTTATCCGGCATACAAAATGGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTTATAAGCGGAGTGCATCAGGGATTtttg                <  2:1270139/100‑1 (MQ=18)
                                                                                   gcggcgTTAACGCCTTATCCAGCCTACAAAATAGTACAAAATCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGAGCATCAGGCATTtttgcc              >  2:704377/1‑100 (MQ=18)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct             <  1:2913962/100‑1 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct             >  1:2748791/1‑100 (MQ=255)
                                                                                     ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt            >  1:1614734/1‑100 (MQ=255)
                                                                                     ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCATGATAAGCGGAGCGCATCAGGCATTtttgcctt            >  1:6083/1‑100 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg          >  2:636789/1‑100 (MQ=255)
                                                                                         ttAACGCCTTATCCGGCCTACAAAATCGTGTAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg          <  2:2555109/98‑1 (MQ=255)
                                                                                         ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc        <  2:475307/100‑1 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt      >  1:2196911/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc     <  1:2420009/100‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GATT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg   <  1:773849/100‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg   <  1:1116001/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGT  >  NC_000913/3705887‑3706081

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: ‑