Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,845,703:1 +G 100% intergenic (‑70/+28) PP_3394 ← / ← PP_3395 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,845,7031.G94.3% 168.6 / 3.2 53intergenic (‑70/+28)PP_3394/PP_33953‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base G (32/18);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 5.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGC  >  NC_002947/3845564‑3845842
                                                                                                                                               |                                                                                                                                           
tCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATc                                                                                                                                     >  6:277481/1‑149 (MQ=255)
    catcatCAGTTGCTGGCTCTGCAAGCCTTCGCGCAGGCCGACTTCATTGATTGTTTTTGTTGGCATCGCGGTCTCTCTGTCTGAATCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATc                                                                                                                                 <  4:7722/149‑1 (MQ=255)
    cATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATc                                                                                                                                 <  4:103630/149‑1 (MQ=255)
        aTCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCa                                                                                                                             <  5:277481/149‑1 (MQ=255)
           aCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGc                                                                                                                           >  7:75544/1‑148 (MQ=255)
             ttGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCg                                                                                                                                <  5:116447/141‑1 (MQ=255)
             ttGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCg                                                                                                                                >  6:116447/1‑141 (MQ=255)
               gCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGa                                                                                                                      >  7:140160/1‑149 (MQ=255)
                   gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGAGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGa                                                                                                                  >  1:14272/1‑149 (MQ=255)
                    cttgcagaGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCCGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGAt                                                                                                                 <  2:170049/142‑1 (MQ=255)
                    ctctGCTAGGCATCGCGCTGGACGAGTTGATTGATTGTTATTGTTGTCATCGCGGTCTCTGTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGGAGTAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCAggatggat                                                                                                                 <  4:184500/149‑1 (MQ=255)
                     tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGAt                                                                                                                     <  8:188374/144‑1 (MQ=255)
                     tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGAt                                                                                                                     >  7:188374/1‑144 (MQ=255)
                       gcagaGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCa                                                                                                               <  6:92604/144‑1 (MQ=255)
                                  gcgcAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGATGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTg                                                                                                   >  1:192728/1‑149 (MQ=255)
                                   cgcAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGc                                                                                                  >  1:337526/1‑149 (MQ=255)
                                       ggCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGAACCATAGTAAGCCACGATGGATCATCCCGCGCCCTGCCAgg                                                                                              >  1:131034/1‑149 (MQ=255)
                                          cGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGAt                                                                                                                 >  6:135083/1‑127 (MQ=255)
                                          cGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGAt                                                                                                                 <  5:135083/127‑1 (MQ=255)
                                          cGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCg                                                                                           >  6:170779/1‑149 (MQ=255)
                                               tcttttattttttttgttGTCATCGCGGTCTCTCTGTCTTAAGCTTTGCGCCAGAATAAGTGCTGTTTCAGGGGAATAAAACACGGTTTTCGGGGGCTTATCCATCTTCAGCCACTATGTATCAGCCCGCGCCCTGCCAGGCCGGTGc                                                                                       <  8:75544/142‑1 (MQ=255)
                                                 attgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACg                                                                                    <  5:274155/149‑1 (MQ=255)
                                                  ttgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCAca                                                                                    >  3:198727/1‑147 (MQ=255)
                                                  ttgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACg                                                                                    >  3:204645/1‑148 (MQ=255)
                                                     attgttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCCCgcag                                                                                 >  2:209272/1‑148 (MQ=255)
                                                     attgttattgttGTCATCGCGGTCTCTCTCTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCCGCCCGCGCCCTGCCAGGCCGGTGCACgcag                                                                                 >  6:284247/1‑148 (MQ=255)
                                                        gttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCAAACTGCGCCACGATGGCTCAGCCCGCGCCCAGCCAGGCCGGTGGACGCAGCaa                                                                              >  6:44534/1‑148 (MQ=255)
                                                              gttttCATGGGGGTCTCTGTGTCTTAAGGTTGCCCCCATAATAAGAGCTGTATCCGGGCAACAAAACCCGTGTTGCGCGGGCTCACCCATCGCCCCCCACTTGGTATCCGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtg                                                                        <  8:140160/148‑1 (MQ=255)
                                                                             tctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccga                                                        >  3:115434/1‑149 (MQ=255)
                                                                             tctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgcc                                                          >  5:301251/1‑147 (MQ=255)
                                                                               tgtgtctgAAGCGTTGCGCCAGACTATGAGCTGGATGACGGGCAGAAACGCCGTGTTTCGGGGGCTTATCCTTCGTCAGCCACGATGGATCTGCCCGCGCCCTGCCAGGCCGGTGc                                                                                       <  7:8731/114‑1 (MQ=255)
                                                                               tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGc                                                                                       >  8:8731/1‑116 (MQ=255)
                                                                               tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                                      >  1:103443/1‑149 (MQ=255)
                                                                               tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATACATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGc                                                                 >  1:163777/1‑138 (MQ=255)
                                                                                      gAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCAGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                               >  1:47454/1‑149 (MQ=255)
                                                                                       aaGCGTTGCGCCTTATTAATAGCTGTTTCTGGGGCAAAAAAGACGTGTTTCGGGGGCTTATCCATCGTCAGCCAGGATGTATAAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACa                                                 <  2:192728/146‑1 (MQ=255)
                                                                                        aGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCga                                             <  5:170779/149‑1 (MQ=255)
                                                                                                 gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtgt                                     >  2:112469/1‑149 (MQ=255)
                                                                                                   caaaCTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGc                                   >  3:16156/4‑149 (MQ=255)
                                                                                                   aaGACTAAGAGCTGGATGACGGGAAGAAAACACGGGTTTCGCGGGCGGATCCATCGTCAGCCACGATGTTTCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGc                                   <  4:204645/148‑1 (MQ=255)
                                                                                                     gACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCg                                 >  5:155126/1‑149 (MQ=255)
                                                                                                             gCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATACATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCa                         >  7:301794/1‑149 (MQ=255)
                                                                                                               tGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCAGATACATAGGGAGCCACGAAGGATAAGCCCCCCCCCAGGCAGGCCCGGGGAAGCCGCAAAAGGGGGTGAAACAGGCGCCGAGCCGCCACCGAGAGGGGgtggcggtggcgaga                        >  2:310563/1‑145 (MQ=255)
                                                                                                                 gttGACGGGCAGAAAAGACGTGTGTCGCGGGCTGAGCCATCGTCATCCACGATGGAGCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCa                     <  1:209272/147‑1 (MQ=255)
                                                                                                                    gACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGa                   <  4:115434/148‑1 (MQ=255)
                                                                                                                         gCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCACGCGCCCTGGCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtggg                                    >  7:160357/1‑124 (MQ=255)
                                                                                                                            aGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAt           >  7:39558/1‑149 (MQ=255)
                                                                                                                             gAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAtt          <  6:155126/149‑1 (MQ=255)
                                                                                                                              aaaaGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAtt          >  6:299530/1‑148 (MQ=255)
                                                                                                                                 aGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGcc                                  >  8:1423/1‑121 (MQ=255)
                                                                                                                                     gtgtTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGa             <  1:214750/138‑1 (MQ=255)
                                                                                                                                     gtgtTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGa             >  2:214750/1‑138 (MQ=255)
                                                                                                                                     gtgtTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgc  <  2:103443/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                           
TCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGC  >  NC_002947/3845564‑3845842

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: