Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,559,065 C→T 10.8% A157A (GCC→GCT)  bcr → MFS family major facilitator transporter, bicyclomycin:cation symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,559,0650CT10.8% 113.5 / 3.6 43A157A (GCC→GCT) bcrMFS family major facilitator transporter, bicyclomycin:cation symporter
Reads supporting (aligned to +/- strand):  ref base C (23/15);  new base T (3/2);  total (26/17)
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.32e-01

ATTCAAGGTTTAACTGGTGGTGGCGTCATCGTGATTGCAAAAGCCTCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAATTACCTAAAACATCTAAAGATGAATTAAAGCAGGTG  >  CP000730/2558925‑2559194
                                                                                                                                            |                                                                                                                                 
aTTCAAGGTTTAACTGGTGGTGGCGTCATCGTGATTGCAAAAGCCTCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGcc                                                                                                                                   <  1:246021/141‑1 (MQ=255)
aTTCAAGGTTTAACTGGTGGTGGCGTCATCGTGATTGCAAAAGCCTCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGcc                                                                                                                                   >  2:54233/1‑141 (MQ=255)
            aCTGGTGGTGGCGTCATCGTGATTGCAAAAGCCTCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCt                                                                                                                       >  1:360226/1‑141 (MQ=255)
             cTGGTGGTGGCGTCATCGTGATTGCAAAAGCCTCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTCGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCtt                                                                                                                      >  1:248531/1‑141 (MQ=255)
                     ggCGTCATCGTGATTGCAAAAACCTCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTa                                                                                                              <  2:26540/141‑1 (MQ=255)
                      gCGTCATCGTGATTGCCAAAGCCTCTGCTGGTGATAAATTTAGTGGCAACGCACTAGCTAAATTTTTAGCATCTTTAATGGTAGTAAATGGCATCATCACTATTCTTGCACCATTAGCTGGTGGATTAGCTTTATCCGTAg                                                                                                             >  1:218988/1‑141 (MQ=255)
                         tCATCGTGATTGCAAAAGCCTCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCaa                                                                                                          >  2:146645/1‑141 (MQ=255)
                            tCGTGATTGCAAAAGCCTCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTTATGGTAGTTAATGGCAGCATCACTATTCTTGCACCATTAGCTGGGGGATTAGCTTTATCCGTAGCAActt                                                                                                       >  2:130326/1‑139 (MQ=255)
                            tCGTGATTGCAAAAGCCTCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACAt                                                                                                       <  2:13564/141‑1 (MQ=255)
                                 aTTGCAAAAGCCTCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGt                                                                                                  >  1:51778/1‑141 (MQ=255)
                                 aTTGCAAAAGCATCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATCTTTAGAATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTGACCGATGGATTAGCTTTATCCGTAGCAACATGTCga                                                                                                  >  2:293719/1‑140 (MQ=255)
                                    gccaaaGCCTCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTAAATGGCATCATCACTATTCTTGCACCATTAGCTGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCt                                                                                               >  1:377669/4‑141 (MQ=255)
                                     cAAAAGCCTCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTa                                                                                              >  1:274095/1‑141 (MQ=255)
                                       aaaGCCTCTGCTGGTGATAAATTTAGTGTCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTAtt                                                                                            <  2:455574/141‑1 (MQ=255)
                                               tgctgGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAAt                                                                                    >  1:228854/1‑141 (MQ=255)
                                                gctgGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAAtt                                                                                   >  1:250827/1‑141 (MQ=255)
                                                gctgGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAAtt                                                                                   >  2:58467/1‑141 (MQ=255)
                                                    gTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTaa                                                                               >  1:158189/1‑141 (MQ=255)
                                                     tGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAAc                                                                              >  2:314864/1‑141 (MQ=255)
                                                         aaaTTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTAtt                                                                          >  1:110597/1‑141 (MQ=255)
                                                          aaTTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATtg                                                                         <  1:428412/141‑1 (MQ=255)
                                                            tttAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATtgtg                                                                       >  2:23430/1‑141 (MQ=255)
                                                              tAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGc                                                                     >  2:106490/1‑141 (MQ=255)
                                                               aGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCTCAATTTTAACTATTGTGGct                                                                    >  2:260280/1‑140 (MQ=255)
                                                                gTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCAc                                                                   <  1:46000/141‑1 (MQ=255)
                                                                 tGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATCTTCACAATTTTAACTATTGTGGCACt                                                                  >  2:272419/1‑141 (MQ=255)
                                                                     aaCGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACtcatc                                                              <  1:122194/141‑1 (MQ=255)
                                                                                 aaaTTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGc                                                  <  1:227995/141‑1 (MQ=255)
                                                                                   aTTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGt                                                <  1:146915/141‑1 (MQ=255)
                                                                                            aTCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCa                                       >  1:158638/1‑141 (MQ=255)
                                                                                            aTCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCa                                       >  1:158749/1‑141 (MQ=255)
                                                                                              cTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAAt                                     <  2:53370/141‑1 (MQ=255)
                                                                                                   aTGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAATTACCt                                <  2:54842/141‑1 (MQ=255)
                                                                                                   aTGGTAGTAAATGGCATCATCACTATTCTTGCACCATTAGCTGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAATTACCt                                <  2:377669/141‑1 (MQ=255)
                                                                                                      gtagtaAATGGCATCATCACTATTCTTGCACCATTAGCTGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAATTACCTaaa                             <  2:144278/141‑1 (MQ=255)
                                                                                                      gTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAATTACCTaaa                             <  1:436426/141‑1 (MQ=255)
                                                                                                        aGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTTTCAATTACCTAAAAc                           >  1:136432/1‑141 (MQ=255)
                                                                                                        aGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAATTACCTAAAAc                           <  1:330452/141‑1 (MQ=255)
                                                                                                        aGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAATTACCTAAAAc                           <  1:146645/141‑1 (MQ=255)
                                                                                                            aaTGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATACGTAGCCACTTGGCGTTCTTCTCCCACAACTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAATTACCTAAAACATCt                       <  1:317187/141‑1 (MQ=255)
                                                                                                            aaTGGCATCATCACTATGCGGGCAGCAGTATCCTTGGGATTAGCGTGAGCCGGAGCTACAGGGCGTTCTATTTTCACAATCTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAATTACCTAAAACATCt                       <  1:384861/141‑1 (MQ=255)
                                                                                                                       caCTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAATTACCTAAAACATCTAAAGATGAAtt            >  2:48862/1‑141 (MQ=255)
                                                                                                                             tCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAATTACCTAAAACATCTAAAGATGAATTAAAGCa      >  2:381989/1‑141 (MQ=255)
                                                                                                                               ttGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAATTACCTAAAACATCTAAAGATGAATTAAAGCAgg    >  2:78484/1‑141 (MQ=255)
                                                                                                                                 gCACCGTTAGCCGGTGGATTCGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTGATTTGCGTCGCTTCTCAATTACCTAAAACATCTAACGATGAATTATAGCAGGTg  >  2:151946/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                 
ATTCAAGGTTTAACTGGTGGTGGCGTCATCGTGATTGCAAAAGCCTCTGCTGGTGATAAATTTAGTGGCAACGCACTCGCTAAATTTTTAGCATCTTTAATGGTAGTTAATGGCATCATCACTATTCTTGCACCATTAGCCGGTGGATTAGCTTTATCCGTAGCAACATGGCGTTCTATTTTCACAATTTTAACTATTGTGGCACTCATCATTTTAATTGGCGTCGCTTCTCAATTACCTAAAACATCTAAAGATGAATTAAAGCAGGTG  >  CP000730/2558925‑2559194

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