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‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑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‑‑AC‑a                                                                              <  2:1529982/100‑1 (MQ=255)
                            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                                <  2:3415184/100‑7 (MQ=25)
                                gCACGCAGAGGAACAGCGCCACCGGGTATACAGATAGTGGACAAAGTGCGCTTTGTTTATGCCGAATAAGGCGTTAACGCCGCATCCGGcataaacaaa                                                                             <  2:1717787/99‑10 (MQ=18)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCg                                                                        <  1:3044625/100‑1 (MQ=255)
                                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGt                                                                       >  1:342496/1‑100 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                      <  1:1370036/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                      >  1:3021453/1‑100 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                     <  1:1974731/100‑1 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCa                                                                    <  1:206555/100‑1 (MQ=255)
                                     cAGACGCACAGCGCCACCGGCTATACAGATTGTTGACAAATTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                     <  2:1077484/99‑1 (MQ=255)
                                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAAt                                                                 <  1:628829/100‑1 (MQ=255)
                                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAAtt                                                                >  2:101148/1‑100 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCAGCCT‑‑AC‑AAAATCGTGCAAATTc                                                               >  2:380018/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCa                                                              >  1:3080347/1‑100 (MQ=255)
                                            acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCaa                                                             >  1:1794866/1‑100 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                            <  2:310388/100‑1 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                            <  2:671520/100‑1 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCTGCCT‑‑AC‑AAAATCGTGCAAATTCAATa                                                           <  1:2219017/100‑1 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATa                                                           >  1:740880/1‑100 (MQ=255)
                                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATAc                                                          >  1:2019660/1‑100 (MQ=255)
                                                cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACa                                                         >  2:3128333/1‑100 (MQ=255)
                                                 gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAt                                                        >  1:134892/1‑100 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAtt                                                       >  1:950880/1‑100 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAtt                                                       <  1:2493535/100‑1 (MQ=255)
                                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTg                                                      <  1:3361132/100‑1 (MQ=255)
                                                     ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGc                                                     <  2:269972/99‑1 (MQ=255)
                                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTc                                          <  1:317356/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAANTGCGNTTTGTTNATGCCGGNTGCGGCGTTAACGNCTTNTCCGGCCT‑‑AC‑AAAATCGTGCAAATNCAANACATTGCA‑GAAT‑TTTc                                          >  2:2221156/1‑99 (MQ=17)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑g                                       <  1:1469346/100‑1 (MQ=255)
                                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTa                                     >  2:3064913/1‑100 (MQ=255)
                                                                  gcataaaCAAAGCGCACTTCGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAgg                                   >  2:1160602/7‑100 (MQ=18)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGcc                                 <  1:2983164/100‑1 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGcc                                 <  2:740877/100‑1 (MQ=255)
                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATa                            >  1:93419/1‑100 (MQ=255)
                                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGAA‑GAAT‑TTTC‑‑GTGGGACTGATatg                          >  2:3579623/1‑97 (MQ=17)
                                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCg                        >  1:115850/1‑100 (MQ=255)
                                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCg                        <  2:1865346/100‑1 (MQ=255)
                                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCgg                       <  1:2445447/100‑1 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAGTCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                     <  2:1587848/99‑1 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                    <  2:617577/100‑1 (MQ=255)
                                                                                  tccgatctaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAgcgc                  <  1:1140327/90‑1 (MQ=18)
                                                                                  gaatttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCT‑‑AC‑AAATTCTTCCAAATTCAGTATATTGCA‑GAAA‑TCAT‑‑GTAGGCCTGATAAGCGTAgcg                   >  1:1752813/11‑100 (MQ=12)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                   >  1:446404/1‑100 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                   >  2:113461/1‑100 (MQ=255)
                                                                                  agttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAAAC‑AAAACTGTGCAAATTCAATATATTGCA‑GGAA‑ACAC‑‑GTAGGCCTGATAAGCGAAg                     >  2:754061/11‑100 (MQ=2)
                                                                                   tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcgc                  >  2:328768/1‑100 (MQ=255)
                                                                                     tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGTGTAGCGCAt                >  1:1271158/8‑100 (MQ=18)
                                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa              >  1:1692890/1‑100 (MQ=255)
                                                                                       tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑‑AC‑AAGTTCGTGCAAATTCAATAAATTGCA‑ATA‑‑TGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  2:1696765/6‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  1:2192925/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑ATCT‑TCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAg             >  2:2763245/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  1:2535929/5‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg            >  1:2786673/4‑100 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCTGGCCT‑‑AC‑AAATTCTTCCAAATTCAATATATTGCA‑GAAA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:1822416/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAACATATTGCA‑‑ATTCTCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:861865/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:1847401/3‑100 (MQ=17)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑ATCT‑TCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc           >  1:1833242/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACTCCGCATCCGGCCT‑‑AC‑AAAAGCACGCAAATTCAATATATTGCA‑GAGA‑TTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:3564957/3‑100 (MQ=11)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACCTTGCTAATTCAATATATTGCAGGGAC‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc           >  1:549715/3‑100 (MQ=18)
                                                                                          aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTG‑‑GTAGGCCTGATGCGCTGCGCTTATCAGGc           <  2:653849/100‑1 (MQ=21)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  1:3553624/1‑100 (MQ=255)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCGT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:2451247/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt         >  1:2154372/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:2813960/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCATGAAAATTCAATACATTGCA‑AGAT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:145228/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:2463804/2‑99 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:13546/2‑100 (MQ=18)
                                                                                            gCCGGATCAGGCGTTAACGCCGCATCCGGCCT‑‑AC‑AAAGTCGAGCAAATTCACTACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGATCGCATCAGGCa          >  2:2439315/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  1:1986198/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgaga          >  1:287876/1‑96 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑ATA‑TGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat        >  2:582636/1‑97 (MQ=21)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCATGAAAATTCAATACATTGCA‑AGAT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCac         >  1:269161/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGCA‑ATCT‑CCGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  2:2319068/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt        >  1:1666181/1‑100 (MQ=21)
                                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt        >  1:1554762/1‑100 (MQ=255)
                                                                                             ccGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACagt        >  1:3240046/1‑98 (MQ=17)
                                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt       >  1:2927220/1‑100 (MQ=255)
                                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  1:1065965/100‑1 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt     >  2:3405976/1‑97 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc  <  2:1505022/100‑1 (MQ=255)
                                                                                                   |                                                                                                          
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTC‑‑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: ‑