Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,970 G→A 29.2% intergenic (‑180/+128) 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,9700GA29.2% 362.7 / 156.2 200intergenic (‑180/+128)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (66/75);  new base A (48/10);  total (115/85)
Fisher's exact test for biased strand distribution p-value = 3.24e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-01

Only 100 of 211 total aligned reads displayed.

CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069
                                                                                                   |                                                                                                          
      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                      <  2:2643177/100‑1 (MQ=255)
      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                      >  1:101470/1‑100 (MQ=255)
      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAATGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                      >  2:2279716/1‑100 (MQ=255)
       tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                     <  2:684624/100‑1 (MQ=255)
        gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                    >  1:1881753/1‑100 (MQ=255)
        gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                    >  1:1489292/1‑100 (MQ=255)
           gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                 <  1:3487683/100‑1 (MQ=255)
            tATTGCTGCCGGATGGCGCTGCACGCAGAGTCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTCATGCCGGATGCGGCGTTAACGc                                                                                                <  2:2110316/100‑1 (MQ=255)
            tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                <  1:2323024/100‑1 (MQ=255)
            tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                <  2:290655/100‑1 (MQ=255)
                gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                            <  1:1709789/100‑1 (MQ=255)
                gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                            >  1:2116179/1‑100 (MQ=255)
                 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                           <  2:494654/100‑1 (MQ=255)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                      <  1:2272090/100‑1 (MQ=255)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                      >  2:1627169/1‑100 (MQ=255)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                      <  1:3293593/100‑1 (MQ=37)
                        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGATTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                    <  2:1992241/100‑1 (MQ=255)
                         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑a                                                                                 <  1:213316/100‑1 (MQ=255)
                          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑Ac                                                                                >  1:161693/1‑100 (MQ=255)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑Ac                                                                                >  2:3460089/1‑99 (MQ=255)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACa                                                                              <  2:1529982/100‑1 (MQ=255)
                            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                                <  2:3415184/100‑7 (MQ=25)
                                gCACGCAGAGGAACAGCGCCACCGGGTATACAGATAGTGGACAAAGTGCGCTTTGTTTATGCCGAATAAGGCGTTAACGCCGCATCCGGcataaacaaa                                                                             <  2:1717787/99‑10 (MQ=18)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCg                                                                        <  1:3044625/100‑1 (MQ=255)
                                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGt                                                                       >  1:342496/1‑100 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTg                                                                      <  1:1370036/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTg                                                                      >  1:3021453/1‑100 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGc                                                                     <  1:1974731/100‑1 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCa                                                                    <  1:206555/100‑1 (MQ=255)
                                     cAGACGCACAGCGCCACCGGCTATACAGATTGTTGACAAATTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGc                                                                     <  2:1077484/99‑1 (MQ=255)
                                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAAt                                                                 <  1:628829/100‑1 (MQ=255)
                                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAAtt                                                                >  2:101148/1‑100 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCAGCCT‑‑ACAAAATCGTGCAAATTc                                                               >  2:380018/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCa                                                              >  1:3080347/1‑100 (MQ=255)
                                            acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCaa                                                             >  1:1794866/1‑100 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAAt                                                            <  2:310388/100‑1 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAAt                                                            <  2:671520/100‑1 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCTGCCT‑‑ACAAAATCGTGCAAATTCAATa                                                           <  1:2219017/100‑1 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATa                                                           >  1:740880/1‑100 (MQ=255)
                                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATAc                                                          >  1:2019660/1‑100 (MQ=255)
                                                cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACa                                                         >  2:3128333/1‑100 (MQ=255)
                                                 gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACAt                                                        >  1:134892/1‑100 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACAtt                                                       >  1:950880/1‑100 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACAtt                                                       <  1:2493535/100‑1 (MQ=255)
                                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTg                                                      <  1:3361132/100‑1 (MQ=255)
                                                     ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGc                                                     <  2:269972/99‑1 (MQ=255)
                                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTc                                          <  1:317356/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAANTGCGNTTTGTTNATGCCGGNTGCGGCGTTAACGNCTTNTCCGGCCT‑‑ACAAAATCGTGCAAATNCAANACATTGCAGAATTTTc                                          >  2:2221156/1‑99 (MQ=17)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                       <  1:1469346/100‑1 (MQ=255)
                                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTa                                     >  2:3064913/1‑100 (MQ=255)
                                                                  gcataaaCAAAGCGCACTTCGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAgg                                   >  2:1160602/7‑100 (MQ=18)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGcc                                 <  1:2983164/100‑1 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGcc                                 <  2:740877/100‑1 (MQ=255)
                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATa                            >  1:93419/1‑100 (MQ=255)
                                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGAAGAATTTTC‑‑GTGGGACTGATatg                          >  2:3579623/1‑97 (MQ=17)
                                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCg                        >  1:115850/1‑100 (MQ=255)
                                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCg                        <  2:1865346/100‑1 (MQ=255)
                                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCgg                       <  1:2445447/100‑1 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAGTCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAg                     <  2:1587848/99‑1 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                    <  2:617577/100‑1 (MQ=255)
                                                                                  tccgatctaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAgcgc                  <  1:1140327/90‑1 (MQ=18)
                                                                                  gaatttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCT‑‑ACAAATTCTTCCAAATTCAGTATATTGCAGAAATCAT‑‑GTAGGCCTGATAAGCGTAgcg                   >  1:1752813/11‑100 (MQ=12)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                   >  1:446404/1‑100 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                   >  2:113461/1‑100 (MQ=255)
                                                                                  agttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAAAACTGTGCAAATTCAATATATTGCAGGAAACAC‑‑GTAGGCCTGATAAGCGAAg                     >  2:754061/11‑100 (MQ=2)
                                                                                   tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcgc                  >  2:328768/1‑100 (MQ=255)
                                                                                     tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGTGTAGCGCAt                >  1:1271158/8‑100 (MQ=18)
                                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa              >  1:1692890/1‑100 (MQ=255)
                                                                                       tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑‑ACAAGTTCGTGCAAATTCAATAAATTGCAATATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  2:1696765/6‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  1:2192925/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAACAGTGCAAATTCAATATATTGCAATCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAg             >  2:2763245/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  1:2535929/5‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg            >  1:2786673/4‑100 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCTGGCCT‑‑ACAAATTCTTCCAAATTCAATATATTGCAGAAATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:1822416/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAACATATTGCAATTCTCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:861865/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:1847401/3‑100 (MQ=17)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAACAGTGCAAATTCAATATATTGCAATCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc           >  1:1833242/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACTCCGCATCCGGCCT‑‑ACAAAAGCACGCAAATTCAATATATTGCAGAGATTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:3564957/3‑100 (MQ=11)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACCTTGCTAATTCAATATATTGCAGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc           >  1:549715/3‑100 (MQ=18)
                                                                                          aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATGCGCTGCGCTTATCAGGc           <  2:653849/100‑1 (MQ=21)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  1:3553624/1‑100 (MQ=255)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCGT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:2451247/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt         >  1:2154372/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:2813960/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCATGAAAATTCAATACATTGCAAGATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:145228/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:2463804/2‑99 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:13546/2‑100 (MQ=18)
                                                                                            gCCGGATCAGGCGTTAACGCCGCATCCGGCCT‑‑ACAAAGTCGAGCAAATTCACTACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGATCGCATCAGGCa          >  2:2439315/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  1:1986198/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACGAAATCGTGCTAATTCAATATATTGTGGGCTCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgaga          >  1:287876/1‑96 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAGTTCGTGCAAATTCAATAAATTGCAATATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat        >  2:582636/1‑97 (MQ=21)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCATGAAAATTCAATACATTGCAAGATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCac         >  1:269161/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATCTCCGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  2:2319068/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAATCATGCAAATCCAAAATATTGCAGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt        >  1:1666181/1‑100 (MQ=21)
                                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt        >  1:1554762/1‑100 (MQ=255)
                                                                                             ccGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACagt        >  1:3240046/1‑98 (MQ=17)
                                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt       >  1:2927220/1‑100 (MQ=255)
                                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  1:1065965/100‑1 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt     >  2:3405976/1‑97 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc  <  2:1505022/100‑1 (MQ=255)
                                                                                                   |                                                                                                          
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069

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: