Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,350,506 G→T R288S (CGT→AGT)  sapA ← peptide ABC transporter substrate‑binding protein SapA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,350,5060GT92.7% 130.7 / ‑5.3 42R288S (CGT→AGT) sapApeptide ABC transporter substrate‑binding protein SapA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base T (18/20);  minor base C (3/0);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTA  >  NZ_CP009273/1350375‑1350640
                                                                                                                                    |                                                                                                                                      
gATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGCGGGGGTAAATGCCAGCTAGGGGGCGTTCCTCCCCAGGAGCAGCGGTAAAACCAAgccc                                                                                                                                 >  2:499283/1‑137 (MQ=255)
 aTGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAATACCAAAGAGGCGAAGGTCAACCCAGGACGCAGCGTTAAACGCAAGagcg                                                                                                                                >  1:411152/1‑139 (MQ=255)
   ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCAGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCggg                                                                                                                              >  1:1042486/1‑139 (MQ=255)
     aTTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAAAGCCAGATAGGGGACGTTTCTCCCCGGACGCAGCGTTAAAATCAAACGCGGgtc                                                                                                                            >  2:172844/1‑139 (MQ=255)
                        attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAGCCAGATAGGCGACGTTCATCCCCAGACGCAGCGGTAAAATCAAGCGCGGGGCGGCACGCAAAAAGGATAgc                                                                                                         >  1:321374/1‑139 (MQ=255)
                               gCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGATTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAAATCAAGCGCGGGTCGTCACGCAAAATGGATAGctggctgg                                                                                                  >  2:769627/1‑139 (MQ=255)
                                       cAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggc                                                                                          <  2:441397/139‑1 (MQ=255)
                                       cAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGATTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggc                                                                                          <  1:769627/139‑1 (MQ=255)
                                              tGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAggccaggcca                                                                                   <  2:368129/139‑1 (MQ=255)
                                              tGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCAGATAGCCGAAGTTCATCCCCGGACGCCGCGTTAAACCCAAGTGCGGGTCGGCCATCAAAATGGATAGCTGGCTGGCAGcacgctcggc                                                                                     >  1:1003731/1‑132 (MQ=255)
                                                  ggACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAAc                                                                               <  1:1256561/139‑1 (MQ=255)
                                                        cGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAAAGCCACAAAGGGGGCGTTTATCCCCGGACGCCGCGGTAAACCCAAGCGCGGGGGGGCACGGAAAAAGGGAAGCGGGCCGGCAGCCGGCCACGGCCGAAAGGCGCa                                                                         >  2:177237/1‑139 (MQ=255)
                                                          ggATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCAtt                                                                       <  2:411152/139‑1 (MQ=255)
                                                                tttAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCg                                                                  <  1:629174/138‑1 (MQ=255)
                                                                              ggCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGa                                                    <  2:159208/138‑1 (MQ=255)
                                                                                     ttAAATGCTAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTc                                            >  1:214126/1‑139 (MQ=255)
                                                                                     ttAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggccaggc                                                                                     <  1:20413/98‑1 (MQ=255)
                                                                                     ttAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggccaggc                                                                                     >  2:20413/1‑98 (MQ=255)
                                                                                     ttAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGt                                             <  2:984987/138‑1 (MQ=255)
                                                                                     ttAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTc                                            <  2:605186/139‑1 (MQ=255)
                                                                                     ttAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTc                                            >  1:984987/1‑139 (MQ=255)
                                                                                          tGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTg                                       <  2:1003731/139‑1 (MQ=255)
                                                                                             cAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcag                                                                                               <  2:519415/80‑1 (MQ=255)
                                                                                             cAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcag                                                                                               >  1:519415/1‑80 (MQ=255)
                                                                                               gATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgcc                                  >  2:8007/1‑139 (MQ=255)
                                                                                                aTAGGCGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGccgccg                                 <  1:710480/139‑1 (MQ=255)
                                                                                                     cGACGTTCATCCCAGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGcc                            >  1:782346/1‑139 (MQ=255)
                                                                                                                   aGGACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCcaccac                >  2:560615/1‑138 (MQ=255)
                                                                                                                     gACGCAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGGTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTAc             >  1:213698/1‑139 (MQ=255)
                                                                                                                         cAGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGc         <  2:16213/139‑1 (MQ=255)
                                                                                                                          aGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcag                                                                                               <  1:824816/52‑1 (MQ=255)
                                                                                                                          aGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcag                                                                                               >  2:824816/1‑52 (MQ=255)
                                                                                                                          aGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcagg                                                                                           >  2:35342/1‑56 (MQ=255)
                                                                                                                          aGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcagg                                                                                           <  1:35342/56‑1 (MQ=255)
                                                                                                                          aGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCGTTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCg        >  2:12952/1‑139 (MQ=255)
                                                                                                                          aGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCa                                                              >  2:841610/1‑85 (MQ=255)
                                                                                                                          aGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCa                                                              >  2:1216593/1‑85 (MQ=255)
                                                                                                                          aGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCa                                                              <  1:841610/85‑1 (MQ=255)
                                                                                                                          aGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCa                                                              <  1:1216593/85‑1 (MQ=255)
                                                                                                                          aGCGTTAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCg        >  2:1009787/1‑139 (MQ=255)
                                                                                                                              ttAAACTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCAt    <  2:214126/139‑1 (MQ=255)
                                                                                                                                 aaCTCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTa  <  1:1009787/138‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
GATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTA  >  NZ_CP009273/1350375‑1350640

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