Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,149,506 T→C S251P (TCC→CCC)  pabC → aminodeoxychorismate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,149,5060TC96.7% 93.1 / ‑3.5 30S251P (TCC→CCC) pabCaminodeoxychorismate lyase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (20/9);  total (20/10)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

GTTTGCTGGCACAATCCTCTTATCAGCTTGTCGAAGTGCAAGCCTCTCTGGAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCG  >  NZ_CP009273/1149375‑1149634
                                                                                                                                   |                                                                                                                                
gTTTGCTGGCACAATCCTCTTATCAGCTTGTCGAAGTGCAAGCCTCTCTGGAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTc                                                                                                                           >  2:86216/1‑139 (MQ=255)
          acaATCCTCTTATCAGCTTGTCGAAGTGCACGCCTCTCTGGAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACg                                                                                                                 <  2:242865/139‑1 (MQ=255)
              tCCTCTTATCAGCTTGTCGAAGTGCAAGCCTCTCTGGAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTta                                                                                                              >  2:588634/1‑138 (MQ=255)
                          cTTGTCGAAGTGCAAGCCTCTCTGGAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAg                                                                                                 >  1:143269/1‑139 (MQ=255)
                                     gCAAGCCTCTCTGGAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTtgt                                                                                      >  1:255264/1‑139 (MQ=255)
                                           ctctctGGAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAgcgc                                                                                >  1:328963/1‑139 (MQ=255)
                                              tctGGAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGccc                                                                              >  1:399446/1‑138 (MQ=255)
                                                     gagTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGt                                                                      >  2:220152/1‑139 (MQ=255)
                                                        tCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGc                                                                                                                     >  2:604549/1‑89 (MQ=255)
                                                        tCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGc                                                                                                                     <  1:604549/89‑1 (MQ=255)
                                                          gTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATg                                                                  >  2:552304/1‑138 (MQ=255)
                                                                     atgagatgGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTa                                                      >  1:122803/1‑139 (MQ=255)
                                                                     atgagatgGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTa                                                      >  1:156212/1‑139 (MQ=255)
                                                                                 tttGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTg                                          >  2:385547/1‑139 (MQ=255)
                                                                                     tAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTa                                       <  2:858710/138‑1 (MQ=255)
                                                                                       aTGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTAtt                                     <  1:552304/138‑1 (MQ=255)
                                                                                       aTGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTg                                    >  1:811140/1‑139 (MQ=255)
                                                                                         gCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCt                                  >  1:834720/1‑139 (MQ=255)
                                                                                          cGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTAtt                                     >  1:632207/1‑135 (MQ=255)
                                                                                          cGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTAtt                                     <  2:632207/135‑1 (MQ=255)
                                                                                              aaTGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCtggtg                              <  1:826846/138‑1 (MQ=255)
                                                                                                tGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTa                           <  1:220152/139‑1 (MQ=255)
                                                                                                     gTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAAt                                              >  1:106711/1‑115 (MQ=255)
                                                                                                     gTGATGCCCGTATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAAt                                              <  2:106711/115‑1 (MQ=255)
                                                                                                               tATGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCt            >  2:814598/1‑139 (MQ=255)
                                                                                                                aTGTGCCTGTGGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTg           >  2:422846/1‑139 (MQ=255)
                                                                                                                         tgGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTggg    <  1:223694/137‑1 (MQ=255)
                                                                                                                         tgGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTggg    <  1:393777/137‑1 (MQ=255)
                                                                                                                         tgGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCg  >  2:817367/1‑139 (MQ=255)
                                                                                                                         tgGCGATGTCCCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCg  >  1:609310/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                
GTTTGCTGGCACAATCCTCTTATCAGCTTGTCGAAGTGCAAGCCTCTCTGGAAGAGTCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCG  >  NZ_CP009273/1149375‑1149634

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