Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,451,709 G→T S310R (AGC→AGA fadJ ← fatty acid oxidation complex subunit alpha FadJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,451,7090GT95.2% 148.0 / ‑0.2 42S310R (AGC→AGAfadJfatty acid oxidation complex subunit alpha FadJ
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (24/16);  total (25/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 = 8.25e-01

TACTTCAGCGCATGATTTATGCCCTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACGCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCAGCGCCTGCGATTGTGGCGTCATCGCCAGTTCGCCAAACGCCCGAGCTTCGGCGTCATAA  >  NZ_CP009273/2451587‑2451844
                                                                                                                          |                                                                                                                                       
tACTTCAGCGCATGATTTATGCCCTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAg                                                                                                                         >  2:943755/1‑139 (MQ=255)
  cTTCAGCGCATGATTTATGCCCTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAggc                                                                                                                       >  1:400609/1‑139 (MQ=255)
  cTTCAGCGCATGATTTATGCCCTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAgg                                                                                                                        <  2:400609/138‑1 (MQ=255)
     cAGCGCATGATTTATGCCCTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAggcgg                                                                                                                     <  1:1064073/138‑1 (MQ=255)
      aGCGCATGATTTATGCCCTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGgcggcg                                                                                                                   <  1:822620/139‑1 (MQ=255)
        cgcATGATTTATGCCCTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCa                                                                                                                 >  1:1164668/1‑139 (MQ=255)
        cgcATGATTTATGCCCTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCa                                                                                                                 >  1:1246521/1‑139 (MQ=255)
        cgcATGATTTATGCCCTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCa                                                                                                                 >  2:479516/1‑139 (MQ=255)
        cgcATGATTTATGCCCTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCa                                                                                                                 >  1:92889/1‑139 (MQ=255)
                       cTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGAtct                                                                                                  >  2:99516/1‑139 (MQ=255)
                       cTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGAtct                                                                                                  <  2:1164668/139‑1 (MQ=255)
                       cTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGAtct                                                                                                  >  2:100223/1‑139 (MQ=255)
                        tGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGAtctt                                                                                                 >  1:1091855/1‑139 (MQ=255)
                            ggTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGAtctttctt                                                                                             <  2:981084/139‑1 (MQ=255)
                            ggTTGATATCTTTAATTAGGAACGGAATCCCCGCTTTACAAGCAGTCACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGAtctttctt                                                                                             <  2:553334/139‑1 (MQ=255)
                                                aCCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCaaaa                                                                          >  2:237672/1‑138 (MQ=255)
                                                aCCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCaaaa                                                                          >  1:122012/1‑138 (MQ=255)
                                                aCCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCaaa                                                                           <  1:195002/137‑1 (MQ=255)
                                                aCCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACGCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCaaaaa                                                                         <  1:152836/139‑1 (MQ=255)
                                                          ccGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTCCTGGCAAAAAAGATACTACg                                                               >  2:1286118/1‑139 (MQ=255)
                                                          ccGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACg                                                               <  1:943755/139‑1 (MQ=255)
                                                            gCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACg                                                               >  2:360451/1‑137 (MQ=255)
                                                                  cAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCAGCGCCt                                                       <  1:340246/139‑1 (MQ=255)
                                                                     gCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCAGCGCCTGc                                                     >  2:986508/1‑138 (MQ=255)
                                                                        gTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCAGCGCCTGCGAtt                                                 <  2:340143/139‑1 (MQ=255)
                                                                                   aaTACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCa                                                             >  1:32446/1‑116 (MQ=255)
                                                                                   aaTACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCa                                                             <  2:32446/116‑1 (MQ=255)
                                                                                     tACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCAGCGCCTGCGATTGTGGCGTCATCGc                                    >  1:903602/1‑139 (MQ=255)
                                                                                     tACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCAGCGCCTGCGATTGTGGCGTCATCGc                                    >  2:794646/1‑139 (MQ=255)
                                                                                     tACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCAGCGCCTGCGATTGTGGCGTCATCGc                                    >  1:638684/1‑139 (MQ=255)
                                                                                     tACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCAGCGCCTGCGATTGTGGCGTCATCGc                                    >  2:1069010/1‑139 (MQ=255)
                                                                                     tACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACGCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCAGCGCCTGCGATTGTGGCGTCATCGc                                    >  1:240873/1‑139 (MQ=255)
                                                                                       ccgccgCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCa                                                             >  1:1149665/1‑112 (MQ=255)
                                                                                       ccgccgCCCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCa                                                             <  2:1149665/112‑1 (MQ=255)
                                                                                          ccgccCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGTCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCAGCGCCTGCGATTGTGGCGTCATCGCCAGtt                               >  1:274269/1‑139 (MQ=255)
                                                                                          ccgccCATCAAGCCACCACCTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGc                                                                                                          <  1:365935/64‑1 (MQ=255)
                                                                                          ccgccCATCAAGCCACCACCTAAAATCCCCACTATGTTTAATGGCGCAGGCGGCGCATCACTGc                                                                                                          >  2:365935/1‑64 (MQ=255)
                                                                                                         ccaccTAAAATCCCCACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCAGCGCCTGCGATTGTGGCGTCATCGCCAGTTCGCCAAACGCCCGAg                <  1:971440/139‑1 (MQ=255)
                                                                                                                     cccACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGAtctttctt                                                                                             >  2:372356/1‑50 (MQ=255)
                                                                                                                     cccACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGAtctttctt                                                                                             <  1:372356/50‑1 (MQ=255)
                                                                                                                     cccACTCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCAGCGCCTGCGATTGTGGCGTCATCGCCAGTTCGCCAAACGCCCGAGCTTCGGCGTCAt    <  1:360451/139‑1 (MQ=255)
                                                                                                                       cACTCTGTTTAATGGCGCAGACGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAAATACTACGCAGCGCCTGCGATTGTGGCGTCATCGCCAGTTCGCCAAACGCCCGAGCGTCGGCGTCATaa  >  1:92666/1‑139 (MQ=255)
                                                                                                                          |                                                                                                                                       
TACTTCAGCGCATGATTTATGCCCTGCGGGTTGATATCTTTAATTCTGACCGGAATCCCCGCTTTACAAGCAGTGACATAAGCAATACCGCCGCCCATCAAGCCACCACCTAAAATCCCCACGCTGTTTAATGGCGCAGGCGGCGCATCACTGCCGGGATCTTTCTTCACGTCCGTACTGGCAAAAAAGATACTACGCAGCGCCTGCGATTGTGGCGTCATCGCCAGTTCGCCAAACGCCCGAGCTTCGGCGTCATAA  >  NZ_CP009273/2451587‑2451844

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