Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,579,168 T→C E339G (GAA→GGA)  ggt ← gamma‑glutamyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,579,1680TC90.9% 82.9 / 2.5 33E339G (GAA→GGA) ggtgamma‑glutamyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (15/15);  total (16/17)
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.08e-01

CAATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTC  >  NZ_CP009273/3579082‑3579248
                                                                                       |                                                                                
cAATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCg                                                                               <  1:119847/90‑1 (MQ=255)
      gATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGt                                                                         >  1:644871/1‑90 (MQ=255)
       agcGGCAATAGATTTGGCATCGGCTTTACTGGGCCGCGCCTGCCACGGGACTATGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTc                                                                        <  2:783600/88‑1 (MQ=255)
       accGGCAATAGATTTGGCATAGGCTTTCTTGGTCAGCGCCTGCCACGGTACTTTGCCAAAATCCGGGTCGCCAAGATATCCCGAGCGGTc                                                                        <  2:544987/88‑1 (MQ=255)
       aTCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTc                                                                        <  1:526972/90‑1 (MQ=255)
             aaTAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCggcgta                                                                  <  1:758588/90‑1 (MQ=255)
                     tGGCATAGGCTTTATTGGTCAGCGCCTGCTACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTAtt                                                          >  2:264963/1‑90 (MQ=255)
                     tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTAtt                                                          >  1:88969/1‑90 (MQ=255)
                     tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTAtt                                                          >  2:88592/1‑90 (MQ=255)
                     tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTAtt                                                          >  2:21527/1‑90 (MQ=255)
                     tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTAtt                                                          >  1:682749/1‑90 (MQ=255)
                     tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTATGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTAtt                                                          >  1:512063/1‑90 (MQ=255)
                         aTAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTtctc                                                      >  2:346723/1‑90 (MQ=255)
                          tAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTcc                                                     <  1:531231/90‑1 (MQ=255)
                               tttATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTc                                                <  1:399028/90‑1 (MQ=255)
                                  aTTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGc                                             <  2:512063/90‑1 (MQ=255)
                                    tGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCa                                           <  2:301210/90‑1 (MQ=255)
                                     ggTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCAt                                          >  1:482852/1‑90 (MQ=255)
                                       tCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGa                                        <  2:644871/90‑1 (MQ=255)
                                        cAGCGCCTGCCACGGTACTTTGGCAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGAt                                       >  1:238363/1‑90 (MQ=255)
                                               tGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTgcatc                                >  2:73207/1‑90 (MQ=255)
                                                        aCTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGGGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGgcg                       <  2:521968/90‑1 (MQ=255)
                                                        aCTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGgcg                       <  2:174689/90‑1 (MQ=255)
                                                         cTTTGACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGgcgc                      >  1:549200/1‑90 (MQ=255)
                                                             gACAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGcc                  >  1:707119/1‑90 (MQ=255)
                                                              aCAAAATCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCa                 <  1:531772/90‑1 (MQ=255)
                                                                  aaTCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTc                                                <  1:828354/55‑1 (MQ=255)
                                                                  aaTCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTc                                                >  2:828354/1‑55 (MQ=255)
                                                                  aaTCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGc             <  2:346186/90‑1 (MQ=255)
                                                                  aaTCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGc             >  2:494425/1‑90 (MQ=255)
                                                                  aaTCCGGGTCGCCAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGc             >  2:202777/1‑90 (MQ=255)
                                                                             ccAAGATATCCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATttctt   <  1:629014/90‑1 (MQ=255)
                                                                             ccAAGATATCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATttcttc  <  1:734460/90‑1 (MQ=255)
                                                                                       |                                                                                
CAATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATATTCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTC  >  NZ_CP009273/3579082‑3579248

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