Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,379,843 C→G 12.5% A569P (GCG→CCG)  aaeB ← p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,379,8430CG12.5% 89.9 / 5.9 32A569P (GCG→CCG) aaeBp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base C (17/11);  new base G (2/2);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTACCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAGATATCACATGGTCTGCTGTGCGGCGCATTTGTCGGTGAAACGCCGATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTG  >  CP009273/3379590‑3380111
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                            
gTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTACCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAGATATCACATGGTCTGCTGTGCGGCGCATTTGTCGGTGAAACGCCGAAAAATCCTCGTTAAACGGGAACGGTGCATCAAGCAGGCGACGGTGCGCGATAATCATAGTCAGcccccg                                                                                                                                                                                                                                                         >  1:200548/1‑273 (MQ=255)
                                                                                         ttCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAGATATCACATGGTCTGCTGTGCGGCGCATTTGTCGGTGAACCGCCGATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGGGGTTCTCTTACAGACGTg                                                                                                                                                                >  1:82827/1‑275 (MQ=255)
                                                                                             tGGTAGATTTCCAGGTCTGACAGCAACTGGCCGAGGTGCCTACGACGTTTATCTGTGCTACGGGCAGGTATCACATGGGCTGCTGTGCGGCGCTTTTGTCGGTGGAACGCCGATAAATCGTCGTTAACCGGGATGGGTGCATCACGCGGGCCGGGGTGCGCGATAAGCATCGTCAGCGCCAGGCGAAATGTCTGCAAATCCCCTGGGAACTTATGCATCAGCAAAAACAGCCGCTGACAAAGTGCCGGGAGGTGGTTCTCTTAACGACGTg                                                                                                                                                                <  2:88507/271‑1 (MQ=255)
                                                                                                                                          ggTTTGTCATCGCTACTCGACGATAGCCCATGGTCTGCTGTTCGGCGCATTTGTCGGTGAAACGCCGATAAATCCTCGTTTACCGGGATCGGTGCAGCACTCCGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAg                                                                                                               <  2:82825/274‑1 (MQ=255)
                                                                                                                                                                     tgactttccatttcagtcagtttctcgacagcgcattagggcaaatcgtcggctgtgtgctcgcgttcaccgttattttGCTGGTGCGGGATAA‑‑ATCGCGCGACAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCaaaaa                                                                                                                                                                                                            <  1:171149‑M1/72‑1 (MQ=255)
                                                                                                                                                                     tgactttccatttcagtcagtgtgacgacagcgcattagggcaaatcgtcggctgtgtgctcgcgttcaccgttattttGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAggtg                                                                                                                                                                                 <  1:230869‑M1/103‑1 (MQ=255)
                                                                                                                                                                         tttccatttcagtcagtttctcgacagcgcattagggcaaatcgtcggctgtgtgctcgcgttcaccgttattttGCTGGTGCGGGATAA‑‑ATCGCGCGACAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAgc                                                                                                                                                                                                        >  2:171145‑M1/76‑151 (MQ=255)
                                                                                                                                                                         tttccatttcagtcagtttctcgacagcgcattagggcaaatcgtcggctgtgtgctcgcgttcaccgttattttGCTGGTGCGCGa                                                                                                                                                                                                                                                                            >  2:71862‑M1/76‑87 (MQ=255)
                                                                                                                                                                         tttccatttcagtcagtttctcgacagcgcattagggcaaatcgtcggctgtgtgctcgcgttcaccgttattttGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTtc                                                                                                                                                                             >  2:230862‑M1/76‑182 (MQ=255)
                                                                                                                                                                         tttccatttcagtcagtttctcgacagcgcattagggaaaatcgtcggctgtgtgctcgcgttcaccgttattttGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCGGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGCCACACTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTg                                                                                                                     >  1:90516‑M1/76‑238 (MQ=255)
                                                                                                                                                                          gTCTGCTGTGTGGCGCATTTGTCGGTGAAACGCCGATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGCTACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTAt                                                                               >  1:327735/1‑275 (MQ=255)
                                                                                                                                                                          gTCTGCTGTGCGGCGCATTTGTCTGTGAAACGCCGATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCgg                                                                                                                                  >  1:235366/1‑224 (MQ=255)
                                                                                                                                                                          gTCTGCTGTGCGGCGCATTTGTCGGTGAAACGCCGATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGAGCGATAATCATCGTCCGCGCCAGGCGAAATTTCg                                                                                                                                                                                                                                              >  2:272873/1‑116 (MQ=255)
                                                                                                                                                                             tgctgTGCGGCGCATTTGTCGGTGAAACGCCGATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTNCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCAGGTCCTGTCGCGCGATTTATccc                                                                            >  1:288677/1‑275 (MQ=255)
                                                                                                                                                                              gctgTGCGGCGCATTTGTCGGTGAAACGCCGATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGCGCGATAATCATCGTCAGCTCCAGGCGAAATTTCGGCaa                                                                                                                                                                                                                                          <  1:272884/116‑1 (MQ=255)
                                                                                                                                                                              gctgTGCGGCGCATTTGTCGGTGAAACGCCGATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGa                                                                                                                              <  2:235359/224‑1 (MQ=255)
                                                                                                                                                                                    cggcgCATTTGTCGGTGAAACGCCGATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGc                                                                                                                                                               >  1:25246/1‑185 (MQ=255)
                                                                                                                                                                                        gcATTTGTCGGTGAAACGCCGATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGCcaca                                                                                                                                                           <  2:25246/185‑1 (MQ=255)
                                                                                                                                                                                                   tGAAACGCCGATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAg                                                                                                                                                                                                                  >  2:336481/1‑119 (MQ=255)
                                                                                                                                                                                                       aCGCCGATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATCCATCAGCaaa                                                                                                                                                                                                              <  1:336493/119‑1 (MQ=255)
                                                                                                                                                                                                            gATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGCGCGATAGTCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTcgc                                                                                        <  1:175271/232‑1 (MQ=255)
                                                                                                                                                                                                                                                aGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGACATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGt                                                                                                                                                                 >  1:111620/1‑123 (MQ=255)
                                                                                                                                                                                                                                                aGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGc                                                                                                                                                                                          >  2:284446/1‑98 (MQ=255)
                                                                                                                                                                                                                                                    GCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCggg                                                                                                                                                                                      >  1:284457‑M1/1‑98 (MQ=255)
                                                                                                                                                                                                                                                    GCTGGTGCGCGATAATCAGCGTCAGCGCCAGGCGAAAGTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGCca                                                                                                                                                             >  2:111618‑M1/1‑123 (MQ=255)
                                                                                                                                                                                                                                                     CTGGTGCGGGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAgg                                                                                                                                                                                   <  1:159804‑M1/100‑1 (MQ=255)
                                                                                                                                                                                                                                                      TGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGc                                                                                                                                                               >  2:298701‑M1/1‑119 (MQ=255)
                                                                                                                                                                                                                                                       GGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTCTTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTacgac                                                                                                                                                                   <  1:191205‑M1/114‑1 (MQ=255)
                                                                                                                                                                                                                                                       GGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGCCCGGTCCTGTCGTGCGATTTAACCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGACTTGCCCGAATGCGATGNTAAGaaactgacgg  <  2:248077‑M1/275‑3 (MQ=255)
                                                                                                                                                                                                                                                         TGCGGGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAggtggt                                                                                                                                                                               >  2:159800‑M1/1‑100 (MQ=255)
                                                                                                                                                                                                                                                           CGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGc                                                                                                                                                               >  2:191201‑M1/1‑114 (MQ=255)
                                                                                                                                                                                                                                                            GCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATc                                                                              <  1:182508‑M1/194‑1 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                            
GTCAACGCATGTTGATACTTATGGAGCATCCCCGCCAGCCGATTTACCGGTTCCGTCACCTGCGGTGGCGCTTGCCAGATGCGTAGCTTTTCCTGGTAGATTTCCAGTTCTTCCAGCAACTGGCCAAAGTACCGACGACGTTTATCATCGCTACGGGCAGATATCACATGGTCTGCTGTGCGGCGCATTTGTCGGTGAAACGCCGATAAATCCTCGTTAACCGGGATCGGTGCATCACGCAGGCGCTGGTGCGCGATAATCATCGTCAGCGCCAGGCGAAATTTCGGCAAATCCCCTGGGAACTTATTCATCAGCAAAAACAGCTGCTGATAAAGTGCCGGGAGGTGGTTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTG  >  CP009273/3379590‑3380111

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: