Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 28.4% 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.4% 553.5 / 196.8 246intergenic (‑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 (74/98);  new base T (62/3);  total (145/101)
Fisher's exact test for biased strand distribution p-value = 3.95e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

Only 100 of 252 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147
                                                                                                       |                                                                                                   
  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAg                                                                                                   <  2:2601910/100‑1 (MQ=255)
  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAg                                                                                                   <  1:3000541/100‑1 (MQ=255)
    ttccgatctggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TGAG‑CCACGTAGGCCTGATAAGCGTAgcg                                                                                                 <  1:2576287/90‑1 (MQ=17)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAgcg                                                                                                 >  1:1975033/1‑100 (MQ=255)
      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCa                                                                                               >  1:2205716/1‑100 (MQ=255)
       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCAt                                                                                              >  2:1291226/1‑100 (MQ=255)
        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATc                                                                                             >  2:251765/7‑100 (MQ=18)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                                                                                             >  2:1916605/1‑100 (MQ=255)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                                                                                             <  1:1275269/100‑1 (MQ=255)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                                                                                             >  2:2665227/1‑100 (MQ=255)
         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCa                                                                                            <  2:3181627/100‑1 (MQ=255)
         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCa                                                                                            <  1:1201295/100‑1 (MQ=255)
         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  1:2562065/6‑100 (MQ=18)
         tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  1:1468470/6‑100 (MQ=21)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAC‑AAAAACCTGCAAATTCAATATATTGCA‑GGAT‑ACGTGTAGGCCTGATAAGCGTAGCGCATCa                                                                                            >  2:1528971/5‑100 (MQ=17)
          gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           >  1:2095719/5‑100 (MQ=17)
           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  1:1107714/4‑100 (MQ=17)
           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACCTGCAAATTCAATATATTGCA‑GGAT‑ACGTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  2:1036689/4‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACTCCGCATCCGGCCT‑AC‑AAAAGCACGCAAATTCAATATATTGCA‑GAGA‑TTGCGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  2:720915/3‑100 (MQ=11)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                         >  2:761071/3‑100 (MQ=21)
             ggCCGGATAAGGCGTTCACTCCGCATCCGGCCT‑AC‑AAAAGCACGCAAATTCAATATATTGCA‑GAGA‑TTGCGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  2:2966949/2‑98 (MQ=14)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:2279719/2‑100 (MQ=17)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  2:3478606/1‑99 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑ATA‑TGACGTAGGCCTGATAAGCGTAGCGCATCAGGcg                                                                                        >  2:2232253/1‑97 (MQ=21)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑CCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                       >  2:1675406/1‑99 (MQ=21)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACCTGCAAATTCAATATATTGCA‑GGAT‑ACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                       >  2:3526736/1‑99 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGA‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                       >  1:2602598/1‑99 (MQ=17)
               ccggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                      >  2:3334806/9‑98 (MQ=18)
               ccggataaggcgTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaat                                                                                      >  1:2325546/9‑98 (MQ=25)
               ccggataaggcgTTCACGCCGCATCCGGCCT‑AC‑AAAATCATGAAAATTCAATACATTGCA‑AGAT‑TTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCacg                                                                                      >  1:318330/9‑98 (MQ=17)
               ccGGATAAGGCGTTTACGCCGCATCTGGCCT‑AC‑AAATTCTTCCAAATTCAATATATTGCA‑GAAA‑TCATGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                                      >  1:1059579/1‑98 (MQ=17)
                cGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                                                     >  2:237122/1‑100 (MQ=255)
                 ggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:1769186/7‑100 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAAGCGTGCAAATTCAATACATTGCATGGGC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                    >  2:1349059/7‑100 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GCGT‑TCACGTAGGCCTGATAAGCTTAGCGCATCAGGCAgttt                                                                                    >  1:2464083/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGGCCT‑AC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                                                   >  2:913768/7‑98 (MQ=21)
                   aTGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                  <  1:1847113/100‑1 (MQ=255)
                    tGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                                                 <  2:1774810/100‑1 (MQ=255)
                     ttccgatctaCGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                                <  1:1726945/91‑1 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               >  2:3567587/1‑100 (MQ=255)
                      aggcgTTCACGCCGCATCCGGCCT‑AC‑AAAACTGTGCAAATTCAATATATTGCA‑GGAA‑ACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTtttgcgt                                                                               >  1:1790180/2‑98 (MQ=18)
                         cgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            <  2:2004902/100‑1 (MQ=255)
                         cgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            >  2:1370108/1‑100 (MQ=255)
                         cgTTAACGCCTTACCCGGCTT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            <  2:3414128/100‑1 (MQ=255)
                            tAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                         >  1:344120/1‑100 (MQ=255)
                             aaCGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                                        <  2:560736/100‑1 (MQ=255)
                               cGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                                                      <  1:1441139/100‑1 (MQ=255)
                                 ccTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                                    >  2:1884261/1‑100 (MQ=255)
                                  cTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                                                   >  1:1691075/1‑100 (MQ=255)
                                   ttATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                                                  >  1:3615893/1‑100 (MQ=255)
                                   ttATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                                                  >  1:2770761/1‑100 (MQ=255)
                                    tATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                                                  >  1:776253/1‑99 (MQ=255)
                                     aTCCGGCAT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                                <  2:1597676/100‑1 (MQ=255)
                                      tCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                                               <  2:530584/100‑1 (MQ=255)
                                       ccGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              >  2:2052772/1‑100 (MQ=255)
                                           ccT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                          <  1:414956/100‑1 (MQ=255)
                                               aC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                       >  1:2746105/1‑100 (MQ=255)
                                                c‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                                                      >  2:1648127/1‑100 (MQ=255)
                                                   aaaTCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                                    >  2:2650189/1‑100 (MQ=255)
                                                       cGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                                                <  2:458910/100‑1 (MQ=255)
                                                        gTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                               >  2:4780/1‑100 (MQ=255)
                                                         tGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  1:1592076/100‑1 (MQ=255)
                                                         tGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  1:389982/100‑1 (MQ=255)
                                                         tGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  2:372538/100‑1 (MQ=255)
                                                          gCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  1:2940908/99‑1 (MQ=255)
                                                          gCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                             <  2:2804138/100‑1 (MQ=255)
                                                           cAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                            >  1:925011/1‑100 (MQ=255)
                                                           cAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                            >  1:2846846/1‑100 (MQ=255)
                                                           cAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                            >  2:1427673/1‑100 (MQ=255)
                                                            aaaTTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                           >  2:2464555/1‑100 (MQ=255)
                                                            aaaTTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                           >  1:1897162/1‑100 (MQ=255)
                                                            aaaTTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                           <  1:2665210/100‑1 (MQ=255)
                                                             aaTTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                                          <  1:816734/100‑1 (MQ=255)
                                                               ttCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                                        <  2:1691082/100‑1 (MQ=255)
                                                                tCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                                       <  1:3153635/100‑1 (MQ=255)
                                                                 cAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                      <  1:3315003/100‑1 (MQ=255)
                                                                     aCATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATAGGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtt                                  <  1:1604974/100‑1 (MQ=255)
                                                                          gCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             <  2:1084901/100‑1 (MQ=255)
                                                                          gCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             >  2:23971/1‑100 (MQ=255)
                                                                           cA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                            >  2:802413/1‑100 (MQ=255)
                                                                            a‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                           <  2:255645/100‑1 (MQ=255)
                                                                            a‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                           >  1:3362561/1‑100 (MQ=255)
                                                                              gAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          >  2:645859/1‑100 (MQ=255)
                                                                              gAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          <  2:3676558/100‑1 (MQ=255)
                                                                                aT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                        <  1:23951/100‑1 (MQ=255)
                                                                                      cGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGttt                   >  2:2728166/1‑100 (MQ=255)
                                                                                       gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                  <  2:3304201/100‑1 (MQ=255)
                                                                                       gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                  <  2:3219230/100‑1 (MQ=255)
                                                                                         aGGCCTGATAAGCGGAGCGCATCATGCATTTTTGCCTTTGCCATCATTCTGGTATTGCGTTTAATTATTCGAGAGAGACGTTTTCGAATTGATGTTTGcc                <  1:3563192/100‑1 (MQ=255)
                                                                                            ccTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTag             >  1:536079/1‑99 (MQ=255)
                                                                                            ccTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTaa             >  2:3387571/1‑100 (MQ=255)
                                                                                              tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  1:3576405/100‑1 (MQ=255)
                                                                                              tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  1:3261325/100‑1 (MQ=255)
                                                                                              tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  2:3461113/100‑1 (MQ=255)
                                                                                              tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  1:3564768/100‑1 (MQ=255)
                                                                                              tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  2:29046/100‑1 (MQ=255)
                                                                                              tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           >  1:2127773/1‑100 (MQ=255)
                                                                                               gATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg          >  2:1147043/1‑100 (MQ=255)
                                                                                                 tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt        <  2:1489810/100‑1 (MQ=255)
                                                                                                    gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa     <  2:2761539/100‑1 (MQ=255)
                                                                                                       |                                                                                                   
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147

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