Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,226,991 T→G 9.8% intergenic (‑26/‑132) CJ019_04921 ← / → fecI_10 hypothetical protein/putative RNA polymerase sigma factor FecI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,226,9910TG 9.8% 149.7 / 2.7 60intergenic (‑26/‑132)CJ019_04921/fecI_10hypothetical protein/putative RNA polymerase sigma factor FecI
Reads supporting (aligned to +/- strand):  ref base T (20/33);  new base G (5/2);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACACTGGCCATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCTTACTACCAGTTGCGCTGTCATGCGCCGGAGTATCA  >  ACB122/5226858‑5227120
                                                                                                                                     |                                                                                                                                 
gACACTGGCCATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTg                                                                                                                                 <  7:306568/136‑1 (MQ=255)
gACACTGGCCATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTg                                                                                                                                 <  3:2006/136‑1 (MQ=255)
    cTGGCCATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCaa                                                                                                                              >  4:160725/1‑135 (MQ=255)
    cTGGCCATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCa                                                                                                                               >  6:363247/1‑134 (MQ=255)
    cTGGCCATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTGTTCGGGGTTGCa                                                                                                                               >  2:38325/1‑134 (MQ=255)
    cTGGCCATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTGTTCGGGGGTGcg                                                                                                                               >  1:129663/1‑133 (MQ=255)
      ggCTATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATa                                                                                                                            >  6:79208/1‑135 (MQ=255)
      ggCCATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATa                                                                                                                            >  6:302593/1‑135 (MQ=255)
        ccATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATa                                                                                                                            <  2:31384/133‑1 (MQ=255)
         cATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTaa                                                                                                                        >  4:263599/1‑136 (MQ=255)
         cATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTaa                                                                                                                        >  8:196903/1‑136 (MQ=255)
         cATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTa                                                                                                                         >  1:117352/1‑135 (MQ=255)
          aTCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGTCCAGGCTGAGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTAGCAATACTAAt                                                                                                                       >  4:161910/1‑136 (MQ=255)
          aTCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACTGGAGGCCGGCGGTTTTGAGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAAt                                                                                                                       >  3:337864/1‑136 (MQ=255)
          aTCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTaa                                                                                                                        >  6:325342/1‑135 (MQ=255)
          aTCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAAt                                                                                                                       >  4:98027/1‑136 (MQ=255)
          aTCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTGTTCGGGGTTGCAATACTaa                                                                                                                        >  1:232922/1‑135 (MQ=255)
          aTCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGCTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTGTTGGGGGTGGCAAGACTaa                                                                                                                        >  5:102783/1‑135 (MQ=255)
              aCGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGGTTTCGGGGGTGCAATACTAATGAtt                                                                                                                   >  4:120247/1‑136 (MQ=255)
                 ccAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCa                                                                                                                <  7:113290/136‑1 (MQ=255)
                      aaGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATtc                                                                                                            <  5:79208/135‑1 (MQ=255)
                      aaGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATtc                                                                                                            <  1:43031/135‑1 (MQ=255)
                            ccgccCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATaa                                                                                                      <  1:101487/135‑1 (MQ=255)
                               cccAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCa                                                                                                   >  4:173695/1‑135 (MQ=255)
                                ccAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCaa                                                                                                  <  5:312410/135‑1 (MQ=255)
                                         cgGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGa                                                                                        <  4:344593/136‑1 (MQ=255)
                                            ccAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGGAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGCCAg                                                                                     >  8:138475/1‑136 (MQ=255)
                                                     gTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTGAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGa                                                                              <  1:132795/134‑1 (MQ=255)
                                                      tGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGGGACGATTg                                                                           <  8:272642/136‑1 (MQ=255)
                                                       gAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGGTGCAATTCGAATGAGTTGAATTATGATAAGCAACGCCGGCGTAAAGGTGCAGAtttc                                                                          >  3:257344/1‑134 (MQ=255)
                                                            aaCGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTGGTCGGGGGTGGAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGCCAGGTGACGATTGCCGaaa                                                                     >  3:16224/1‑136 (MQ=255)
                                                               gCCCAGCCGGTAACTGTACGGGAGGCCGGCGGTTTTGGGCGTCATGGAGTAGGTCCTGGCGTTTTCGGGGTTTCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGGTGACAGGTGACGATTGCCGAAACAg                                                                  <  8:50044/136‑1 (MQ=255)
                                                               gCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAg                                                                  <  3:98027/136‑1 (MQ=255)
                                                                cccAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAg                                                                  <  1:70333/135‑1 (MQ=255)
                                                                cccAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAg                                                                  <  5:363247/135‑1 (MQ=255)
                                                                cccAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAg                                                                  <  5:97975/135‑1 (MQ=255)
                                                                 ccAGCCTGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGAc                                                                <  4:337864/136‑1 (MQ=255)
                                                                 ccAGCCGGTAAGTGAAGGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGGGTTTTGGGGGTTGCAATACTAATGATTTCTATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACTGAc                                                                <  8:390237/136‑1 (MQ=255)
                                                                 ccAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGAc                                                                <  7:196903/136‑1 (MQ=255)
                                                                 ccAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGAc                                                                <  7:26254/136‑1 (MQ=255)
                                                                  cAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGAc                                                                <  1:208960/135‑1 (MQ=255)
                                                                       ggTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGTAACAGACGGCCAc                                                          >  8:95622/1‑136 (MQ=255)
                                                                            cTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTGTTAGGGGTTTGAATACTAATGATTTCAATTATCAAAAGGAAATCCGGCTGACAGGTGTCGATTGGCGAAACAGACGGCCCCCCaaa                                                     >  7:388116/1‑134 (MQ=255)
                                                                              gaccGGGAGGCCGGCGGTTTTGCGCGTCATGGAGAAGGGCCGGGCGTTTTCGGGGGTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTTCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATAc                                                    <  1:323429/132‑1 (MQ=255)
                                                                               aaCGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATAc                                                    <  4:380914/134‑1 (MQ=255)
                                                                                aCGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGGCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATAc                                                    <  5:229041/133‑1 (MQ=255)
                                                                                          ggcggTTTTGCGCGTCATGGCGAAGGTACGGGCGTTTTAGGGGTTGCAATAATAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACGGACGGCCACCCATACAATGCGAACaa                                         >  6:287519/1‑134 (MQ=255)
                                                                                             ggTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCtt                                    <  3:134720/136‑1 (MQ=255)
                                                                                             ggTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCtt                                    <  8:26216/136‑1 (MQ=255)
                                                                                                      cgTCATGGCGAAGGTGCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCTTACTACCAg                            <  6:102783/135‑1 (MQ=255)
                                                                                                        tCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCTTACTACCAGtt                          <  6:135204/135‑1 (MQ=255)
                                                                                                                 aGGGCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCTTACTACCAGTTGCGCTGTCAt                <  3:144183/136‑1 (MQ=255)
                                                                                                                    tCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCTTACTACCAGTTGCGCTGTCATgcg             <  3:263599/136‑1 (MQ=255)
                                                                                                                    tCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCTTACTACCAGTTGCGCTGTCATgcg             >  8:140329/1‑136 (MQ=255)
                                                                                                                    gccGGGCGTTTTCGGGGGTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGGGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCTTACTACCAGTTGCGCTGTCATgcg             <  4:16224/135‑1 (MQ=255)
                                                                                                                     ccGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATACGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCTTACTACCAGTTGCGCTGTCATgcg             <  1:369205/135‑1 (MQ=255)
                                                                                                                     ccGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCTTACTACCAGTTGCGCTGTCATgcg             <  2:232922/135‑1 (MQ=255)
                                                                                                                       gggCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCTTACTACCAGTTGCGCTGTCATGCGcc           <  6:84838/135‑1 (MQ=255)
                                                                                                                               tCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCTTACTACCAGTTGCGCTGTCATGCGCCGGAGTATc   <  7:15602/135‑1 (MQ=255)
                                                                                                                                cGGGGGTGCAATACGAAAGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGGGACGCTTGCCGCAACAGACGGCCACCCATACAAAGCGAACAATTCTTACTACCCGTTGCGCTGTCATGCGCCGGAGTATCa  >  2:295728/1‑135 (MQ=255)
                                                                                                                                     |                                                                                                                                 
GACACTGGCCATCGACGCCAGCAAGTAGCCGCCCAGCAGCGCGGCCAGGCTGCGTGAGGCAACGCCCAGCCGGTAACTGAACGGGAGGCCGGCGGTTTTGCGCGTCATGGCGAAGGTCCGGGCGTTTTCGGGGTTGCAATACTAATGATTTCAATTCTCATAAGCAAAGCCTGCTGACAGGTGACGATTGCCGAAACAGACGGCCACCCATACAATGCGAACAATTCTTACTACCAGTTGCGCTGTCATGCGCCGGAGTATCA  >  ACB122/5226858‑5227120

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