Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,571,178 G→C 5.3% R361R (CGG→CGC)  yibH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,571,1780GC 5.3% 306.6 / 4.8 94R361R (CGG→CGC) yibHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (44/45);  new base C (2/3);  total (46/48)
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.08e-01

CTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAT  >  Exported/2571109‑2571248
                                                                     |                                                                      
cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                                       <  1:3002484/71‑1 (MQ=255)
  tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCgg                                                                        >  1:3681008/1‑68 (MQ=255)
   aCCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaaa                                                                     <  1:3466162/70‑1 (MQ=255)
    cccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa                                                                      <  1:3303731/68‑1 (MQ=255)
    cccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGt                                                                   <  1:679108/71‑1 (MQ=255)
     ccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtg                                                                  <  1:2989367/71‑1 (MQ=255)
      cGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgc                                                                 <  1:2982458/71‑1 (MQ=255)
       gACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaaa                                                                     >  1:2365831/1‑66 (MQ=255)
       gACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgct                                                                >  1:2541110/1‑71 (MQ=255)
        aCGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgctg                                                               <  1:845885/71‑1 (MQ=255)
        aCGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgct                                                                >  1:1965061/1‑70 (MQ=255)
        aCGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgc                                                                 <  1:656072/69‑1 (MQ=255)
         cGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgctgc                                                              >  1:2875994/1‑71 (MQ=255)
         cGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgctgc                                                              >  1:2705032/1‑71 (MQ=255)
          ggCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgct                                                                >  1:2781492/1‑68 (MQ=255)
            cATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct                                                             >  1:1883652/1‑69 (MQ=255)
            cATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct                                                             >  1:849546/1‑69 (MQ=255)
            cATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                                           <  1:911040/71‑1 (MQ=255)
            cATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                                           <  1:1094787/71‑1 (MQ=255)
            cATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                                           <  1:1133596/71‑1 (MQ=255)
               cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGATAAGaa                                                        >  1:295521/1‑71 (MQ=255)
                 aCGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa                                                                      <  1:3390957/55‑1 (MQ=255)
                 aCGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                                       >  1:3214490/1‑70 (MQ=255)
                 aCGCCCAGGTGGCGGTTTACTCAGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                        >  1:619580/1‑69 (MQ=255)
                  cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                        <  1:2378183/68‑1 (MQ=255)
                  cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                                                         <  1:2791841/67‑1 (MQ=255)
                  cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                                     <  1:108415/71‑1 (MQ=255)
                   gCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAg                                                                    >  1:419597/1‑55 (MQ=255)
                   gCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAc                                                    <  1:3565690/71‑1 (MQ=255)
                    cccAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAc                                                    >  1:2782610/1‑70 (MQ=255)
                     ccAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGAc                                                    >  1:1570520/1‑69 (MQ=255)
                      cAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGt                                                                   >  1:1414710/1‑53 (MQ=255)
                      cAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAg                                                 <  1:594198/71‑1 (MQ=255)
                      cAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAg                                                 <  1:3737795/71‑1 (MQ=255)
                       aGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtg                                                                  <  1:1215904/53‑1 (MQ=255)
                       aGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtg                                                                  <  1:3546402/53‑1 (MQ=255)
                       aGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtg                                                                  <  1:626346/53‑1 (MQ=255)
                       aGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgctg                                                               >  1:2395968/1‑56 (MQ=255)
                       aGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGc                                                <  1:610570/71‑1 (MQ=255)
                        ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgc                                                                 <  1:3712722/53‑1 (MQ=255)
                         gTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgct                                                                >  1:1151875/1‑53 (MQ=255)
                          tGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTgg                                             <  1:2285061/71‑1 (MQ=255)
                           ggCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                            >  1:3401522/1‑71 (MQ=255)
                              ggTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATg                                          >  1:1734108/1‑70 (MQ=255)
                                tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                                           >  1:2166353/1‑51 (MQ=255)
                                tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                       <  1:3574804/71‑1 (MQ=255)
                                tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                       >  1:1907160/1‑71 (MQ=255)
                                tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                       >  1:1845279/1‑71 (MQ=255)
                                tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                       <  1:1830674/71‑1 (MQ=255)
                                 ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                       >  1:450177/1‑70 (MQ=255)
                                 ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                       >  1:2069407/1‑70 (MQ=255)
                                 ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                       >  1:219965/1‑70 (MQ=255)
                                  tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                        >  1:709563/1‑52 (MQ=255)
                                  tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt                                      <  1:1355367/70‑1 (MQ=255)
                                  tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt                                      <  1:301701/70‑1 (MQ=255)
                                  tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt                                      <  1:3172591/70‑1 (MQ=255)
                                   aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTAt                                    >  1:223161/1‑71 (MQ=255)
                                     tCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtg                                                                  <  1:3621087/39‑1 (MQ=255)
                                       cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCt                                  <  1:3068778/69‑1 (MQ=255)
                                         aCCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAt                              <  1:463010/71‑1 (MQ=255)
                                          ccATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAtt                             >  1:223946/1‑71 (MQ=255)
                                           cATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAttt                            <  1:2292099/71‑1 (MQ=255)
                                           cATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAttt                            <  1:2421005/71‑1 (MQ=255)
                                              ttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                              <  1:3234450/50‑1 (MQ=255)
                                              ttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                              <  1:2435110/50‑1 (MQ=255)
                                                cAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                          >  1:970297/1‑36 (MQ=255)
                                                cAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                          >  1:3145693/1‑36 (MQ=255)
                                                cAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                          >  1:1210522/1‑36 (MQ=255)
                                                cAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                          >  1:2631515/1‑36 (MQ=255)
                                                cAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                        >  1:1251983/1‑70 (MQ=255)
                                                    cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                   <  1:2262314/71‑1 (MQ=255)
                                                     aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAcc                                                   >  1:3072032/1‑38 (MQ=255)
                                                     aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAt                              <  1:2325307/59‑1 (MQ=255)
                                                     aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                    >  1:935101/1‑69 (MQ=255)
                                                     aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                    >  1:181178/1‑69 (MQ=255)
                                                     aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                  <  1:3363742/71‑1 (MQ=255)
                                                     aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                  <  1:2841017/71‑1 (MQ=255)
                                                     aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                  <  1:2756632/71‑1 (MQ=255)
                                                     aTGTTTCGGTGACGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                  <  1:799871/71‑1 (MQ=255)
                                                      tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgag                 <  1:2031544/71‑1 (MQ=255)
                                                       gTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagg                >  1:409872/1‑71 (MQ=255)
                                                       gTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATAATTgagg                >  1:1648363/1‑71 (MQ=255)
                                                        tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggg               <  1:3096963/71‑1 (MQ=255)
                                                        tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagg                <  1:3398433/70‑1 (MQ=255)
                                                         ttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCt                                  <  1:2646435/51‑1 (MQ=255)
                                                          tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCttt                                >  1:489480/1‑52 (MQ=255)
                                                            ggTGATGCGGAAAGTGCTGCTAAGAATGAACAGCTGGATGCATTATCTTTATTTGGATCATTgagggggag           <  1:2521871/71‑1 (MQ=255)
                                                              tGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggggaggg         >  1:847273/1‑71 (MQ=255)
                                                              tGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggggagg          <  1:296740/70‑1 (MQ=255)
                                                                 tGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATc                       >  1:121726/1‑54 (MQ=255)
                                                                 tGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATc                       >  1:3400938/1‑54 (MQ=255)
                                                                 tGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATc                       >  1:184713/1‑54 (MQ=255)
                                                                  gCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGagggggaggggg       >  1:1914582/1‑69 (MQ=255)
                                                                   cGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg              >  1:518087/1‑61 (MQ=255)
                                                                   cGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATATTTATTTGGATCATTgagggg              >  1:1325974/1‑61 (MQ=255)
                                                                     gAAAGTACTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt  >  1:2584/1‑71 (MQ=255)
                                                                     |                                                                      
CTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAT  >  Exported/2571109‑2571248

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