Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,453,506 C→T 13.6% intergenic (‑17/+96) PP_3069 ← / ← ppiC‑I outer membrane autotransporter/peptidyl‑prolyl cis‑trans isomerase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,453,5060CT13.6% 107.9 / 8.6 44intergenic (‑17/+96)PP_3069/ppiC‑Iouter membrane autotransporter/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base C (19/19);  new base T (4/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

CACAGGTGTAGGTGTCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATC  >  NC_002947/3453357‑3453653
                                                                                                                                                     |                                                                                                                                                   
cacaGGTGTAGGTGTCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTgc                                                                                                                                                     <  1:526590/150‑1 (MQ=255)
 acaGGTGTAGGTGTCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTgcg                                                                                                                                                    >  1:614747/1‑150 (MQ=255)
     gtgtaggtgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGcacc                                                                                                                                                >  1:161265/1‑150 (MQ=255)
      tgtaggtgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGcacca                                                                                                                                               <  2:555908/150‑1 (MQ=255)
          ggtgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGTGCACCAcgc                                                                                                                                            >  2:1191484/1‑149 (MQ=255)
          ggtgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGTGCACCAcgc                                                                                                                                            <  1:1191484/149‑1 (MQ=255)
           gtgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGATCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgt                                                                                                                                          >  1:799211/1‑150 (MQ=255)
            tgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCgtg                                                                                                                                         <  2:737323/150‑1 (MQ=255)
            tgtCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATAAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGTGCACCACGCgtg                                                                                                                                         >  1:533074/1‑150 (MQ=255)
                    ccggggccacgTGAGAAGGTACAAGCCGCGTGTGCCGCTGGCCATGCGATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAg                                                                                                                                 <  2:11434/147‑1 (MQ=255)
                    cccGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCGCAg                                                                                                                                 <  1:454681/150‑1 (MQ=255)
                    cccGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATAAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGTGCACCACGCGTGGTGCGCAg                                                                                                                                 <  2:533074/150‑1 (MQ=255)
                                    aGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGta                                                                                                                 <  2:22089/150‑1 (MQ=255)
                                          aaGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGcgca                                                                                                                                  >  2:57888/1‑124 (MQ=255)
                                          aaGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGcgca                                                                                                                                  <  1:57888/127‑1 (MQ=255)
                                             ccGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGt                                                                                                        <  1:873762/150‑1 (MQ=255)
                                               gcgTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGc                                                                                                      <  1:775687/150‑1 (MQ=255)
                                                      cAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTCGCCGCACCAGGAACGGCTCCAGTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGa                                                                                               <  1:448592/150‑1 (MQ=255)
                                                        gCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGAGGTTCCATAGGGGCCCGCGTCACGAcc                                                                                             >  2:69181/1‑150 (MQ=255)
                                                                              gcGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGt                                                                       >  2:418818/1‑150 (MQ=255)
                                                                               cGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGAACACAACTGGCAAAAATGTCAGGTg                                                                      <  1:342437/150‑1 (MQ=255)
                                                                                         cccAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCt                                                            >  1:420163/1‑150 (MQ=255)
                                                                                          ccAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTa                                                           >  1:664470/1‑150 (MQ=255)
                                                                                          ccAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTa                                                           >  1:650562/1‑150 (MQ=255)
                                                                                                tCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGt                                                     <  2:420163/150‑1 (MQ=255)
                                                                                                   gCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAg                                                  >  2:991521/1‑150 (MQ=255)
                                                                                                        aGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTgc                                                                                                                                                     >  2:905198/1‑46 (MQ=255)
                                                                                                        aGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTgc                                                                                                                                                     <  1:905198/46‑1 (MQ=255)
                                                                                                        aGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTcgcg                                             >  2:44642/1‑150 (MQ=255)
                                                                                                          gCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGt                                                                                                                  <  2:809074/79‑1 (MQ=255)
                                                                                                          gCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGt                                                                                                                  >  1:809074/1‑79 (MQ=255)
                                                                                                           cGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTa                                          <  1:991521/150‑1 (MQ=255)
                                                                                                                 cGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGt                                    <  1:69181/150‑1 (MQ=255)
                                                                                                                  gCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCggg                                                                                                                       >  2:145224/1‑66 (MQ=255)
                                                                                                                  gCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCggg                                                                                                                       <  1:145224/66‑1 (MQ=255)
                                                                                                                   cATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGttt                                  <  2:650562/150‑1 (MQ=255)
                                                                                                                    aTCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTg                                 >  2:963528/1‑150 (MQ=255)
                                                                                                                               cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGt                                                     <  1:349586/119‑1 (MQ=255)
                                                                                                                               cTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGt                                                     >  2:349586/1‑119 (MQ=255)
                                                                                                                                    ttCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCg                 >  2:1144313/1‑150 (MQ=255)
                                                                                                                                        ccAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAAtt             >  1:437142/1‑150 (MQ=255)
                                                                                                                                         cAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAAttt            >  2:588223/1‑150 (MQ=255)
                                                                                                                                                 ccTGTGCACCACGCGTGGTGCACAGGGTTGCCGGGGATCTATAGAGGAGTCCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGa    >  2:920624/1‑150 (MQ=255)
                                                                                                                                                 ccTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGa    <  1:44642/150‑1 (MQ=255)
                                                                                                                                                  cTGTGCACCACGCGTGGTGCGCAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTa                     >  1:639058/1‑132 (MQ=255)
                                                                                                                                                   tGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTa                                          >  1:963460/1‑110 (MQ=255)
                                                                                                                                                   tGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTa                                          <  2:963460/110‑1 (MQ=255)
                                                                                                                                                   tGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATc  >  1:137893/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CACAGGTGTAGGTGTCGTTACCCGGGCCGGGTGAGAAGGTACAAGCCGCGTGTGCAGCGGGCGAGGCCATGCCCAGCAGCGGGAGCGAACCCAGGGTCAGCATCAGGCGTAGACGCATCAGGAAAACCTCCATTCACCAGGGAAACCTGCGCACCACGCGTGGTGCACAGGGTTACCGGGGTTGTATAGAGGAGTGCGTCACGACCACAACTGGCAAAAATGTCAGGTGGCAAAGGGCTAGCGCGTTAGTCGCGGTAGAACGTTTGCACCAGGTGGTAGCCGAATTTGCTCTTGATC  >  NC_002947/3453357‑3453653

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: ‑