Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 32.5% 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,0480GT32.5% 495.4 / 255.2 272intergenic (‑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 (80/98);  new base T (81/3);  total (170/102)
Fisher's exact test for biased strand distribution p-value = 3.08e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

Only 100 of 282 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147
                                                                                                          |                                                                                                   
   gcTTTGTTTATGCCGGATGCGTCGTTAATGCCTTTTTTTGTTT‑AC‑AAAATCGCGCCAACTCAAAACATTG‑CA‑GAAG‑TTTC‑‑ATAGGCCTGATAAGCGCAg                                                                                                   <  2:1717930/99‑1 (MQ=12)
       tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAC‑AAAAACCTGCAAATTCAATATATTG‑CA‑GGAT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCa                                                                                               >  1:3159563/8‑100 (MQ=12)
        ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                            >  2:3139033/7‑100 (MQ=18)
        ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCAAATTCAACATATTG‑C‑‑‑AAT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATc                                                                                             >  1:1249574/7‑100 (MQ=18)
        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATc                                                                                             >  1:3864568/7‑100 (MQ=21)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                             >  1:875368/1‑100 (MQ=255)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                             >  2:2122604/1‑100 (MQ=255)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                             <  1:1811/100‑1 (MQ=255)
          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           <  1:2980829/100‑1 (MQ=255)
          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           <  2:1140647/100‑1 (MQ=255)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           <  1:2324286/96‑1 (MQ=18)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  2:2324274/5‑100 (MQ=18)
           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  1:890024/4‑100 (MQ=18)
           tATGCCGGCTGCGGCGTTAACGCCTTCTCCGGCCT‑AC‑TAAATCGTGCAAATTCAATACATTG‑CA‑TAAT‑TTTC‑‑GTATTGCTGATAAGCGGAGCGCATCAgg                                                                                          <  2:2712831/100‑1 (MQ=17)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                          <  2:150553/100‑1 (MQ=255)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                          <  2:3018724/100‑1 (MQ=255)
            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:3274975/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAAGCGTGCAAATTCAATACATTG‑CATGGGC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:2734088/3‑100 (MQ=17)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:1742172/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTG‑CA‑GTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                         >  1:2884127/3‑100 (MQ=21)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:2163579/3‑100 (MQ=21)
            aggCCGGATAAGGCGTTAACGCCGCATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                                         >  1:1418089/3‑100 (MQ=18)
             ggtcggataaggcgTTCACGCCGCATCCGACCC‑AC‑AAAGTCTTGCAAAATCAATAAATTG‑CAGGAAC‑GAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:1896532/11‑100 (MQ=21)
             ggCCGGATAAGGGGCTTACGCCGCATCCGGCAT‑AC‑AAGAACGTGCAAATTCAACAAATTGTGATGCAT‑CCT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  2:746507/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:1631156/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAGATCGTGCAAATTCAATATATTG‑CA‑TGAG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:135189/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:2571781/2‑100 (MQ=21)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACATGCAAATTCAATAGATTG‑CA‑GAGA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:1747833/2‑100 (MQ=17)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GACT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                       >  2:1059684/1‑100 (MQ=255)
              gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCAAATTCAACATATTG‑CA‑ATTC‑TCTT‑‑GTAGGTCTGATAAGCGTAGCGCATCAGGctg                                                                                       >  1:563498/1‑98 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  2:3093156/1‑99 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTG‑CA‑GTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                                       >  1:1109645/1‑99 (MQ=25)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTG‑CA‑GTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                                       >  1:1590760/1‑99 (MQ=21)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                       >  2:592406/1‑99 (MQ=21)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAGCATGAAAATTCAATATATTG‑CA‑GGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  1:3892116/1‑99 (MQ=17)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAACTGTGCAAATTCAATATATTG‑CA‑GGAA‑ACAC‑‑GTAGGCCTGAT‑AGCGAAGCGCATCAGGCagt                                                                                      >  2:1406779/1‑98 (MQ=18)
              gCCGGATAAGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                       >  2:202301/1‑100 (MQ=255)
               ccggataaggcgTTCACGCCGCATCCGGCCT‑AC‑AAAAGCATGAAAATTCAATATATTG‑CA‑GGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                      >  1:3074817/9‑98 (MQ=18)
               ccggataaggcgTTCACGCCGCATCCGGCCT‑AC‑‑AAATCATGCAAATCCAAAATATTG‑CA‑GGAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                                                     >  1:3738540/9‑100 (MQ=21)
               ccGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAACGCATCAGGc                                                                                         >  2:1849155/1‑97 (MQ=17)
                cggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt                                                                                     >  1:3495562/8‑97 (MQ=17)
                cggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAACAGTGCAAATTCAATATATTG‑CA‑ATCT‑TCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaagc                                                                                     >  1:1637391/8‑97 (MQ=21)
                cggataaggcgTTCACGCCGCATCCGGCCT‑AC‑AAAATTGTGTAAATTCAATAACTTG‑AT‑CGAT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                                                                                     >  2:670402/8‑100 (MQ=21)
                 ggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:2632075/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCTGCCT‑AC‑AAAAGCACGCAAATTCAATATATTG‑CA‑GAGA‑TGGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttt                                                                                    >  1:2198275/7‑100 (MQ=14)
                 ggataaggcgTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:3946351/7‑100 (MQ=18)
                   aTGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                  <  2:3235920/100‑1 (MQ=255)
                     gcggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                                <  1:2015218/100‑1 (MQ=255)
                     gcggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                                <  1:985615/100‑1 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               >  1:36715/1‑100 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               >  2:35405/1‑100 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               >  2:1493828/1‑100 (MQ=255)
                        gcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt                                                                             >  2:2926462/1‑100 (MQ=255)
                         cgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            <  2:1612480/100‑1 (MQ=255)
                         cgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            <  2:3595180/100‑1 (MQ=255)
                           ttAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                          <  1:3227671/100‑1 (MQ=255)
                             aaCGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                                        >  1:2086203/1‑100 (MQ=255)
                              aCGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       <  2:830338/100‑1 (MQ=255)
                                gCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                                                     <  2:1244451/100‑1 (MQ=255)
                                 ccTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                                    <  1:818842/100‑1 (MQ=255)
                                    tATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                                                 <  1:2345365/100‑1 (MQ=255)
                                      tccttccT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                                               <  2:3884239/95‑1 (MQ=255)
                                        cGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                                             <  1:2476189/100‑1 (MQ=255)
                                                c‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGTGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttn                                                      <  2:1931116/100‑2 (MQ=255)
                                                c‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                                                      >  1:2366901/1‑100 (MQ=255)
                                                c‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                                                      <  1:2419988/100‑1 (MQ=255)
                                                  aaaaTCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                     <  1:1657849/100‑1 (MQ=255)
                                                   aaaTCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                                    <  1:4210/100‑1 (MQ=255)
                                                   aaaTCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                                    <  1:2193478/100‑1 (MQ=255)
                                                    aaTCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                   <  1:634616/100‑1 (MQ=255)
                                                         tGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  2:3824882/100‑1 (MQ=255)
                                                         tGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  2:2227926/100‑1 (MQ=255)
                                                         tGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  2:2086192/100‑1 (MQ=255)
                                                           cAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                            >  2:3373384/1‑100 (MQ=255)
                                                            aaaTTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                           <  1:1974633/100‑1 (MQ=255)
                                                            aaaTTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                           <  1:2705240/100‑1 (MQ=255)
                                                                 cAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                      <  1:1327043/100‑1 (MQ=255)
                                                                      cATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATGAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                 >  1:2393459/1‑100 (MQ=255)
                                                                       aTTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtttt                                >  1:3035376/1‑100 (MQ=255)
                                                                        ttG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                               >  2:3787056/1‑100 (MQ=255)
                                                                          g‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGATGTTTTCGa                             >  1:1576522/1‑100 (MQ=255)
                                                                          g‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             <  2:349183/100‑1 (MQ=255)
                                                                          g‑CA‑GAAT‑TGTC‑‑GTAGGCCTGATAAGCGGGGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCTATAGAGACGTTTTCGa                             <  2:127637/100‑1 (MQ=255)
                                                                            cA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                            >  1:3116603/1‑100 (MQ=255)
                                                                             a‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                           <  2:181879/100‑1 (MQ=255)
                                                                             a‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                           <  2:455336/100‑1 (MQ=255)
                                                                                  t‑tttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt                       >  1:3590499/1‑100 (MQ=255)
                                                                                          gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                  <  1:910787/100‑1 (MQ=255)
                                                                                          gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                  <  1:3404951/100‑1 (MQ=255)
                                                                                             gtccTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               <  2:2398520/98‑1 (MQ=255)
                                                                                             ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTTCTATCAGTCTTATATGGCTTTTAATTATTCGATAGAGACGTTTTCTAAGTGATGTTTGCCt               >  2:3201867/1‑100 (MQ=38)
                                                                                             ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               <  1:1172942/100‑1 (MQ=255)
                                                                                                 tGATAAGCGGAGCACATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           >  1:479381/1‑100 (MQ=255)
                                                                                                  gATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTCTTGCTTTTAATTCTTCGATAGAGACGTTTTCGAAGTGACGTTTGCCTAATgg          >  1:1885206/1‑100 (MQ=38)
                                                                                                    tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt        <  1:1490492/100‑1 (MQ=255)
                                                                                                      aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa      >  2:1152792/1‑100 (MQ=255)
                                                                                                       gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa     >  2:1975992/1‑100 (MQ=255)
                                                                                                        cGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAc    <  2:147475/100‑1 (MQ=255)
                                                                                                         ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTGATGGATCAAcc   <  1:3933044/100‑1 (MQ=255)
                                                                                                         ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc   <  1:101270/100‑1 (MQ=255)
                                                                                                          |                                                                                                   
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147

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: ‑