Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,010,375 C→G 6.8% R82P (CGC→CCC)  pobA ← 4‑hydroxybenzoate 3‑monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,010,3750CG 6.8% 192.4 / 2.9 59R82P (CGC→CCC) pobA4‑hydroxybenzoate 3‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base G (2/2);  total (29/30)
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.78e-01

GCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCGGATGCGGCCTAGTACGTACTCGGCA  >  NC_002947/4010235‑4010515
                                                                                                                                             |                                                                                                                                            
gCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCa                                                                                                                                      <  3:89861/149‑1 (MQ=255)
    aaCGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAAc                                                                                                                                    <  2:55823/147‑1 (MQ=255)
     aCGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGGCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACcag                                                                                                                                 <  4:222711/149‑1 (MQ=255)
          ggcggcTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGt                                                                                                                             <  7:197299/148‑1 (MQ=255)
              gcTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAAc                                                                                                                        <  3:213469/149‑1 (MQ=255)
                tGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAAccc                                                                                                                      >  5:38876/1‑149 (MQ=255)
                tGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAAccc                                                                                                                      >  3:235416/1‑149 (MQ=255)
                      aaTGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGGCGCCGCCagcgtc                                                                                                                                 >  1:60785/1‑127 (MQ=255)
                      aaTGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGGCGCCGCCagcgtc                                                                                                                                 >  1:65689/1‑127 (MQ=255)
                      aaTGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGGCGCCGCCagcgtc                                                                                                                                 <  2:65689/132‑6 (MQ=255)
                      aaTGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGGCGCCGCCagcgtc                                                                                                                                 <  2:60785/132‑6 (MQ=255)
                      aaTGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGt                                                                                                                <  2:268741/149‑1 (MQ=255)
                       aTGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGTGGACAGGCGGCGCCCGCCAACCAGCAGTTAAACCCCCTCGt                                                                                                                >  1:110686/1‑149 (MQ=255)
                             ggCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAg                                                                                                         <  4:235416/149‑1 (MQ=255)
                               cgcACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGc                                                                                                       >  5:79659/1‑149 (MQ=255)
                               cgcACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGGCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGc                                                                                                       <  2:125455/149‑1 (MQ=255)
                                         gCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTc                                                                                             >  5:37726/1‑149 (MQ=255)
                                          cTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTcc                                                                                            <  6:38876/149‑1 (MQ=255)
                                                  gCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCa                                                                                                  <  5:88417/135‑1 (MQ=255)
                                                  gCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCa                                                                                                  >  6:88417/1‑135 (MQ=255)
                                                  gCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTc                                                                                    <  6:79659/149‑1 (MQ=255)
                                                  gCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTc                                                                                    >  7:182490/1‑149 (MQ=255)
                                                    cTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTc                                                                                             >  3:100571/1‑138 (MQ=255)
                                                    cTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTc                                                                                             <  4:100571/138‑1 (MQ=255)
                                                     tGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGc                                                                                 >  2:326568/1‑149 (MQ=255)
                                                        aTCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCAc                                                                              >  2:329350/1‑149 (MQ=255)
                                                            ggTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGc                                                                            >  6:84263/1‑147 (MQ=255)
                                                            ggTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCa                                                                          <  8:182490/149‑1 (MQ=255)
                                                                aCGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACcagca                                                                                                                               <  1:147762/92‑1 (MQ=255)
                                                                aCGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACcagca                                                                                                                               >  2:147762/1‑92 (MQ=255)
                                                                   ggTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTcgc                                                                   <  1:329350/149‑1 (MQ=255)
                                                                        cTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCgcagc                                                              >  5:220511/1‑149 (MQ=255)
                                                                        cTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCgcag                                                               >  1:260826/1‑148 (MQ=255)
                                                                            ggTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGt                                                          >  3:269719/1‑149 (MQ=255)
                                                                              tCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCg                                                        >  8:136750/1‑149 (MQ=255)
                                                                                     gTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGc                                                 <  5:84263/149‑1 (MQ=255)
                                                                                           cATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGAc                                                      >  5:178028/1‑138 (MQ=255)
                                                                                           cATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGAc                                                      <  6:178028/138‑1 (MQ=255)
                                                                                           cATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGct                                           >  3:92236/1‑149 (MQ=255)
                                                                                           cATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGc                                            >  3:24429/1‑148 (MQ=255)
                                                                                                          gccgGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACGAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGgcgc                            <  4:24429/149‑1 (MQ=255)
                                                                                                            cgGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCgg                          <  6:37726/149‑1 (MQ=255)
                                                                                                              gTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCGGAt                        <  1:262671/149‑1 (MQ=255)
                                                                                                                 aGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTGCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCGGATGCg                     <  2:110686/149‑1 (MQ=255)
                                                                                                                  gggCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTcgcg                                                                  >  7:171416/1‑103 (MQ=255)
                                                                                                                  gggCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTcgcg                                                                  <  8:171416/103‑1 (MQ=255)
                                                                                                                       ttGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCGGATGCGGCCTAg               >  8:103228/1‑149 (MQ=255)
                                                                                                                         gagaTCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCACCCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCGGATGCGGCCTAgta             >  5:61790/1‑149 (MQ=255)
                                                                                                                         gagaTCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCGGATGCGGCCTAgt              <  7:136750/148‑1 (MQ=255)
                                                                                                                             tCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCGGATGCGGCCTAgtacgta         >  3:233013/1‑149 (MQ=255)
                                                                                                                              ccAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGc                                                                                                              <  3:261931/47‑1 (MQ=255)
                                                                                                                              ccAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGc                                                                                                              >  4:261931/1‑47 (MQ=255)
                                                                                                                                aGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCGGATGCGGCCTAGTACGTACTc      <  5:72793/149‑1 (MQ=255)
                                                                                                                                    gCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCGGATGCg                     >  6:59568/1‑130 (MQ=255)
                                                                                                                                    gCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCGGATGCg                     <  5:59568/130‑1 (MQ=255)
                                                                                                                                    gCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCGGATGCGGCCTAGTACGTACTCgg    >  1:15270/1‑147 (MQ=255)
                                                                                                                                    gCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCGGATGCGGCCTAGTACGTACTCGGCa  <  2:260826/149‑1 (MQ=255)
                                                                                                                                       tGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCg                                                                                   >  7:307899/1‑66 (MQ=255)
                                                                                                                                       tGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCg                                                                                   <  8:307899/66‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                            
GCTGAACGTTGGCGGCTGAATAAATGATCGGCGCACCACTGGCTTCGCGGGCCTGCATCAGGTCACGGGTGACTTCGGTCTGGCCGTAGACCATCACCGTCTTGCCGCCGGTCAGGGCTTTGAGATCCAGACGCTGGCGGCGCCCGCCAACCAGCAGTTCAACCCCCTCGTGCACCAGGCCTTCACGGTCCATGCGCTCGGCCACGCCAGCCTCGCGCAGCAGGTCGACCGTGCCTTGCTCTAGCACCCCGGCGCGGATGCGGCCTAGTACGTACTCGGCA  >  NC_002947/4010235‑4010515

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