Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 24.8% intergenic (‑258/+50) 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,706,0480GT24.8% 611.8 / 182.0 284intergenic (‑258/+50)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (101/107);  new base T (67/2);  total (175/109)
Fisher's exact test for biased strand distribution p-value = 2.61e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146
                                                                                                            |                                                                                                  
tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCgg                                                                                                    <  1:3151016/100‑1 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 <  1:2370531/100‑1 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 >  2:2706405/1‑100 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                >  1:1099675/1‑100 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                >  1:1172899/1‑100 (MQ=255)
      atttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCa                                                                                              >  1:2628907/9‑100 (MQ=18)
       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                                                             >  2:1263790/1‑100 (MQ=255)
        gTTTATGCCGGATGCGGCGTTAACGCCTTCTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                            >  2:3658973/1‑100 (MQ=255)
         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                           >  1:691795/6‑100 (MQ=14)
          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATAGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                          <  2:2414482/100‑1 (MQ=255)
          gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAATCTTC‑‑‑TG‑‑GTAGGCCTGATAAGCATCGCGCATCAg                                                                                          >  1:19611/5‑100 (MQ=18)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCA‑‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                          >  1:167847/5‑100 (MQ=18)
           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACCAAACTTTGCAAATTCAATATATTGC‑‑AATCT‑‑‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                         >  2:2095046/4‑100 (MQ=18)
            nggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  2:3164560/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAAC‑‑‑AC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc                                                                                        >  1:225326/3‑100 (MQ=18)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  2:3561352/1‑100 (MQ=255)
             tGCCGGATGCGGCGTTAACGCCTTATACGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAAGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  2:290824/1‑100 (MQ=255)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATGTTGCAGAAACC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:2504201/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATCC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:2267474/2‑100 (MQ=21)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  2:1052906/1‑100 (MQ=255)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  1:1622948/1‑100 (MQ=255)
              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                                                     >  1:3468936/1‑100 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCA‑‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                      >  2:2806682/1‑99 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCC‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                      >  1:2550234/1‑99 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTA‑‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                      >  2:1964773/1‑99 (MQ=17)
               caggataaggcgTTCACGCCGCATCCTGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATG‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                                                     >  1:2072162/9‑100 (MQ=14)
                cggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                                                                                   >  1:2950735/8‑99 (MQ=18)
                cggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCA‑‑‑T‑‑GTAAACCTGATAAGCGTAGCGCATCAGGCAgtt                                                                                    >  1:2114005/8‑100 (MQ=18)
                cggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                                                                                    >  2:1017215/8‑100 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  2:3331247/7‑100 (MQ=17)
                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAC‑‑‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCAatt                                                                                    >  1:1546523/7‑100 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑AATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                                                                                    >  1:440205/7‑97 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCA‑‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  2:393529/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCA‑‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:1605903/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  1:2016312/7‑96 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTAT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                                                   >  2:59395/7‑96 (MQ=25)
                 ggataaggcgTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                                                                                   >  2:534952/7‑100 (MQ=18)
                   aTGCGGCGTTAACTCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:2408202/100‑1 (MQ=255)
                    ttcggcgTTAACGACTTATCCGGCCTACAAAATCGTGCAAATACAATACATTGCAGAATTT‑‑‑TC‑‑GTAAGCCTGATAAGCGGAGCGCATCAGTCATTtttgc                                                                                >  1:2650609/3‑100 (MQ=18)
                    taaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCC‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg                                                                                 >  1:1830812/4‑95 (MQ=18)
                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                              <  1:2861154/100‑1 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTACAGAATTA‑‑‑TC‑‑ATAAACATAATAAGCGGAGCACATCAGGCATTtttgcct                                                                              >  1:1370461/1‑100 (MQ=18)
                       ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                                             <  2:1382847/100‑1 (MQ=255)
                       ggcaaacatcacttcgaaaacgtctctatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  1:3102815‑M1/74‑100 (MQ=255)
                       ggcaaacatcacttcgaaaacgtctctatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  1:1333199‑M1/74‑100 (MQ=255)
                          gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAGTTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                                          <  1:1265379/100‑1 (MQ=255)
                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                         <  1:1200613/100‑1 (MQ=255)
                            acatcacttcgaaaacgtctctatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 >  1:1786730‑M1/69‑100 (MQ=255)
                             aaCGCCTTATCCGACCTACAAAATCGTGCAAATTCAATTCATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                                       >  2:2408127/1‑100 (MQ=255)
                               cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATAAGGCATTTTTGCCTTTGCCATCa                                                                     >  1:1015696/1‑100 (MQ=255)
                                gCCTTATCCGGTCTACAAAATCGTGCAAATTCAATACATTACAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGTGCGTATCAGGCATTTTTGCCTTTGCCATCAg                                                                    <  2:915478/100‑1 (MQ=21)
                                 ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                                   <  2:1455428/100‑1 (MQ=255)
                                     aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                               >  1:2584491/1‑100 (MQ=255)
                                     aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                               <  2:3762462/100‑1 (MQ=255)
                                     aTCCGGCCTAAAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                               <  1:802008/100‑1 (MQ=255)
                                      tCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                                              >  1:622140/1‑100 (MQ=255)
                                        cGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                                            >  1:1044028/1‑100 (MQ=255)
                                         ggCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                           <  1:564194/100‑1 (MQ=255)
                                           ccTACAAAATCGTGCAAATTCTATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                         >  2:3190842/1‑100 (MQ=255)
                                           ccTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                         >  1:2723415/1‑100 (MQ=255)
                                              aCAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       <  2:1205512/99‑1 (MQ=255)
                                                  aaaTCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                                   <  1:3313972/100‑1 (MQ=255)
                                                   aaTCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                  <  1:887696/100‑1 (MQ=255)
                                                   aaTCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                  <  2:2324906/100‑1 (MQ=255)
                                                       gTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                              >  2:2149316/1‑100 (MQ=255)
                                                        tGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:3714728/100‑1 (MQ=255)
                                                        tGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:3215207/100‑1 (MQ=255)
                                                        tGCAAATTCAATACACTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:3868206/100‑1 (MQ=255)
                                                         gCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            <  1:3402754/100‑1 (MQ=255)
                                                          cAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatagatcgg                                           >  2:2603597/1‑93 (MQ=255)
                                                          cAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                           <  2:1989175/100‑1 (MQ=255)
                                                           aaaTTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATAAGTCTTGTATGGCTTTTAATTATTCGAt                                          >  1:484995/1‑100 (MQ=255)
                                                                 aaTACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc                                    >  1:775810/1‑100 (MQ=255)
                                                                       ttGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                              >  2:1691914/1‑100 (MQ=255)
                                                                         gCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            >  2:1556844/1‑100 (MQ=255)
                                                                         gCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            <  2:775801/100‑1 (MQ=255)
                                                                         gCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            >  1:1472838/1‑100 (MQ=255)
                                                                         gCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            >  1:2394447/1‑100 (MQ=255)
                                                                         gCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            <  2:3025539/100‑1 (MQ=255)
                                                                          cAGAATTT‑‑‑TC‑‑GTAGGCATGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                           >  1:1457587/1‑100 (MQ=255)
                                                                           aGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  1:668154/100‑1 (MQ=255)
                                                                           aGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  1:3616515/100‑1 (MQ=255)
                                                                             aaTTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                        <  2:3462058/100‑1 (MQ=255)
                                                                                ttt‑‑‑tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt                      >  1:591142/1‑100 (MQ=255)
                                                                                ttt‑‑‑tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt                      >  1:279712/1‑100 (MQ=255)
                                                                                  tt‑‑‑tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATg                     <  2:3312742/100‑1 (MQ=255)
                                                                                    t‑‑‑tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGt                    <  2:2194126/100‑1 (MQ=255)
                                                                                            gTAGGCTTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 <  1:2662858/100‑1 (MQ=255)
                                                                                            gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 <  1:2459902/100‑1 (MQ=255)
                                                                                              aGGCTTGATAAGCGGAGCGGATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  2:3110103/100‑1 (MQ=255)
                                                                                               ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGAGAtcggaagag              >  2:799123/1‑92 (MQ=38)
                                                                                                  CTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           >  1:287892‑M1/1‑100 (MQ=255)
                                                                                                  CTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           <  1:2294402‑M1/100‑1 (MQ=255)
                                                                                                   tGATAAGCGGATCGCATCAGGCATTTTTTTCTTTTTCTTCTCTCTTATATGACTTTTAATTATTCGATACAGACGTTTTCTCAGTGATGTTTGCCTAAt           >  2:1805848/1‑99 (MQ=255)
                                                                                                    GATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         <  1:3045446‑M1/100‑1 (MQ=255)
                                                                                                      TAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt       >  1:39681‑M1/1‑100 (MQ=255)
                                                                                                       aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc      <  1:1691929/100‑1 (MQ=255)
                                                                                                        aGCGGAGCGCATCAGGCATTTTTGTCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     <  1:3584991/100‑1 (MQ=255)
                                                                                                        aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATTTTTGCCTAATGGATCa     >  2:2489706/1‑100 (MQ=255)
                                                                                                        aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     >  2:3583350/1‑100 (MQ=255)
                                                                                                            |                                                                                                  
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146

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