Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,166,086 C→G 9.6% F358L (TTC→TTG baeS → sensory histidine kinase in two‑component regulatory system with BaeR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,166,0860CG 9.6% 204.6 / 9.7 73F358L (TTC→TTGbaeSsensory histidine kinase in two‑component regulatory system with BaeR
Reads supporting (aligned to +/- strand):  ref base C (33/33);  new base G (4/3);  total (37/36)
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.30e-01

GAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGCCTGCAAATCTCTGCCGGGCAGC  >  AP009048/2166017‑2166156
                                                                     |                                                                      
gAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTc                                                                        <  1:935868/70‑1 (MQ=255)
 aaaCTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTc                                                                        <  1:2483291/69‑1 (MQ=255)
  aaCTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCaa                                                                      >  1:4618909/1‑70 (MQ=255)
  aaCTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCaa                                                                      >  1:2192413/1‑70 (MQ=255)
  aaCTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAAt                                                                     <  1:3347698/71‑1 (MQ=255)
       cAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACt                                                                 >  1:6239856/1‑70 (MQ=255)
       cAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACt                                                                 >  1:5360162/1‑70 (MQ=255)
          ttttCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACt                                                             >  1:3967901/1‑71 (MQ=255)
          ttttCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACt                                                             >  1:5318924/1‑71 (MQ=255)
            ttCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTg                                                            >  1:4730050/1‑70 (MQ=255)
            ttCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTg                                                            >  1:7030763/1‑70 (MQ=255)
            ttCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTg                                                            >  1:5121087/1‑70 (MQ=255)
             tCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTgg                                                           <  1:3668325/70‑1 (MQ=255)
             tCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTgg                                                           <  1:5399495/70‑1 (MQ=255)
               ccTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTGAATAACt                                                                 <  1:1014351/62‑1 (MQ=255)
               ccTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTGAATAACt                                                                 <  1:4955852/62‑1 (MQ=255)
               ccTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTGAATAACt                                                                 <  1:1132138/62‑1 (MQ=255)
                 tGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTGAATAACtgc                                                               >  1:7824385/1‑60 (MQ=255)
                 tGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTGAATAACtgc                                                               >  1:1522469/1‑60 (MQ=255)
                 tGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTGAATAACtgc                                                               >  1:1166638/1‑60 (MQ=255)
                 tGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaaaa                                                       >  1:4771210/1‑70 (MQ=255)
                 tGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaa                                                         <  1:7904657/68‑1 (MQ=255)
                  gCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAAc                                                                  >  1:4057548/1‑58 (MQ=255)
                   ccAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAg                                                    >  1:1695420/1‑71 (MQ=255)
                    cAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGc                                                   <  1:7597285/71‑1 (MQ=255)
                      gACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCt                                                 >  1:1091781/1‑71 (MQ=255)
                        cAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTc                                                                        >  1:35355/1‑46 (MQ=255)
                           tATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaa                                                         <  1:7170704/58‑1 (MQ=255)
                             ttACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTa                                           <  1:5664580/70‑1 (MQ=255)
                             ttACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTa                                           >  1:3376327/1‑70 (MQ=255)
                              tACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaa                                                         <  1:2755632/55‑1 (MQ=255)
                              tACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaa                                                         <  1:6097063/55‑1 (MQ=255)
                                  gTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACt                                       >  1:6815972/1‑69 (MQ=255)
                                     tttGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACt                                                                 <  1:6726948/40‑1 (MQ=255)
                                     tttGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaaaa                                                       <  1:4218827/50‑1 (MQ=255)
                                     tttGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaaa                                                        >  1:6126082/1‑49 (MQ=255)
                                     tttGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACa                                   <  1:628714/70‑1 (MQ=255)
                                      ttGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGATACACTGACAgc                                 >  1:5134065/1‑71 (MQ=255)
                                       tGGCGATCGCGACCGTTTAATGCAGTTATTGAATAACTTa                                                               >  1:7412472/1‑40 (MQ=255)
                                       tGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTa                                                               >  1:2446147/1‑40 (MQ=255)
                                       tGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAgc                                 <  1:5702678/70‑1 (MQ=255)
                                        ggCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTTGAAAACAGCCTGCGCTACACTGACAgcg                                <  1:874031/70‑1 (MQ=255)
                                         gCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAgcgg                               <  1:2335044/70‑1 (MQ=255)
                                          cGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaaaa                                                       >  1:7420997/1‑45 (MQ=255)
                                          cGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaaaa                                                       >  1:35130/1‑45 (MQ=255)
                                          cGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaaaa                                                       >  1:2528989/1‑45 (MQ=255)
                                          cGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaaaa                                                       >  1:3806206/1‑45 (MQ=255)
                                          cGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAgcggc                              >  1:5217066/1‑70 (MQ=255)
                                           gATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaaaa                                                       <  1:1175674/44‑1 (MQ=255)
                                           gATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaaaa                                                       <  1:3757411/44‑1 (MQ=255)
                                           gATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaa                                                         <  1:6864803/42‑1 (MQ=255)
                                           gATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaa                                                         <  1:5054033/42‑1 (MQ=255)
                                           gATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGaa                                                         <  1:2695043/42‑1 (MQ=255)
                                           gATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAg                                  >  1:3855255/1‑65 (MQ=255)
                                           gATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGcggcgg                            <  1:6847980/71‑1 (MQ=255)
                                              cgcgACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACa                                                     <  1:5824865/43‑1 (MQ=255)
                                              cgcgACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCa                          <  1:4593389/70‑1 (MQ=255)
                                              cgcgACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAg                         >  1:7338904/1‑71 (MQ=255)
                                                cgACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGcc                       <  1:323712/71‑1 (MQ=255)
                                                 gACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGcc                                                  <  1:3073831/43‑1 (MQ=255)
                                                 gACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGCCt                      >  1:1963251/1‑71 (MQ=255)
                                                   ccGTTTAATGCAGTTATTCAATAACTTACTGGaaaa                                                       >  1:74614/1‑36 (MQ=255)
                                                      tttAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGcc                       <  1:1825160/65‑1 (MQ=255)
                                                      tttAATGCAGTTATTCAATAACTTACTGGAAAACAGCCCGCGCTACACTGACAGCGGCGGCAGcc                       <  1:3264339/65‑1 (MQ=255)
                                                       ttAATGCAGTTATTCAATAACTTACTGGAAAACAGc                                                   >  1:2035617/1‑36 (MQ=255)
                                                       ttAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTaca                                         >  1:1414423/1‑46 (MQ=255)
                                                       ttAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTaca                                         >  1:961388/1‑46 (MQ=255)
                                                       ttAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGCCTGCAAAt                >  1:3589579/1‑71 (MQ=255)
                                                        tAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGCCTg                     >  1:3236194/1‑65 (MQ=255)
                                                         aaTGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGCCTGCAAAtc               <  1:2508034/70‑1 (MQ=255)
                                                           tGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGCCTGCAAATctct            <  1:6532683/71‑1 (MQ=255)
                                                                ttATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAg                         >  1:4698342/1‑53 (MQ=255)
                                                                  aTTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGCCTGCAAATCTCTGCCggg      <  1:2765074/70‑1 (MQ=255)
                                                                   ttCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGCCTGCaa                  <  1:2262789/57‑1 (MQ=255)
                                                                   ttCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGCCTGCAAATCTCTGCCGGGc     <  1:2552376/70‑1 (MQ=255)
                                                                   ttCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGCCTGCAAATCTCTGCCGGGc     <  1:7357825/70‑1 (MQ=255)
                                                                     cAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGCCTGCAAATCTCTGCCGGGCAgc  >  1:348559/1‑71 (MQ=255)
                                                                     cAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGCCTGCAAATCTCTGCCGGGCAg   >  1:3073350/1‑70 (MQ=255)
                                                                     |                                                                      
GAAACTGCAATTTTCCCTGCCAGACAGTATTACCGTATTTGGCGATCGCGACCGTTTAATGCAGTTATTCAATAACTTACTGGAAAACAGCCTGCGCTACACTGACAGCGGCGGCAGCCTGCAAATCTCTGCCGGGCAGC  >  AP009048/2166017‑2166156

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