Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,443,731 C→T A111V (GCC→GTC)  yjfF → sugar ABC transporter permease YjfF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,443,7310CT97.7% 152.5 / ‑3.6 44A111V (GCC→GTC) yjfFsugar ABC transporter permease YjfF
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (23/20);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

CGGTGATCGCCTTTACTGGCGTGTTTCTGGCAAAAGTGATTGGCGATTTCGGCCTCTCGCCGCTGCTGGCGTTTCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGCCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAA  >  NZ_CP009273/4443599‑4443865
                                                                                                                                    |                                                                                                                                      
cggtGATCGCCTTTACTGGCGTGTTTCTGGCAAAAGTGATTGGCGATTTCGGCCTCTCGCCGATGCTGGCGTTTCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCAGTTATGGGGCTTCTGATAGACGTCCTGac                                                                                                                                  >  2:864563/1‑138 (MQ=255)
 gttGATCGCCTTTCCTGGCGTGTTTCTGGAAAAAGTGATTGGCGATTTCGGCCTCTCGCCGCTGCTGGCGTTTCCGCTAGTGCTGGTGTTGGGCTGTGCCTTCGTCGCATTTATGGGGCTTCTGATTGACGTCCTGaa                                                                                                                                  <  2:910980/136‑1 (MQ=255)
     aTCGCCTTTACTGGCGTGTTTCTGGCAAAAGTGATTGGCGATTTCGGCCTCTCGCCGCTGCTGGCGTTTCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATcc                                                                                                                             <  2:631652/139‑1 (MQ=255)
     aTCGCCTTTACTGGCGTGTTTCTGGCAAAAGTGATTGGCGATTTCGGCCTCTCGCCGCTGCTGGCGTTTCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATcc                                                                                                                             <  2:1095616/139‑1 (MQ=255)
     aTCGCCTTTACTGGCGTGTTTCTGGCAAAAGTGATTGGCGATTTCGGCCTCTCGCCGCTGCTGGCGTTTCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATc                                                                                                                              <  1:424128/138‑1 (MQ=255)
                  gCGTGTTTCTGGCAAAAGTGATTGGCGATTTCGGCCTCTCGCCGCTGCTGGCGTTTCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCAt                                                                                                                >  1:972054/1‑139 (MQ=255)
                    gtgtTTCTGGCAAAAGTGATTGGCGATTTCGGCCTCTCGCCGCTGCTGGCGTTTCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTa                                                                                                              >  2:1250326/1‑139 (MQ=255)
                       tttCTGGCAAAAGTGATTGGCGATTTCGGCCTCTCGCCGCTGCTGGCGTTTCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGc                                                                                                           >  1:922366/1‑139 (MQ=255)
                       tttCTGGCAAAAGTGATTGGCGATTTCGGCCTCTCGCCGCTGCTGGCGTTTCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGc                                                                                                           >  1:539190/1‑139 (MQ=255)
                             gCAAAAGTGATTGGCGATTTCGGCCTCTCGCCGCTGCTGGCGTTTCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCg                                                                                                      <  2:655516/138‑1 (MQ=255)
                             gCAAAAGTGATTGGCGATTTCGGCCTCTCGCCGCTGCTGGCGTTTCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGCCCTGAAGATCCCGGCATTTATCATTACGCTTGCg                                                                                                      <  1:875121/138‑1 (MQ=255)
                                                                         tCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTAc                                                                                                             <  2:1136580/87‑1 (MQ=255)
                                                                         tCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTAc                                                                                                             >  1:1136580/1‑87 (MQ=255)
                                                                             cTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGa                                                     <  1:46278/139‑1 (MQ=255)
                                                                             cTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGa                                                     <  1:341556/139‑1 (MQ=255)
                                                                                gTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGAt                                                    <  1:99022/137‑1 (MQ=255)
                                                                                gTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGAt                                                    <  2:1174487/137‑1 (MQ=255)
                                                                                 tGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTACGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTcc                                                 >  1:908174/1‑139 (MQ=255)
                                                                                     ggTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGAt                                              >  1:1177157/1‑138 (MQ=255)
                                                                                      gTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATa                                             >  1:664616/1‑138 (MQ=255)
                                                                                       tGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATaaa                                           <  2:924795/139‑1 (MQ=255)
                                                                                       tGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATaaa                                           >  1:273822/1‑139 (MQ=255)
                                                                                          tGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCAtt                                                                                                                      <  2:14881/61‑1 (MQ=255)
                                                                                          tGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCAtt                                                                                                                      >  1:14881/1‑61 (MQ=255)
                                                                                            ggCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCAt                                       <  2:1070998/138‑1 (MQ=255)
                                                                                                gtgCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGa                                                             >  1:539306/1‑112 (MQ=255)
                                                                                                gtgCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGa                                                             <  2:539306/112‑1 (MQ=255)
                                                                                                gtgCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCAt                                  >  1:840810/1‑139 (MQ=255)
                                                                                                 tgCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCAtt                                 >  1:655046/1‑139 (MQ=255)
                                                                                                 tgCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCAtt                                 >  1:1187809/1‑139 (MQ=255)
                                                                                                 tgCCTTCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCAtt                                 >  2:268074/1‑139 (MQ=255)
                                                                                                      tCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGc                                                                                                                         <  1:66457/46‑1 (MQ=255)
                                                                                                      tCGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGc                                                                                                                         >  2:66457/1‑46 (MQ=255)
                                                                                                       cGGCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGa                            >  2:44648/1‑138 (MQ=255)
                                                                                                        ggCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGaca                          >  1:793251/1‑139 (MQ=255)
                                                                                                        ggCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGaca                          <  2:857541/139‑1 (MQ=255)
                                                                                                        ggCGCATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGaca                          <  1:444515/139‑1 (MQ=255)
                                                                                                          cgcATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGACACg                        <  2:426640/139‑1 (MQ=255)
                                                                                                          cgcATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGACACg                        <  1:650923/139‑1 (MQ=255)
                                                                                                          cgcATTTATGGGGCTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGACGATTCCGATAAACCATCCCATTTATGACACg                        <  1:1105638/139‑1 (MQ=255)
                                                                                                                       cTTCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTg           >  1:775860/1‑139 (MQ=255)
                                                                                                                         tCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCg         >  1:1204580/1‑139 (MQ=255)
                                                                                                                         tCTGATCGACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCg         >  2:730331/1‑139 (MQ=255)
                                                                                                                                gACGTCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGaaaa  >  1:1277828/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                      
CGGTGATCGCCTTTACTGGCGTGTTTCTGGCAAAAGTGATTGGCGATTTCGGCCTCTCGCCGCTGCTGGCGTTTCCGCTAGTGCTGGTGATGGGCTGTGCCTTCGGCGCATTTATGGGGCTTCTGATCGACGCCCTGAAGATCCCGGCATTTATCATTACGCTTGCGGGGATGTTCTTTTTGCGCGGCGTCAGCTATCTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAA  >  NZ_CP009273/4443599‑4443865

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