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.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑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)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                              <  2:499872/100‑1 (MQ=255)
                                        ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                            >  1:1141105/1‑100 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCTGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTTCAAATTCAATACATTGCAGAAtttt                                                       >  2:117613/1‑100 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATCACGGATGCAGCGTTAACTCCTTATCCGGCCTACAAAATCGTGCAAATTCATTACATTGCAGAAtttt                                                       <  2:85273/100‑1 (MQ=17)
                                            taCAGATTGTTGAAAAAGTGCGCTTTGTTTATGCCGGACGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtttt                                                       >  2:1871446/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                   >  1:2857584/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                   >  1:1702126/1‑100 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCATATTCAATACATTGCAGAATTTTC‑‑GTa                                                 <  1:1843893/100‑1 (MQ=255)
                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCTGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAgg                                               <  1:1201207/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGTTTTATCCTGTTTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGcc                                             <  1:983872/100‑1 (MQ=18)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTg                                           >  2:759821/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGa                                          >  2:1712624/1‑100 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGAt                                         <  2:2811126/100‑1 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCAGGATGCGGCGTTAAAGCCTTATCAGGCATACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGAt                                         >  2:2140109/1‑100 (MQ=18)
                                                           aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAg                                      <  1:2655613/100‑1 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGc                                     >  1:2827266/1‑100 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCgg                                   <  1:330353/100‑1 (MQ=255)
                                                                cgcTTTGTTGATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAg                                 <  2:2570162/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                >  2:373331/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAGAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                               <  2:1944875/100‑1 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                <  1:405163/99‑1 (MQ=255)
                                                                   tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcgc                              >  1:746255/1‑100 (MQ=255)
                                                                     tttgtaggCCGGATAAGTCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTACAGAGATTGC‑‑GTAGGCCTGATAAGCGTAg                                 >  2:1462092/8‑95 (MQ=11)
                                                                     tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                            >  1:2836239/8‑100 (MQ=18)
                                                                     tGTTTCTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAATTTC‑‑GTACGCCTGATAAGCGGAgcgc                              >  2:1938451/1‑98 (MQ=255)
                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAATGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                            >  2:1671444/1‑100 (MQ=255)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATc                           >  2:2399942/7‑100 (MQ=18)
                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                           <  2:2326732/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGAATCa                          >  1:577411/1‑100 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGACTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTATC‑‑GTAGGCCTGATAAACGGAGCGCATCa                          >  2:2268461/1‑100 (MQ=255)
                                                                       tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                          >  1:1798739/6‑100 (MQ=18)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                         >  2:3020255/5‑100 (MQ=18)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAg                         <  1:2399970/96‑1 (MQ=18)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                        >  2:3026071/4‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                       >  1:412280/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                       >  1:396157/3‑100 (MQ=17)
                                                                           ggccagataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGACAAGCGTAGCGCATCAGGCa                      >  2:2877078/11‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                      >  1:85464/2‑100 (MQ=14)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                      >  1:509184/2‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTAT‑‑GTAGACCTGATAAGCATAGCGCATCAGGCa                      >  1:1212855/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGT‑‑ATAGGCCTGATAAGCGTAGCGCATCAGGCa                      >  1:321438/2‑100 (MQ=17)
                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                     >  1:370551/1‑100 (MQ=255)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                     >  1:2262768/1‑99 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                    >  2:945169/9‑98 (MQ=21)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                    >  1:146806/9‑98 (MQ=17)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                    >  1:2028258/1‑100 (MQ=255)
                                                                             ccGGATCCGGCGTTAACGCCTTATCCAGCCTACAAAATTGTGCAAAATCAATTGATTGCACAGCTAGC‑‑GTAGGCCTGATGCGCTACGCTCATCAGGc                       <  2:1483766/97‑1 (MQ=18)
                                                                              cggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                    >  2:1357700/8‑97 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                 >  1:2576818/7‑98 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                  >  1:1458984/7‑100 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                  >  2:1008095/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                  >  2:1559778/7‑100 (MQ=14)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                  >  1:131475/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                  >  2:1708710/7‑96 (MQ=25)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                   >  1:3034549/7‑99 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                  >  2:1239319/7‑96 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                 >  1:1066030/7‑98 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                  >  2:2269240/7‑96 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCATACAAATCATGCAAATCCAAAATATTGCAGGAGATAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACAttta                  >  1:2756742/7‑98 (MQ=18)
                                                                                 aTGGGGCGTTATCGCTTTATAGGGCGAAAAAAATCGTGGAAATTCACTACATTGCAGAATTTTC‑‑GGAGGGCTAATAAGTGGAGCGCATCAGGCATTtttg                <  2:1974961/100‑1 (MQ=16)
                                                                                 aTGCGGGGTTAACTCCTTATCCGGCATACAAAATGGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTTATAAGCGGAGTGCATCAGGGATTtttg                <  2:1270139/100‑1 (MQ=18)
                                                                                   gcggcgTTAACGCCTTATCCAGCCTACAAAATAGTACAAAATCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGAGCATCAGGCATTtttgcc              >  2:704377/1‑100 (MQ=18)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct             <  1:2913962/100‑1 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct             >  1:2748791/1‑100 (MQ=255)
                                                                                     ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt            >  1:1614734/1‑100 (MQ=255)
                                                                                     ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCATGATAAGCGGAGCGCATCAGGCATTtttgcctt            >  1:6083/1‑100 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg          >  2:636789/1‑100 (MQ=255)
                                                                                         ttAACGCCTTATCCGGCCTACAAAATCGTGTAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg          <  2:2555109/98‑1 (MQ=255)
                                                                                         ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc        <  2:475307/100‑1 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt      >  1:2196911/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc     <  1:2420009/100‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGATTTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg   <  1:773849/100‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg   <  1:1116001/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑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: