Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,429,428 T→A 10.0% V380E (GTA→GAA)  ykpA → ABC superfamily ATP binding cassette transporter binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,429,4280TA10.0% 113.4 / 4.6 40V380E (GTA→GAA) ykpAABC superfamily ATP binding cassette transporter binding protein
Reads supporting (aligned to +/- strand):  ref base T (17/19);  new base A (2/2);  total (19/21)
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.19e-01

GTATTAGATAATATTTCATTCACAATGAATCCAAATGATAAAGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAGACAATACGCACCAGAAGATGAGCAAACTGAAACATTTTTACGCGGCTTCTTAGGCCGTA  >  CP000730/1429289‑1429568
                                                                                                                                           |                                                                                                                                            
gTATTAGATAATATTTCATTCACAATGAATCCAAATGATAAAGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTa                                                                                                                                             >  1:130255/1‑141 (MQ=255)
       ataataTTTCATTCACAATGAATCCAAATGATAAAGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGt                                                                                                                                      <  1:135777/141‑1 (MQ=255)
         aataTTTCATTCACAATGAATCCAAATGATAAAGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCa                                                                                                                                    >  2:161947/1‑141 (MQ=255)
            atTTCATTCACAATGAATCCAAATGATAAAGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTa                                                                                                                                 >  2:64467/1‑141 (MQ=255)
             tttCATTCACAATGAATCCAAATGATAAAGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTaa                                                                                                                                <  2:354622/141‑1 (MQ=255)
               tcattcaCAATGAATCCAAATGATAAAGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGt                                                                                                                              >  1:127323/1‑141 (MQ=255)
                cattcaCAATGAATCCAAATGATAAAGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGtt                                                                                                                             >  2:201583/1‑141 (MQ=255)
                  ttcaCAATGAATCCAAATGATAAAGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGGCGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTAc                                                                                                                           >  2:425583/1‑141 (MQ=255)
                     acaATGAATCCAAATGATAAAGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTc                                                                                                                        <  1:13960/141‑1 (MQ=255)
                      caATGAATCCAAATGATAAAGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTcc                                                                                                                       >  2:369446/1‑141 (MQ=255)
                              ccAAATGATAAAGCGATTTTAATTGGAGATAGTGAAATTGCAAAAACAACATTACTTAAAATATTAGCTGGCGAAATGGAACCAGACGAAGGTTCATTTAAATGGGGTGTTACTACATCATTAAGTTACTTCCCTAAAGAt                                                                                                               >  1:57663/1‑141 (MQ=255)
                                aaaTGATAAAGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATaa                                                                                                             <  1:260207/141‑1 (MQ=255)
                                      tAAAGCAATTTTAATTGGGGATAGTGAAATTGCGACAACCACATTGCTTCAAATATTAGCCGGTGAAATGGAACCACACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGa                                                                                                       <  1:42708/141‑1 (MQ=255)
                                        aaGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAAt                                                                                                     <  2:101733/141‑1 (MQ=255)
                                          gCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATTGCTCAGAATTc                                                                                                   <  2:204369/141‑1 (MQ=255)
                                                   aTTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGc                                                                                          <  2:197769/141‑1 (MQ=255)
                                                     tGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGACCCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCCGAATTCTTTCAGGGCGt                                                                                        >  2:237511/1‑141 (MQ=255)
                                                      ggggATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGtt                                                                                       <  1:37001/141‑1 (MQ=255)
                                                       gggATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTa                                                                                      >  1:137772/1‑141 (MQ=255)
                                                          aTAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAata                                                                                   <  1:175153/141‑1 (MQ=255)
                                                             gTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGa                                                                                >  2:323938/1‑141 (MQ=255)
                                                                    tGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGt                                                                         <  2:41612/141‑1 (MQ=255)
                                                                             cacaTTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGtt                                                                <  1:265255/141‑1 (MQ=255)
                                                                                aTTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAg                                                             >  1:181344/1‑141 (MQ=255)
                                                                                aTTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAg                                                             >  2:201124/1‑141 (MQ=255)
                                                                                         aaTATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAGACAATACGc                                                    <  1:291834/141‑1 (MQ=255)
                                                                                          atatTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAGACAATACGCa                                                   <  1:174101/141‑1 (MQ=255)
                                                                                          atatTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAGACAATACGCa                                                   >  2:183294/1‑141 (MQ=255)
                                                                                           agtcaGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAGACAATACGCAc                                                  <  1:347749/137‑1 (MQ=255)
                                                                                            atTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAGACAATACGCAcc                                                 <  2:281281/141‑1 (MQ=255)
                                                                                              tAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGAAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCa                                                                                                         <  2:402624/83‑1 (MQ=255)
                                                                                              tAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGAAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCa                                                                                                         >  1:402624/1‑83 (MQ=255)
                                                                                              tAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGAAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCa                                                                                                         >  1:402141/1‑83 (MQ=255)
                                                                                              tAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGAAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCa                                                                                                         <  2:402141/83‑1 (MQ=255)
                                                                                                       aaaTGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAGACAATACGCACCAGAAGATGAGc                                      <  2:57495/141‑1 (MQ=255)
                                                                                                         aTGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAGACAATACGCACCAGAAGATGAGCaa                                    >  2:88229/1‑141 (MQ=255)
                                                                                                           ggAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAGACAATACGCACCAGAAGATGAGCAAAc                                  <  2:311174/141‑1 (MQ=255)
                                                                                                                       aGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAGACAATACGCACCAGAAGATGAGCAAACTGAAACAttttt                      >  2:129100/1‑141 (MQ=255)
                                                                                                                              tataAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTGGTAGATTGGTTAAGACAATACGCACCAGAAGATGAGCAAACTGAAACATTTTTACGCGGc               <  1:178613/141‑1 (MQ=255)
                                                                                                                                       ggTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAGACAATACGCACCAGAAGATGAGCAAACTGAAACATTTTTACGCGGCTTCTTAGGc      >  2:184195/1‑141 (MQ=255)
                                                                                                                                           tAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAGACAATACGCACCAGAAGATGAGCAAACTGAAACATTTTTACGCGGCTTCTTAGGCCGTa  <  1:24488/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
GTATTAGATAATATTTCATTCACAATGAATCCAAATGATAAAGCAATTTTAATTGGGGATAGTGAAATTGCGAAAACCACATTGCTTAAAATATTAGCCGGTGAAATGGAACCAGACGAAGGTTCATATAAATGGGGTGTAACAACGTCATTAAGTTACTTCCCTAAAGATAACTCAGAATTCTTTGAGGGCGTTAATATGAACCTTGTAGATTGGTTAAGACAATACGCACCAGAAGATGAGCAAACTGAAACATTTTTACGCGGCTTCTTAGGCCGTA  >  CP000730/1429289‑1429568

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