Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 28.2% 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,0480GT28.2% 601.4 / 228.5 313intergenic (‑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 (76/141);  new base T (64/18);  total (153/160)
Fisher's exact test for biased strand distribution p-value = 2.59e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 321 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147
                                                                                                         |                                                                                                   
tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCgg                                                                                                     <  1:3476098/100‑1 (MQ=255)
 gcgcTTTGTTTATGCCGGATTCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGa                                                                                                    >  2:2524683/1‑100 (MQ=255)
 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGa                                                                                                    <  1:542938/100‑1 (MQ=255)
 gcgcTTGGTGGATGCCGTAGGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGCAT‑‑TTTC‑‑GTAGGTATGATAAGCgg                                                                                                     <  2:3865933/99‑1 (MQ=17)
  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAg                                                                                                   <  2:685038/100‑1 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  >  1:3551034/1‑100 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  <  1:1844845/100‑1 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 >  1:2337018/1‑100 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 >  1:2234460/1‑100 (MQ=255)
     gttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAA‑‑ACAC‑‑GTAGGCCTGATAAGCGAAgcgc                                                                                                >  1:28312/10‑100 (MQ=18)
      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                                               <  1:635910/100‑1 (MQ=255)
      atttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAA‑‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCa                                                                                               >  1:253299/9‑100 (MQ=18)
       tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGAT‑‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                                              >  1:3894808/8‑100 (MQ=18)
         tttATGCCGGATGCGTCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                                            <  1:4070090/100‑1 (MQ=255)
         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                                            >  1:3654130/1‑100 (MQ=255)
         tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                            >  2:586815/6‑100 (MQ=18)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCT‑‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  2:452644/5‑100 (MQ=21)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAT‑‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  2:2685580/5‑100 (MQ=17)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAGTCGTATAAATTCAATGTATTGCAGTTT‑‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAg                                                                                           >  1:2081561/5‑100 (MQ=21)
          gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           >  2:368802/5‑100 (MQ=18)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                          <  1:92434/100‑1 (MQ=255)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                          <  2:2747110/100‑1 (MQ=255)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTAC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           >  2:2564800/1‑99 (MQ=255)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                                                                                         >  2:2992827/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGA‑‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                                                                                         >  1:3151707/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGAT‑‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:3217352/3‑100 (MQ=18)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                         >  1:1037291/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:1277254/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTT‑‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                                        >  1:1462990/2‑100 (MQ=21)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAAT‑‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:3000158/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGAT‑‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:4063685/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAG‑‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                                       >  2:2929945/2‑100 (MQ=21)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  1:2663455/1‑99 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAG‑‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  2:1063557/1‑99 (MQ=17)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAG‑‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  1:1226572/1‑99 (MQ=17)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAA‑‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCag                                                                                       >  2:383007/1‑99 (MQ=18)
              gCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                       >  1:2356037/1‑100 (MQ=18)
                tggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTT‑‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatc                                                                                     <  1:3426685/93‑4 (MQ=25)
                cggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAG‑‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                                                                                     >  2:3637288/8‑100 (MQ=17)
                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAA‑‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:2875162/7‑100 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGT‑‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:2021854/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:875970/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTT‑‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                                                    >  2:3104334/7‑96 (MQ=25)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAG‑‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                                                   >  2:2104867/7‑98 (MQ=21)
                 ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                                                                                    >  1:878592/1‑100 (MQ=255)
                  gataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAA‑‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAgtttt                                                                                   >  1:3978452/6‑100 (MQ=18)
                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                  <  1:2634783/100‑1 (MQ=255)
                     gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                                <  1:694945/100‑1 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCACTACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATGcttgcct                                                                               >  2:4120221/1‑100 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               >  1:39669/1‑100 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               >  1:959821/1‑100 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            <  1:3804816/100‑1 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            >  2:1152345/1‑100 (MQ=255)
                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                          <  1:1509777/100‑1 (MQ=255)
                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                          <  1:794578/100‑1 (MQ=255)
                                gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATctg                                                                     >  2:2496787/1‑98 (MQ=255)
                                    tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                                                 >  1:3398643/1‑100 (MQ=255)
                                     aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                                <  2:800362/100‑1 (MQ=255)
                                     aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                                <  1:2306532/100‑1 (MQ=255)
                                         ggCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                            <  2:2202202/100‑1 (MQ=255)
                                           ccTACAAAATCGTGCTAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                                                           <  2:2418874/99‑1 (MQ=255)
                                           ccTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                          <  2:1336299/100‑1 (MQ=255)
                                           ccTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                          <  2:4040525/100‑1 (MQ=255)
                                               cAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                                                      <  1:2601241/100‑1 (MQ=255)
                                                      cGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTTATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                                                <  1:1184777/100‑1 (MQ=255)
                                                      cGTGAAAATGCAATAAATTGCATCAT‑‑CTTC‑‑TTAGCCTTGTTAAGCGGAGCTCATCAGGCATTTTTGCCTTTCCCATCAGTCCTGTATGGCTTTTAatt                                                  <  2:1721596/98‑1 (MQ=255)
                                                       gTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                               <  2:2337005/100‑1 (MQ=255)
                                                        tGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              >  2:1454734/1‑100 (MQ=255)
                                                        tGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              >  1:1909469/1‑100 (MQ=255)
                                                         gCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                             <  2:3551011/100‑1 (MQ=255)
                                                          cAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                            <  2:1530810/100‑1 (MQ=255)
                                                            aaTTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAtn                                          <  1:3333189/100‑2 (MQ=255)
                                                             aTTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                                         <  1:1107641/100‑1 (MQ=255)
                                                              ttCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                           <  2:1251386/97‑1 (MQ=255)
                                                                  aTACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg                                    <  2:1852365/100‑1 (MQ=255)
                                                                     cATTGCAGACT‑‑TTTC‑‑GTCGGCCTGATAAGCGGCGCGCATCAGTCATTTTTGCCTTTGCCATCAGTCCTGTACGGCCTTTAATTATTCGATAGAGATGttt                                 >  2:1966019/1‑100 (MQ=255)
                                                                       ttGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                               >  1:3267176/1‑100 (MQ=255)
                                                                        tGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCg                              <  2:3223579/100‑1 (MQ=255)
                                                                           aGAAT‑‑TTTC‑‑GTAGGCCTGTTAAGCGGAGCGCATCAGGCATTTTTGGCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                           <  2:1796508/100‑1 (MQ=255)
                                                                             gAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          >  2:1766039/1‑100 (MQ=255)
                                                                             gAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          <  1:1939843/100‑1 (MQ=255)
                                                                             gAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          >  1:2952919/1‑100 (MQ=255)
                                                                                t‑‑tttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt                       <  1:1424498/100‑1 (MQ=255)
                                                                                   tttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCctacgaaaattctgcaatgtattgaagttgcaagattttgtagtccggataaggcgttaacgcc                         <  2:1032666‑M2/97‑65 (MQ=255)
                                                                                   tttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATg                      <  2:3398623/100‑1 (MQ=255)
                                                                                     tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGtt                    <  2:2619780/100‑1 (MQ=255)
                                                                                           aGGCCTGATAAGCGGAGCGCATCAGGCctacgaaaattctgcaatgtattgaatttgcacgattttgtaggccggataaggcgttaacgccgcatccggc                <  2:1208489‑M2/100‑74 (MQ=255)
                                                                                            ggCCTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgccgcatccggca               <  1:3105598‑M2/100‑75 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCctacgaaaattctgcaatgtattgaatttgcacgattttgtaggccggataaggcgttaacgccgcatccggca               >  2:1433885‑M2/1‑26 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCTATCAGTCTTGTATGGCTTTTAATTATTCGATATAGACGTTTTCGAAGTGATGTTTGCCt               >  2:519829/1‑100 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               >  2:2270031/1‑100 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               <  1:4091775/100‑1 (MQ=255)
                                                                                             gCCTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgccgcatccggcaa              <  1:1383917‑M2/100‑76 (MQ=255)
                                                                                              ccTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgccgcatccggcaag             <  2:919671‑M2/100‑77 (MQ=255)
                                                                                              ccTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgatttggtaggccggataaggcgttcacgccgcatccggcaag             <  2:3292919‑M2/100‑77 (MQ=255)
                                                                                                tGATAAGCGTAGCGCATCAGGCctacatggtcctgcaatatattgaattggcaagatttttgtaggccggataaggcgttcacgccgcatccggcatgaa           <  1:1615250‑M2/100‑79 (MQ=255)
                                                                                                tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  2:1626204/100‑1 (MQ=255)
                                                                                                   tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt        <  1:2785223/100‑1 (MQ=255)
                                                                                                     aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa      >  1:3369503/1‑100 (MQ=255)
                                                                                                        ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc   <  1:1795122/100‑1 (MQ=255)
                                                                                                         |                                                                                                   
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147

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