Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 31.7% intergenic (‑194/+114) 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,9840TG31.7% 438.0 / 254.3 239intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (77/86);  new base G (69/7);  total (146/93)
Fisher's exact test for biased strand distribution p-value = 1.53e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

Only 100 of 244 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083
                                                                                                   |                                                                                                          
ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                            <  2:702400/100‑1 (MQ=255)
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                          <  1:1195709/100‑1 (MQ=255)
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                          >  2:3191822/1‑100 (MQ=255)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                    >  1:672778/1‑100 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                   <  1:3470915/100‑1 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                   <  1:282890/100‑1 (MQ=255)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                  <  1:3693956/100‑1 (MQ=255)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                  <  2:2060781/100‑3 (MQ=37)
           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                                 <  2:3008974/100‑1 (MQ=37)
                 tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                                          >  2:3562546/1‑100 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                        >  1:1866644/1‑100 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                        <  2:488648/100‑1 (MQ=255)
                      gatctggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                     <  1:100918/95‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                     <  2:1625288/100‑1 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                    <  2:2933578/100‑1 (MQ=255)
                        agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaa                                                                                   >  1:1928904/1‑100 (MQ=255)
                         gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                                  >  2:73991/1‑100 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                              >  1:620944/1‑100 (MQ=255)
                                 gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                          <  2:1315500/100‑1 (MQ=255)
                                  cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACa                                                                         >  1:3186095/1‑100 (MQ=255)
                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                       <  2:430284/100‑1 (MQ=255)
                                     cACCGGGTATACAGATTGTTGACAAAATGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCTGACCTACAAGATCGTGCAAATTCCATAAAtt                                                                       >  1:665051/1‑99 (MQ=17)
                                      aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                     <  2:444333/100‑1 (MQ=255)
                                       ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                                    <  2:819002/100‑1 (MQ=255)
                                            tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AAtt                                                           >  1:2698589/1‑100 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AAtttt                                                         <  2:637278/100‑1 (MQ=255)
                                               aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTc                                                        >  1:1134631/1‑100 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑Gt                                                    >  1:238974/1‑100 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑Gt                                                    <  1:3673704/100‑1 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTa                                                   <  2:1100419/100‑1 (MQ=255)
                                                   aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAg                                                  <  1:2145517/100‑1 (MQ=255)
                                                     tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAgcc                                                >  1:2342165/1‑98 (MQ=255)
                                                     cataaaCAAAGCGCACTTTGGTTATGCAGCATGCGGCGCTAACACCGTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTa                                                   >  1:2492197/6‑97 (MQ=17)
                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGcc                                               <  2:3522025/100‑1 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTg                                             <  2:1159066/100‑1 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCATTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTg                                             >  2:2688905/1‑100 (MQ=255)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGa                                            <  1:298523/100‑1 (MQ=255)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGa                                            >  2:765150/1‑100 (MQ=255)
                                                          aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGAt                                           >  1:2834820/1‑100 (MQ=255)
                                                          aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGAt                                           >  2:1519232/1‑100 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATa                                          <  1:423645/100‑1 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGAGGCGGCGTTACCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCgg                                     <  2:3364961/100‑1 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                  <  1:1257887/100‑1 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                  <  2:211256/100‑1 (MQ=255)
                                                                       tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCAt                              >  2:2814988/8‑100 (MQ=18)
                                                                       tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAG‑‑‑GGACTAT‑‑GTAGGCCTGATAAGCATAGCGCAt                              >  2:2718370/8‑100 (MQ=18)
                                                                        ttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAG‑‑‑AAATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATc                             >  2:1585578/7‑100 (MQ=17)
                                                                        ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAG‑‑‑AGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATc                             >  1:1985931/7‑100 (MQ=17)
                                                                         tgtaggCCGGATAAGGCGTTCACGCCGCATCCAGCCTACAAAATCATGAAAATTCAATACATTGCAA‑‑‑GATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                            >  2:3651936/6‑100 (MQ=17)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑GACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:2946257/5‑100 (MQ=18)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAG‑‑‑AGATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:1787920/5‑100 (MQ=17)
                                                                           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAG‑‑‑AAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                          >  2:3749263/4‑100 (MQ=18)
                                                                           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                          >  2:3227214/1‑100 (MQ=255)
                                                                           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                          >  2:2254131/1‑100 (MQ=255)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAG‑‑‑AAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:1266113/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAG‑‑‑GGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                         >  1:3522269/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAG‑‑‑GAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgaga                        >  2:2876384/3‑97 (MQ=21)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAG‑‑‑AAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:3286393/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑GACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:52625/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:258154/2‑100 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑GACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  2:614621/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAG‑‑‑TTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                       >  1:998545/1‑99 (MQ=21)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAG‑‑‑TTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                       >  2:3051095/1‑99 (MQ=25)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAG‑‑‑AGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:2804639/1‑99 (MQ=17)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAG‑‑‑AGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  2:3699756/9‑98 (MQ=17)
                                                                               ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:1906813/1‑100 (MQ=255)
                                                                               ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAG‑‑‑GAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:3843047/1‑98 (MQ=18)
                                                                                cggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                     >  2:1682055/8‑97 (MQ=18)
                                                                                cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                     >  2:2374711/1‑100 (MQ=255)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAG‑‑‑AAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:588582/7‑100 (MQ=18)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCAT‑‑‑GGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCActtt                    >  1:2643936/7‑100 (MQ=17)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAG‑‑‑AGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:934387/7‑100 (MQ=17)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                     >  2:929691/7‑97 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑GACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:664762/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATC‑‑‑GATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:3099603/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATC‑‑‑GATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:3758047/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAA‑‑‑GATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                    >  2:39541/7‑96 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAG‑‑‑GAACGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                    >  1:1013021/7‑96 (MQ=21)
                                                                                 ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                    >  1:1435319/1‑100 (MQ=255)
                                                                                  gataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAG‑‑‑AGATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattta                   >  2:1215837/6‑99 (MQ=17)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:866431/100‑1 (MQ=255)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:436968/100‑1 (MQ=255)
                                                                                     gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                <  1:763975/100‑1 (MQ=255)
                                                                                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               <  2:2436025/100‑1 (MQ=255)
                                                                                       ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCaaaaatgcctg              <  2:152752/100‑11 (MQ=255)
                                                                                       ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt              >  2:3453735/1‑100 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:2666780/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:1056723/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:1461450/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:2184453/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            >  2:1743360/1‑100 (MQ=255)
                                                                                          gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           <  2:3180371/100‑1 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         >  2:2007269/1‑100 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         <  1:1558964/100‑1 (MQ=255)
                                                                                             aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGTCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt        >  1:916708/1‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  2:1961149/100‑1 (MQ=255)
                                                                                               cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGTTTTTTCTATCa      >  1:1495488/1‑100 (MQ=17)
                                                                                                gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     >  1:202999/1‑100 (MQ=255)
                                                                                                gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  1:1743347/100‑1 (MQ=255)
                                                                                                   ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  <  1:3294168/100‑1 (MQ=255)
                                                                                                   |                                                                                                          
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083

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