Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,820,144 Δ1 bp intergenic (+24/+18) b4542 → / ← oraA hypothetical protein/OraA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,820,1430G.99.0% 419.9 / ‑3.0 95intergenic (+23/+19)b4542/oraAhypothetical protein/OraA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (61/33);  total (62/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 102 total aligned reads displayed.

CCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCAC  >  NC_000913/2820069‑2820214
                                                                          |                                                                       
 caGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgta                                                                       <  1:102206/76‑3 (MQ=255)
 caGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgta                                                                       <  1:5303740/76‑3 (MQ=255)
 caGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgta                                                                       <  2:2897182/76‑3 (MQ=255)
  aGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgtaa                                                                      >  1:1497595/1‑73 (MQ=255)
  aGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgtaa                                                                      <  1:2111950/76‑4 (MQ=255)
  aGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgtaa                                                                      >  2:1815274/1‑73 (MQ=255)
    cTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTaaaa                                                                   <  1:5646949/76‑1 (MQ=255)
     tCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAAt                                                                  <  1:4609226/76‑1 (MQ=255)
     tCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAAt                                                                  >  2:748212/1‑76 (MQ=255)
     tCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAAt                                                                  >  2:2102323/1‑76 (MQ=255)
     tCCAGCGTGTCTTAACCGCCGGGCTGGGAACTGAAAAGTGGGAATAAGATAAGCTTTCTTGACTGGGAAGTAAAAt                                                                  >  1:220698/1‑76 (MQ=255)
     tCCAACGTGTCTGAACAGCCGAGATGGGAAATGAAAAGTGGGAATAAGATAAGCTTTCTTGACTGGCAAGTaaaac                                                                  >  1:810043/1‑75 (MQ=255)
     cccAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAAt                                                                  <  2:5078866/75‑1 (MQ=255)
        aGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATAcc                                                               >  1:1085849/1‑76 (MQ=255)
        aGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATAcc                                                               <  1:1782490/76‑1 (MQ=255)
         gCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCg                                                              >  1:1724315/1‑76 (MQ=255)
          ggtgtCTAAGCCGCCGTTCTGGTAACGGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATCCCGt                                                             <  2:2215166/75‑1 (MQ=255)
          cGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATCAGTTATCTTGACTGGGAAGTAAAATACCGt                                                             <  1:1815274/76‑1 (MQ=255)
          cGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGt                                                             <  1:2628768/76‑1 (MQ=255)
          cGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGt                                                             >  2:459442/1‑76 (MQ=255)
          cGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGt                                                             <  2:2954139/76‑1 (MQ=255)
           gtgtCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTa                                                            >  2:375093/1‑76 (MQ=255)
           gtgtCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTa                                                            >  2:2877056/1‑76 (MQ=255)
            tgtCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTAt                                                           >  1:38823/1‑76 (MQ=255)
             gtCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATg                                                          <  1:222880/76‑1 (MQ=255)
              ncTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTATTGACTGGGAAGTACAATACCGTATga                                                         >  1:4346832/2‑75 (MQ=255)
               cTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCg                                                        <  1:5092124/76‑1 (MQ=255)
               cTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCg                                                        >  1:1121306/1‑76 (MQ=255)
               cTTAACCGCCGGGCTGGCAACTGAAAAGTGGGAATAAGATAAGCTTTCTTGACTGGGAAGTAACATACCGTATGCg                                                        >  2:4499672/1‑76 (MQ=255)
               cTTAACCGCCGGGCTGGAAACTGAAAAGTGCGAATAAGATAAGCTTTCTTGACTCGGAAGTAAAATACCGTATGCg                                                        >  2:4514408/1‑76 (MQ=255)
                 tAACCGCCGGGCCGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGtt                                                      <  1:1261665/76‑1 (MQ=255)
                 tAACCGCCCGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGtt                                                      >  2:2036907/1‑76 (MQ=255)
                    ccgccgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAg                                                   <  2:20266/76‑1 (MQ=255)
                     cgccgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAgc                                                  >  2:2615361/1‑75 (MQ=255)
                      gccgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTc                                                 >  1:3371076/1‑76 (MQ=255)
                      gccgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTc                                                 >  1:2099255/1‑76 (MQ=255)
                      gccgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTc                                                 <  2:2021296/76‑1 (MQ=255)
                        cgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCgg                                               >  2:5000716/1‑76 (MQ=255)
                        cgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCgg                                               >  2:3467490/1‑76 (MQ=255)
                         gggCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGc                                              >  2:1324167/1‑76 (MQ=255)
                         gggCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGc                                              >  1:3205301/1‑76 (MQ=255)
                         gggCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGc                                              >  2:2407045/1‑76 (MQ=255)
                         gggCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGc                                              >  2:736390/1‑76 (MQ=255)
                             tGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCaaaa                                          <  1:5534843/76‑1 (MQ=255)
                              ggTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAAt                                         >  1:774894/1‑76 (MQ=255)
                              ggTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAAt                                         >  2:5110599/1‑76 (MQ=255)
                                tAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAAttt                                       >  1:3979575/1‑76 (MQ=255)
                                  aCTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCg                                     >  1:267523/1‑76 (MQ=255)
                                  aCTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCg                                     >  1:1968125/1‑76 (MQ=255)
                                  aCTGAAAAGTGGGAAGAAGCTAAGTTTTCTTGACTGTGAAGTAAAATACCGTCTGCGTTCCGTCGACAAAAattac                                     >  1:295186/1‑74 (MQ=255)
                                   cTGAAAAGTGGGACTAACATAAGTTTTCTTGACTGGGCAGTAACATACCGTATGCGTTCAGTCGGCCAAATTTCGc                                    <  1:3853727/76‑1 (MQ=255)
                                   cTGAAAAGGGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGc                                    >  1:4811332/1‑76 (MQ=255)
                                        aaGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAt                               >  2:2625458/1‑76 (MQ=255)
                                        aaGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAt                               >  1:3279933/1‑76 (MQ=255)
                                        aaGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAt                               >  2:3683244/1‑76 (MQ=255)
                                         aGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAtc                              <  2:3961058/76‑1 (MQ=255)
                                          gTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAtct                             <  2:3469952/76‑1 (MQ=255)
                                           tGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAtctc                            >  2:2553826/1‑76 (MQ=255)
                                           tGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAtctc                            >  1:141297/1‑76 (MQ=255)
                                            gggAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTcc                           <  1:5683535/76‑1 (MQ=255)
                                             ggAATAAGATAAGTTTTCTTGACTGGGAAGTAAAGTACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAcctcct                          >  2:5004064/1‑76 (MQ=255)
                                             ggAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCt                          >  2:5125625/1‑76 (MQ=255)
                                               aaTAAGATAAGGTTTCTTGACTGGGAAGCAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTgg                        >  2:1560205/1‑76 (MQ=255)
                                                 taagataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGat                      <  1:1987450/76‑1 (MQ=255)
                                                 taagataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGat                      >  1:2650900/1‑76 (MQ=255)
                                                  aagataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGata                     <  1:32144/76‑1 (MQ=255)
                                                  aagataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGata                     <  1:4957536/76‑1 (MQ=255)
                                                  aagataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGata                     >  2:2911820/1‑76 (MQ=255)
                                                   agataagTTTTCTTGACTGGGAAGTAAAATACCGTATGTGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGatat                    <  1:3381155/76‑1 (MQ=255)
                                                    gataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATc                   >  2:1805155/1‑76 (MQ=255)
                                                    gataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATc                   >  1:4409443/1‑76 (MQ=255)
                                                     ataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCt                  >  2:4019093/1‑76 (MQ=255)
                                                     ataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCt                  >  1:5697458/1‑76 (MQ=255)
                                                        cgTTTTCTTGACTGGGAAGTCAAATACCGTATGCGTTCAGTCGGCACAATTTCGCCAAATCTCCTGGATATCTTcc               <  2:501538/75‑1 (MQ=255)
                                                         gTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCa              <  2:859453/76‑1 (MQ=255)
                                                         gTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCa              >  1:3086440/1‑76 (MQ=255)
                                                          ttttCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCAt             <  1:2859746/76‑1 (MQ=255)
                                                           tttCTTGACTGGGAAGGAAAATACAGTATGAGTTCAGGCGGAAAAATTGAGCCAAATCGCCTGGGTAGATTCCAta            >  1:1043701/1‑75 (MQ=255)
                                                           tttCTTGACTGGGAAGTAAAATACCGTATGCGTTCGGTCGGGAAAATTTCGCCAAATCTCCTGGATATCTTCCATc            <  2:2644550/76‑1 (MQ=255)
                                                           tttCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATc            >  2:4963080/1‑76 (MQ=255)
                                                           tttCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATc            >  1:1935510/1‑76 (MQ=255)
                                                           tttCTTAACTGGGAAGTAAAATACCGGATGCGGTCAGTCGGCAAAATGTCGACAAATCTTCTGGATATCGTCCATc            >  1:2710044/1‑76 (MQ=255)
                                                            ttCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCa           >  2:4341606/1‑76 (MQ=255)
                                                              cTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGa         <  1:5256091/76‑1 (MQ=255)
                                                               ttGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGCATATCTTCCAGCAGAt        >  2:3560497/1‑76 (MQ=255)
                                                                tGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATa       <  2:1724315/76‑1 (MQ=255)
                                                                tGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATa       >  1:1224544/1‑76 (MQ=255)
                                                                tGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATa       >  1:1365777/1‑76 (MQ=255)
                                                                tGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATa       <  1:1456678/76‑1 (MQ=255)
                                                                tGACTGGGAAGTAAAATACCGTATGCGCTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATa       >  1:2535130/1‑76 (MQ=255)
                                                                 gACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAg      <  2:3135208/76‑1 (MQ=255)
                                                                   cTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGcc    <  1:2911820/76‑1 (MQ=255)
                                                                   cTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGcc    >  2:5007414/1‑76 (MQ=255)
                                                                   cTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGcc    >  1:5380562/1‑76 (MQ=255)
                                                                    tGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCa   <  1:4213867/76‑1 (MQ=255)
                                                                    tGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCa   <  1:5196901/76‑1 (MQ=255)
                                                                    tGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCa   >  2:5400632/1‑76 (MQ=255)
                                                                    tGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCa   >  2:163342/1‑76 (MQ=255)
                                                                     gggAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCAc  >  1:1143798/1‑76 (MQ=255)
                                                                     gggAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCAc  >  2:867511/1‑76 (MQ=255)
                                                                          |                                                                       
CCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCAC  >  NC_000913/2820069‑2820214

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