Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 674,442 T→G T141P (ACT→CCT)  ybeT ← tetratricopeptide repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273674,4420TG90.9% 64.5 / ‑0.1 22T141P (ACT→CCT) ybeTtetratricopeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (9/11);  total (10/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 = 9.91e-01

GTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGA  >  NZ_CP009273/674355‑674524
                                                                                       |                                                                                  
gttgttTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGgat                                                                                  >  1:125518/1‑90 (MQ=255)
  tgttTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGatat                                                                                <  1:89587/90‑1 (MQ=255)
   gttTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATAtt                                                                               >  2:330667/1‑90 (MQ=255)
        gCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagca                                                                          <  1:65691/90‑1 (MQ=255)
         cGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagcat                                                                         <  2:465837/90‑1 (MQ=255)
               gTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTgagcatgagcat                                                                   <  2:125518/90‑1 (MQ=255)
                         tttttGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTc                                                         >  1:301460/1‑90 (MQ=255)
                                gaTTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAgctgct                                                  <  1:333886/90‑1 (MQ=255)
                                                  ttAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAAtt                                <  1:220960/90‑1 (MQ=255)
                                                         cTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTAc                         >  1:133239/1‑90 (MQ=255)
                                                          tGAGCCATTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                        >  2:262119/1‑90 (MQ=255)
                                                          tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                        >  2:19000/1‑90 (MQ=255)
                                                          tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                        >  2:295321/1‑90 (MQ=255)
                                                           gAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACgg                       >  2:202483/1‑90 (MQ=255)
                                                                aTTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTc                  <  1:262119/90‑1 (MQ=255)
                                                                 tttCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCg                 >  1:68756/1‑90 (MQ=255)
                                                                        aGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGa          <  2:117879/90‑1 (MQ=255)
                                                                              aaCCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAAtt                                <  2:152035/62‑1 (MQ=255)
                                                                              aaCCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAAtt                                >  1:152035/1‑62 (MQ=255)
                                                                              aaCCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTg           <  1:154363/83‑1 (MQ=255)
                                                                              aaCCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTg           >  2:154363/1‑83 (MQ=255)
                                                                               aCCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAg   <  1:437700/90‑1 (MQ=255)
                                                                                ccTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGa  <  2:301460/90‑1 (MQ=255)
                                                                                       |                                                                                  
GTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGA  >  NZ_CP009273/674355‑674524

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