Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,370,345 T→C 7.5% A13A (GCT→GCC PP_3847 → LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,370,3450TC 7.5% 146.5 / 2.8 53A13A (GCT→GCCPP_3847LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (25/24);  new base C (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

CACGTGTATGTTGGCTCAATCGCCAAACTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGGGTCGTGCTGATATTCAGCAGTGGTCGCCCGCAGGTGCTCGCTGAAT  >  NC_002947/4370199‑4370493
                                                                                                                                                       |                                                                                                                                                    
cACGTGTATGTTGGCTCAATCGCCAAACTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCtt                                                                                                                                                     <  7:29496/148‑1 (MQ=255)
   gtgtATGTTGGCTCAATCGCCAAACTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCtt                                                                                                                                                     <  8:197537/145‑1 (MQ=255)
   gtgtATGTTGGCTCAATCGCCAAACTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCtt                                                                                                                                                     >  7:197537/1‑145 (MQ=255)
   gtgtATGTTGGCTCAATCGCCAAACTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCa                                                                                                                                                 >  4:301679/1‑149 (MQ=255)
            ggCTCAATCGCCAAACTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCAt                                                                                                                                        >  6:284826/1‑149 (MQ=255)
                 aaTCGCCAAACTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCa                                                                                                                                   <  3:301679/149‑1 (MQ=255)
                     gCCAAACTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGGAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGAc                                                                                                                                >  2:181362/1‑148 (MQ=255)
                     gCCAAACTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGAc                                                                                                                                >  3:243819/1‑148 (MQ=255)
                      ccAAACTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACt                                                                                                                               >  6:295101/1‑148 (MQ=255)
                        aaaCTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCCTGGTGCCAAGCCAtgtc                                                                                                                                     >  7:47911/1‑137 (MQ=255)
                        aaaCTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCCTGGTGCCAAGCCAtgtc                                                                                                                                     <  8:47911/140‑4 (MQ=255)
                         aaCTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATc                                                                                                                           >  7:139334/1‑149 (MQ=255)
                           cTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGa                                                                                                                         <  5:284826/149‑1 (MQ=255)
                            tGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGACAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGACCAT‑‑TCCACTGCAGGA‑‑‑CATGACTTGGCACCAGGCCATCGGCACACTAATCTAg                                                                                                                        >  6:211936/1‑149 (MQ=255)
                                      ccGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGAcc                                                                                                              >  8:176558/1‑149 (MQ=255)
                                      ccGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATAACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGAcc                                                                                                              >  2:108595/1‑149 (MQ=255)
                                        gCCAGGTCTGAAACTAAGTATGTAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGTAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑GATGACTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGc                                                                                                            <  5:211936/149‑1 (MQ=255)
                                                tGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTgg                                                                                                                                                   <  8:59181/102‑1 (MQ=255)
                                                   aaCTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTct                                                                                                  <  7:130764/148‑1 (MQ=255)
                                                     cTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCCTGGTGCCAAGCCAtgtc                                                                                                                                     >  2:311365/1‑108 (MQ=255)
                                                     cTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCCTGGTGCCAAGCCAtgtc                                                                                                                                     <  1:311365/111‑4 (MQ=255)
                                                        aGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCg                                                                                                                          >  3:216286/1‑119 (MQ=255)
                                                        aGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCg                                                                                                                          <  4:216286/119‑1 (MQ=255)
                                                        aGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTctggctg                                                                                             <  6:181341/148‑1 (MQ=255)
                                                              gagaCCCGTAACACGGAAAGTAATTACTTCGATGGGTTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTgg                                                                                      <  3:110235/149‑1 (MQ=255)
                                                               agagCCGTAAAACGGTAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTgg                                                                                      <  1:181362/148‑1 (MQ=255)
                                                                   ccGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTa                                                                                                                              <  5:131149/104‑1 (MQ=255)
                                                                   ccGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTa                                                                                                                              >  6:131149/1‑104 (MQ=255)
                                                                   ccGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGt                                                                                 >  8:70545/1‑149 (MQ=255)
                                                                      tAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATTGGTAGAAACGGGAGCATTCTACACTGCAGGATCCCCCTGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGc                                                                                   <  1:108595/149‑1 (MQ=255)
                                                                           cGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGga                                                                         <  6:187973/149‑1 (MQ=255)
                                                                            ggAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGgat                                                                        <  8:139334/149‑1 (MQ=255)
                                                                                 gTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGgatagata                                                                   <  5:295101/149‑1 (MQ=255)
                                                                                  tAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGga                                                                         >  4:39298/1‑142 (MQ=255)
                                                                                  tAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGga                                                                         <  3:39298/142‑1 (MQ=255)
                                                                                         ttCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATc                                                                                                                           <  5:317473/85‑1 (MQ=255)
                                                                                         ttCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATc                                                                                                                           >  6:317473/1‑85 (MQ=255)
                                                                                            gATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGgataga                                                                     >  6:10101/1‑136 (MQ=255)
                                                                                            gATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTg                                                         <  7:176558/148‑1 (MQ=255)
                                                                                              tGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTctgg                                                                                                <  8:289260/107‑1 (MQ=255)
                                                                                              tGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTctgg                                                                                                >  7:289260/1‑107 (MQ=255)
                                                                                                       cccGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGTACTTGACCGCCCTAATTTct                                                                                                  >  6:345072/1‑96 (MQ=255)
                                                                                                       cccGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGTACTTGACCGCCCTAATTTct                                                                                                  <  5:345072/96‑1 (MQ=255)
                                                                                                                gTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGGGTCGTGCTGata                                    >  3:195375/1‑149 (MQ=255)
                                                                                                                 tAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGa                                                                                                                                 <  2:315326/55‑1 (MQ=255)
                                                                                                                 tAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGa                                                                                                                                 >  1:315326/1‑55 (MQ=255)
                                                                                                                      aCGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGGGTCGTGCTGATATTcag                               >  1:55800/1‑148 (MQ=255)
                                                                                                                        gggAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGTCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGTGTCGTGCTGATATTcagaag                            >  4:59604/1‑149 (MQ=255)
                                                                                                                           aGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGGGTCGTGCTGATATTcagcag                            <  1:29737/146‑1 (MQ=255)
                                                                                                                           aGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGGGTCGTGCTGATATTcagcag                            >  2:29737/1‑146 (MQ=255)
                                                                                                                               t‑‑tACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGGGTCGTGCTGATATTCAGCAGTGGTCGc                     >  3:141288/1‑149 (MQ=255)
                                                                                                                                                   tGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGGGTCGTGCTGATATTCAGCAGTGGTCGCCCGCAGGTGCTCGCt      <  4:195375/149‑1 (MQ=255)
                                                                                                                                                     gCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGGGTCGTGCTGATATTca                                >  5:192761/1‑121 (MQ=255)
                                                                                                                                                     gCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGGGTCGTGCTGATATTca                                >  3:151912/1‑121 (MQ=255)
                                                                                                                                                     gCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGGGTCGTGCTGATATTca                                <  4:151912/121‑1 (MQ=255)
                                                                                                                                                     gCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGGGTCGTGCTGATATTca                                <  6:192761/121‑1 (MQ=255)
                                                                                                                                                       ttGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGGGTCGTGCTGATATTCAGCAGTGGTCGCCCGCACGTGCTCGCTGAAt  >  8:223354/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
CACGTGTATGTTGGCTCAATCGCCAAACTGAACAGCACCCGCCAGGGCTGAAACTAAGTATGGAGAGCCGTAAAACGGAAAGTAATTACTTCGATGGGGTAACCCCGTATGGGTAGAAACGGGAGCAT‑‑TACACTGCAGGA‑‑‑CATGGCTTGGCACCAGGCCATCGGCAGACTAATCGAGGCACTTGACCGCCCTAATTTCTGGCTGGCCTTGGTGCGTCAACTGGATAGATACATTCCGGTGGACAGTTGGGTCGTGCTGATATTCAGCAGTGGTCGCCCGCAGGTGCTCGCTGAAT  >  NC_002947/4370199‑4370493

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: