Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,420,184 T→A 8.9% M1M (ATG→TTG) †‡ PP_3033 ← CI/C2 family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,420,1840TA 8.9% 162.3 / 4.6 56M1M (ATG→TTG) †‡PP_3033CI/C2 family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (24/27);  new base A (3/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GGCTTGAGTCGTTTTATGAACCCCGACCATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGAAGCTGTTGGAGGTATCGGTGCT  >  NC_002947/3420041‑3420331
                                                                                                                                               |                                                                                                                                                   
ggCTTGAGTCGTTTTATGAACCCCGACCATCTTCAGGCGTTTGTGGACTTCCACGCCTCTGCCTCGGCTTCTAACACCTGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGt                                                                                                                                                 <  7:259628/148‑1 (MQ=255)
   ttGAGTCGTTTTATGAACCCCGACCATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGAt                                                                                                                                             >  7:257606/1‑149 (MQ=255)
       gTCGTTTTATGAACCCCGACCATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGt                                                                                                                                          <  3:301149/148‑1 (MQ=255)
           ttttATGAACCCCGACCATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGca                                                                                                                                      >  5:11682/1‑148 (MQ=255)
               aTGAACCCCGACCATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGtt                                                                                                                                         <  6:378467/141‑1 (MQ=255)
               aTGAACCCCGACCATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGtt                                                                                                                                         >  5:378467/1‑141 (MQ=255)
                tGAACCCCGACCATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACagag                                                                                                                                 <  1:90580/148‑1 (MQ=255)
                 gAACCCCGACCATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGAt                                                                                                                                             <  8:7430/135‑1 (MQ=255)
                 gAACCCCGACCATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGAt                                                                                                                                             >  7:7430/1‑135 (MQ=255)
                    ccccGACCATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGACGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTgc                                                                                                                            >  7:253952/1‑149 (MQ=255)
                     cccGACCATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTgcg                                                                                                                           >  4:124389/1‑149 (MQ=255)
                        gACCATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCaa                                                                                                                        >  2:277532/1‑149 (MQ=255)
                          ccATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATa                                                                                                                      <  2:89038/149‑1 (MQ=255)
                          ccATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATa                                                                                                                      <  1:277532/149‑1 (MQ=255)
                          ccATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATa                                                                                                                      >  8:267653/1‑149 (MQ=255)
                          ccATCTTCAGGCGTGTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATa                                                                                                                      <  7:234741/149‑1 (MQ=255)
                           cATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAg                                                                                                                     >  7:235307/1‑149 (MQ=255)
                                 cAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAg                                                                                                                 <  4:301469/147‑1 (MQ=255)
                                  aGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGtt                                                                                                               <  2:274164/148‑1 (MQ=255)
                                     cGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATc                                                                                                           <  6:108994/149‑1 (MQ=255)
                                       tttGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGCCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCta                                                                                                         <  8:253952/149‑1 (MQ=255)
                                           tgGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTataat                                                                                                     <  3:124389/149‑1 (MQ=255)
                                           tgGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTACACAGAGCTTGCGCAATAGTCAGTTGATCTataa                                                                                                      >  3:128462/1‑148 (MQ=255)
                                              aCGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATc                                                                                                  <  6:11682/149‑1 (MQ=255)
                                                   cACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGAACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTc                                                                                             >  6:135734/1‑149 (MQ=255)
                                                    aCGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTACACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTcc                                                                                            <  4:128462/149‑1 (MQ=255)
                                                        ctctGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGt                                                                                        <  7:267653/149‑1 (MQ=255)
                                                             gcTCGGCTTCTAACACCGGCATCGGCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGAt                                                                                    <  6:378917/147‑1 (MQ=255)
                                                               tCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATa                                                                                   >  3:334362/1‑147 (MQ=255)
                                                                 ggCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAg                                                                                <  8:257606/148‑1 (MQ=255)
                                                                  gCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAg                                                                                >  5:367403/1‑147 (MQ=255)
                                                                           acCGGCATCGTCGAGGGCCTCAAGTAGTCGTGCGCTGAAGGCTGCGCGTAGATCATTTTTATCAAACATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAAATCCGGCTTGATAAAGGTGATTGAAAt                                                                     >  8:188167/1‑149 (MQ=255)
                                                                                      cgAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCg                                                          <  6:134910/149‑1 (MQ=255)
                                                                                               cATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGc                                                             <  4:38202/137‑1 (MQ=255)
                                                                                               cATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGc                                                             >  3:38202/1‑137 (MQ=255)
                                                                                                   tAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGaa                                             >  5:58701/1‑149 (MQ=255)
                                                                                                           cgcTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATTGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATa                                      <  5:135734/148‑1 (MQ=255)
                                                                                                           cgcTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTcc                                                >  1:250030/1‑138 (MQ=255)
                                                                                                           cgcTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTcc                                                <  2:250030/138‑1 (MQ=255)
                                                                                                                       cgcgTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCtttt                          <  4:334362/148‑1 (MQ=255)
                                                                                                                        gcgTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATa                                      >  7:337502/1‑135 (MQ=255)
                                                                                                                        gcgTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATa                                      <  8:337502/135‑1 (MQ=255)
                                                                                                                               tCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGt                                                                                                                    >  8:251338/1‑50 (MQ=255)
                                                                                                                               tCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGAAGCTGtt                 >  4:69209/1‑149 (MQ=255)
                                                                                                                                 aTTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGAAGCTGTTg                <  8:235307/148‑1 (MQ=255)
                                                                                                                                   ttttATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGAAGCTGTTGGAg             >  7:308388/1‑149 (MQ=255)
                                                                                                                                   aaccaTCAACCCATGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGAAGCTGTTGGAg             >  3:27586/5‑149 (MQ=255)
                                                                                                                                    tttATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGAAGCTGTTGGAg             <  6:367403/148‑1 (MQ=255)
                                                                                                                                    tttATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGAAGCTGTTGGAg             <  6:58701/148‑1 (MQ=255)
                                                                                                                                    accaTCAACCCATGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGaa                       <  1:116209/135‑1 (MQ=255)
                                                                                                                                    accaTCAACCCATGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGaa                       >  2:116209/4‑138 (MQ=255)
                                                                                                                                     ccaTCAACCCATGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGAAGCTGTTGGa              <  2:87972/144‑1 (MQ=255)
                                                                                                                                     ccaTCAACCCATGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGAAGCTGTTGGa              >  1:87972/3‑146 (MQ=255)
                                                                                                                                       aTCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGAAGCTGTTGGAGGTAt         >  6:42094/1‑149 (MQ=255)
                                                                                                                                          aaCCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGAAGCTGTTGGAGGTATCgg      >  7:166403/1‑149 (MQ=255)
                                                                                                                                              aTGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGAAGCTGTTGGAGGTATCGGtgct  >  3:101954/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGCTTGAGTCGTTTTATGAACCCCGACCATCTTCAGGCGTTTGTGGACGTCCACGCCTCTGCCTCGGCTTCTAACACCGGCATCGTCGAGGGCCTCATGTAGTCGTGCGCTGAAGGCTGCGCGTAGCTCATTTTTATCAACCATGGGTTGATGGTTGCACAGAGCTTGCGCAATAGTCAGTTGATCTATAATATCAACTCCGGGTTGATAAAGGAGATTGAAATGCTGCGGCCCGAAGAATTTCCGAATGCGATAGCGTTCGCTTTTGAAGCTGTTGGAGGTATCGGTGCT  >  NC_002947/3420041‑3420331

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