Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,048,354 C→G 13.0% T258R (ACG→AGG)  PP_5300 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,048,3540CG13.0% 149.9 / 13.9 54T258R (ACG→AGG) PP_5300LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (22/25);  new base G (4/3);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

ACACCTCCGGAGCACCTGGTTGGCCGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATGCGCTGCCCGGGGTGGATGCAC  >  NC_002947/6048209‑6048503
                                                                                                                                                 |                                                                                                                                                     
acacCTCCGGAGCACCTGGTTGGCCGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGGGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACggtg                                                                                                                                                   >  2:809475/1‑150 (MQ=255)
acacCTCCGGAGCACCTGGTTGGCCGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACggtg                                                                                                                                                   >  1:1105238/1‑150 (MQ=255)
acacCTCCGGAGCACCTGGTTGGCCGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACggtg                                                                                                                                                   <  2:754844/150‑1 (MQ=255)
   ccgcctgAGCACCGGTGTGGCCGCAACCGGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTc                                                                                                                                                >  2:16464/7‑150 (MQ=255)
       cGGAGCACCTGGTTGGCCGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGgc                                                                                                                                            <  1:94494/150‑1 (MQ=255)
       cGGAGCACCTGGTTGGCCGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGgc                                                                                                                                            <  2:1022962/150‑1 (MQ=255)
         gAGCACCTGGTTGGCCGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCGTAGGGTGGTcctgcat                                                                                                                                          <  1:538131/150‑7 (MQ=255)
            cACCTGGTTGGCCGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCTGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGcaac                                                                                                                                       >  2:182981/1‑150 (MQ=255)
            cACCTGGTTGGCCGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGcaac                                                                                                                                       >  1:983110/1‑150 (MQ=255)
             aCCTGGTTGGCCGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGcaacc                                                                                                                                      <  1:705552/150‑1 (MQ=255)
               cTGGTTGGCCGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGcaaccaa                                                                                                                                    >  2:340517/1‑150 (MQ=255)
                tGGTTGGCCGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCTGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGcaaccaac                                                                                                                                   <  1:182981/150‑1 (MQ=255)
                       ccGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACAccc                                                                                                                            <  1:340517/150‑1 (MQ=255)
                       ccGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGCCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGCATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCACCAACcccc                                                                                                                            >  1:610479/1‑150 (MQ=255)
                       ccGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCACCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGAACCACCCTACGGTGGTCTGGCGCAACCAAAAacactc                                                                                                                            >  2:183047/1‑148 (MQ=255)
                             ccTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCt                                                                                                                      <  1:464831/150‑1 (MQ=255)
                                          tcctATGTAGTCTGCGGCCAGCCACAGTGGCGCTAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTAc                                                                                                         <  2:697558/150‑1 (MQ=255)
                                          tcctATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTAc                                                                                                         <  2:1105238/150‑1 (MQ=255)
                                          tcctATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTAc                                                                                                         <  2:98107/150‑1 (MQ=255)
                                           cctATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGcaaccaa                                                                                                                                    <  2:364058/122‑1 (MQ=255)
                                           cctATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGcaaccaa                                                                                                                                    >  1:364058/1‑122 (MQ=255)
                                            ctATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCa                                                                                                       >  1:978039/1‑150 (MQ=255)
                                            ctATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCa                                                                                                       <  2:983110/150‑1 (MQ=255)
                                                 tAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTgca                                                                                                  <  2:978039/150‑1 (MQ=255)
                                                                cacaGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTAc                                                                                                                                                       <  1:568802/82‑1 (MQ=255)
                                                                cacaGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTAc                                                                                                                                                       >  2:568802/1‑82 (MQ=255)
                                                                   aggggCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCg                                                                                <  1:16464/147‑1 (MQ=255)
                                                                           gaccgcCCCAGCTAGCGCCCCGCCAGCGTGCCGTGGATCGCACCAGACGGCTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTg                                                                        <  2:610479/147‑1 (MQ=255)
                                                                           gAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTg                                                                        <  1:809475/150‑1 (MQ=255)
                                                                                ccTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAAt                                                                   >  1:1192845/1‑150 (MQ=255)
                                                                                 cTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATg                                                                  >  1:794420/1‑150 (MQ=255)
                                                                                  tCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTgca                                                                                                  <  1:687747/117‑1 (MQ=255)
                                                                                  tCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTgca                                                                                                  >  2:687747/1‑117 (MQ=255)
                                                                                    aGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAAAACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGTATGTCGACCATCGCCCAGTTGGTAATggcc                                                               <  1:183047/150‑1 (MQ=255)
                                                                                     gCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGt                                                                                         <  1:783785/123‑1 (MQ=255)
                                                                                     gCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGt                                                                                         >  2:783785/1‑123 (MQ=255)
                                                                                             cccGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGc                                                      >  2:391325/1‑150 (MQ=255)
                                                                                               cGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGt                                                    <  2:652915/150‑1 (MQ=255)
                                                                                                                        gacgacTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTg                           <  1:432637/150‑1 (MQ=255)
                                                                                                                              tCGCTGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTAGCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATgcg                     >  2:866678/1‑150 (MQ=255)
                                                                                                                                   gccagACCACCGTAGGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATgcgc                    >  1:305217/5‑146 (MQ=255)
                                                                                                                                   gccagACCACCGTAGGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATgcgc                    <  2:305217/142‑1 (MQ=255)
                                                                                                                                   gccagACCACCGTAGGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATGCGCTGcc                >  1:619868/5‑150 (MQ=255)
                                                                                                                                    ccagACCACCGTAGGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATGCGCTGccc               >  2:201273/4‑150 (MQ=255)
                                                                                                                                    ccagACCACCGTAGGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATGCGCTGccc               >  1:1186412/4‑150 (MQ=255)
                                                                                                                                    ccagACCACCCTACGTTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTgg                                                  >  1:1131677/4‑115 (MQ=255)
                                                                                                                                    ccagACCACCCTACGTTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTgg                                                  <  2:1131677/112‑1 (MQ=255)
                                                                                                                                       gACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCCTGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATGCGCTGCCCggg            >  2:41133/1‑150 (MQ=255)
                                                                                                                                         ccaccGTAGGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATGCGCTGCCCGGGGt          <  2:619868/150‑1 (MQ=255)
                                                                                                                                         ccaccCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATg                       <  2:1175205/137‑1 (MQ=255)
                                                                                                                                         ccaccCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATg                       >  1:1175205/1‑137 (MQ=255)
                                                                                                                                         ccaccCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATGCGCTGCCCGGGGt          <  2:1192845/150‑1 (MQ=255)
                                                                                                                                          caccCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATGCGCTGCCCGGGGTg         >  2:584813/1‑150 (MQ=255)
                                                                                                                                          caccCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATGCGCTGCCCGGGGTg         >  2:309673/1‑150 (MQ=255)
                                                                                                                                            cccTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATGCGCTGCCCGGGGTGGa       >  1:867830/1‑150 (MQ=255)
                                                                                                                                                aCGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATGCGCTGCCCGGGGTGGGTGCa   <  2:794420/150‑1 (MQ=255)
                                                                                                                                                 cGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATGCGCTGCCCGGGGTGGATGCAc  <  1:619781/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
ACACCTCCGGAGCACCTGGTTGGCCGCAACCTGGGTTCGACCTCCTATGTAGTCTGCGGCCAGCCACAGTGGCGCGAGCGCCTCAGCGAGTCCCCCGCCAGCGTGCCGTGGATCGCACCAGACGACTCGATGCAGGACCACCCTACGGTGGTCTGGCGCAACCAACAACACCCGGGCCTGAGCCCGCGCTACCAGTGCAGCAGCATGTCGACCATCGCCCAGTTGGTAATGGCCGGCCTGGGCGTGGCAGCACTGCCGGACTACATGGTGCATGCGCTGCCCGGGGTGGATGCAC  >  NC_002947/6048209‑6048503

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