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,4420TG94.1% 49.7 / ‑3.4 17T141P (ACT→CCT) ybeTtetratricopeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (10/6);  total (11/6)
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.89e-01

GTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCTT  >  NZ_CP009273/674355‑674528
                                                                                       |                                                                                      
gttgttTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGgat                                                                                      <  2:661590/90‑1 (MQ=255)
         cGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagcat                                                                             <  2:615862/90‑1 (MQ=255)
               gTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTgagcatgagcat                                                                       <  2:298431/90‑1 (MQ=255)
                         tttttGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTc                                                             >  1:763272/1‑90 (MQ=255)
                                gaTTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAgctgct                                                      <  2:313386/90‑1 (MQ=255)
                                    cATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCt                                                  <  2:763272/90‑1 (MQ=255)
                                             tAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCaga                                         >  2:15979/1‑90 (MQ=255)
                                                    aaTACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGc                                  >  2:559079/1‑90 (MQ=255)
                                                          tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                            >  2:804685/1‑90 (MQ=255)
                                                          tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                            >  2:548809/1‑90 (MQ=255)
                                                          tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                            >  2:441540/1‑90 (MQ=255)
                                                          tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                            >  2:283161/1‑90 (MQ=255)
                                                          tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                            >  1:332825/1‑90 (MQ=255)
                                                          tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                            >  1:312641/1‑90 (MQ=255)
                                                                        aGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGa              <  1:804685/90‑1 (MQ=255)
                                                                           tCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGAccc           >  2:415655/1‑90 (MQ=255)
                                                                              aaCCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCa        >  1:595222/1‑90 (MQ=255)
                                                                                    aaGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCtt  <  1:182756/90‑1 (MQ=255)
                                                                                       |                                                                                      
GTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCTT  >  NZ_CP009273/674355‑674528

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