Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 392,218 T→C S42P (TCC→CCC)  sbmA → peptide antibiotic transporter SbmA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273392,2180TC96.0% 77.1 / ‑3.8 25S42P (TCC→CCC) sbmApeptide antibiotic transporter SbmA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (12/12);  total (13/12)
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.87e-01

TTTCTCTCGGCCTTCGTTTGGGCATTGATTGCCGTTATCTTCTGGCAAGCCGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTTCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTACGCTTACTACATTGTTTGCGTAGGACTT  >  NZ_CP009273/392131‑392304
                                                                                       |                                                                                      
tttCTCTCGGCCTTCGTTTGGGCATTGATTGCCGTTATCTTCTGGCAAGCCGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTccc                                                                                      <  1:565624/90‑1 (MQ=255)
             tcatttGGGCATTGATTGCCGTTATCTTCTGGCAAGCCGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCa                                                                                 >  2:472410/4‑82 (MQ=255)
             tcatttGGGCATTGATTGCCGTTATCTTCTGGCAAGCCGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCa                                                                                 <  1:472410/79‑1 (MQ=255)
                   gggCATTGATTGCCGTTATCTTCTGGCAAGCCGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTTCCGGGCAGATCCCGATTAgcg                                                                   >  1:126176/1‑90 (MQ=255)
                       attgattgCCGTTATCTTCTGGCAAGCCGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGcgccgc                                                               >  1:834274/1‑90 (MQ=255)
                         tgatttCCGTTATCTTCTGGCAAGCCGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCcgcgc                                                             <  2:836050/90‑1 (MQ=255)
                               ccGTTATCTTCTGGCAAGCCGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCt                                                       <  1:169380/90‑1 (MQ=255)
                               ccGTTATCTTCTGGCAAGCCGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCt                                                       >  1:941237/1‑90 (MQ=255)
                                           ggCAAGCCGGTGGGGGTGACTGGGTGGCGCGTATCCCCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGAtt                                           <  2:343873/90‑1 (MQ=255)
                                              aaGCCGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTcc                                        >  1:299523/1‑90 (MQ=255)
                                              aaGCCGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTcc                                        >  2:695475/1‑90 (MQ=255)
                                               aGCCGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCt                                       >  1:807351/1‑90 (MQ=255)
                                                  cGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGAt                                    >  2:685020/1‑90 (MQ=255)
                                                        gggTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTa                              <  2:299523/90‑1 (MQ=255)
                                                        gggTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTa                              <  1:459963/90‑1 (MQ=255)
                                                          gTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTACg                            <  2:834274/90‑1 (MQ=255)
                                                          gTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTACg                            <  1:695475/90‑1 (MQ=255)
                                                          gTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTACg                            <  1:506352/90‑1 (MQ=255)
                                                          gTGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGATCGTTGGATTTCCTGATTTTTTACg                            <  2:340542/90‑1 (MQ=255)
                                                           tGACTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTACGc                           >  2:552353/1‑90 (MQ=255)
                                                              cTGGGTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTACGCTta                        >  1:1045964/1‑90 (MQ=255)
                                                                 ggTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTACGCTtacta                     >  2:895846/1‑90 (MQ=255)
                                                                  gTGGCGCGTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTACGCTtactac                    >  2:1010187/1‑90 (MQ=255)
                                                                       gcgTATCACCGGCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTACGCTTACTACATTGt               >  1:93174/1‑90 (MQ=255)
                                                                                 ggCGCTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTACGCTTACTACATTGTTTGCGTAGGa     <  1:572677/90‑1 (MQ=255)
                                                                                    gcTCCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTACGCTTACTACATTGTTTGCGTAGGACtt  <  1:729956/90‑1 (MQ=255)
                                                                                       |                                                                                      
TTTCTCTCGGCCTTCGTTTGGGCATTGATTGCCGTTATCTTCTGGCAAGCCGGTGGGGGTGACTGGGTGGCGCGTATCACCGGCGCTTCCGGGCAGATCCCGATTAGCGCCGCGCGTTTCTGGTCGTTGGATTTCCTGATTTTTTACGCTTACTACATTGTTTGCGTAGGACTT  >  NZ_CP009273/392131‑392304

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