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.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑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)
                 tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                          >  2:3562546/1‑100 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                        >  1:1866644/1‑100 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                        <  2:488648/100‑1 (MQ=255)
                      gatctggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                     <  1:100918/95‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                     <  2:1625288/100‑1 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                    <  2:2933578/100‑1 (MQ=255)
                        agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaa                                                                                   >  1:1928904/1‑100 (MQ=255)
                         gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                                  >  2:73991/1‑100 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                              >  1:620944/1‑100 (MQ=255)
                                 gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                          <  2:1315500/100‑1 (MQ=255)
                                  cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACa                                                                         >  1:3186095/1‑100 (MQ=255)
                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                       <  2:430284/100‑1 (MQ=255)
                                     cACCGGGTATACAGATTGTTGACAAAATGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCTGACCTAC‑AAGATCGTGCAAATTCCATAAAtt                                                                       >  1:665051/1‑99 (MQ=17)
                                      aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                     <  2:444333/100‑1 (MQ=255)
                                       ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                                    <  2:819002/100‑1 (MQ=255)
                                            tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑Att                                                           >  1:2698589/1‑100 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑Atttt                                                         <  2:637278/100‑1 (MQ=255)
                                               aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTc                                                        >  1:1134631/1‑100 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑Gt                                                    >  1:238974/1‑100 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑Gt                                                    <  1:3673704/100‑1 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTa                                                   <  2:1100419/100‑1 (MQ=255)
                                                   aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAg                                                  <  1:2145517/100‑1 (MQ=255)
                                                     tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAgcc                                                >  1:2342165/1‑98 (MQ=255)
                                                     cataaaCAAAGCGCACTTTGGTTATGCAGCATGCGGCGCTAACACCGTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTa                                                   >  1:2492197/6‑97 (MQ=17)
                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGcc                                               <  2:3522025/100‑1 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTg                                             <  2:1159066/100‑1 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCATTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTg                                             >  2:2688905/1‑100 (MQ=255)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGa                                            <  1:298523/100‑1 (MQ=255)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGa                                            >  2:765150/1‑100 (MQ=255)
                                                          aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGAt                                           >  1:2834820/1‑100 (MQ=255)
                                                          aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGAt                                           >  2:1519232/1‑100 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATa                                          <  1:423645/100‑1 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGAGGCGGCGTTACCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCgg                                     <  2:3364961/100‑1 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                  <  1:1257887/100‑1 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                  <  2:211256/100‑1 (MQ=255)
                                                                       tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑ATCTCCGTGTAGGCCTGATAAGCGTAGCGCAt                              >  2:2814988/8‑100 (MQ=18)
                                                                       tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAG‑‑‑G‑GACTAT‑‑GTAGGCCTGATAAGCATAGCGCAt                              >  2:2718370/8‑100 (MQ=18)
                                                                        ttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCAG‑‑‑A‑AATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATc                             >  2:1585578/7‑100 (MQ=17)
                                                                        ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAG‑‑‑A‑GTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATc                             >  1:1985931/7‑100 (MQ=17)
                                                                         tgtaggCCGGATAAGGCGTTCACGCCGCATCCAGCCTAC‑AAAATCATGAAAATTCAATACATTGCAA‑‑‑G‑ATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                            >  2:3651936/6‑100 (MQ=17)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑G‑ACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:2946257/5‑100 (MQ=18)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCAG‑‑‑AGATCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:1787920/5‑100 (MQ=17)
                                                                           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAG‑‑‑A‑AACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                          >  2:3749263/4‑100 (MQ=18)
                                                                           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                          >  2:3227214/1‑100 (MQ=255)
                                                                           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                          >  2:2254131/1‑100 (MQ=255)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAG‑‑‑A‑AACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:1266113/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAG‑‑‑G‑GACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                         >  1:3522269/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAG‑‑‑G‑AGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgaga                        >  2:2876384/3‑97 (MQ=21)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAG‑‑‑A‑AACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:3286393/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑G‑ACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:52625/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑ATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:258154/2‑100 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑G‑ACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  2:614621/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAG‑‑‑T‑TTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                       >  1:998545/1‑99 (MQ=21)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAG‑‑‑T‑TTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                       >  2:3051095/1‑99 (MQ=25)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCAG‑‑‑AGATTAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:2804639/1‑99 (MQ=17)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCAG‑‑‑A‑GATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  2:3699756/9‑98 (MQ=17)
                                                                               ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:1906813/1‑100 (MQ=255)
                                                                               ccGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCAG‑‑‑G‑AGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:3843047/1‑98 (MQ=18)
                                                                                cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑ATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                     >  2:1682055/8‑97 (MQ=18)
                                                                                cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                     >  2:2374711/1‑100 (MQ=255)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAG‑‑‑A‑AACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:588582/7‑100 (MQ=18)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCAT‑‑‑G‑GGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCActtt                    >  1:2643936/7‑100 (MQ=17)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAG‑‑‑A‑GTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:934387/7‑100 (MQ=17)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGCAATTCA‑ACTT‑‑‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                     >  2:929691/7‑97 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑G‑ACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:664762/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGATC‑‑‑G‑ATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:3099603/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGATC‑‑‑G‑ATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:3758047/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCAA‑‑‑G‑ATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                    >  2:39541/7‑96 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAG‑‑‑G‑AACGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                    >  1:1013021/7‑96 (MQ=21)
                                                                                 ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                    >  1:1435319/1‑100 (MQ=255)
                                                                                  gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCAG‑‑‑AGATCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattta                   >  2:1215837/6‑99 (MQ=17)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:866431/100‑1 (MQ=255)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:436968/100‑1 (MQ=255)
                                                                                     gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                <  1:763975/100‑1 (MQ=255)
                                                                                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               <  2:2436025/100‑1 (MQ=255)
                                                                                       ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCaaaaatgcctg              <  2:152752/100‑11 (MQ=255)
                                                                                       ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt              >  2:3453735/1‑100 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:2666780/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:1056723/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:1461450/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:2184453/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            >  2:1743360/1‑100 (MQ=255)
                                                                                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           <  2:3180371/100‑1 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         >  2:2007269/1‑100 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         <  1:1558964/100‑1 (MQ=255)
                                                                                             aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGTCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt        >  1:916708/1‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  2:1961149/100‑1 (MQ=255)
                                                                                               cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGTTTTTTCTATCa      >  1:1495488/1‑100 (MQ=17)
                                                                                                gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     >  1:202999/1‑100 (MQ=255)
                                                                                                gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  1:1743347/100‑1 (MQ=255)
                                                                                                   ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  <  1:3294168/100‑1 (MQ=255)
                                                                                                   |                                                                                                          
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑ATTTTC‑‑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: ‑