Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 176,296 C→G R191R (CGC→CGG dgt → dGTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273176,2960CG76.0% 51.8 / 15.0 25R191R (CGC→CGGdgtdGTPase
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base G (10/9);  total (13/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 = 1.00e+00

ACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTG  >  NZ_CP009273/176209‑176381
                                                                                       |                                                                                     
acgggacgggGAAGAGCCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGGCt                                                                                     <  1:184453/90‑1 (MQ=255)
  gggacgggGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGGCTgg                                                                                   <  1:495940/90‑1 (MQ=255)
      cgggGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGGCTGGTGCa                                                                               >  2:80415/1‑90 (MQ=255)
      cgggGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCa                                                                               <  1:521530/90‑1 (MQ=255)
        gggAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGGCTGGTGCATa                                                                             <  2:279539/90‑1 (MQ=255)
             gaaCCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTg                                                                        >  2:91925/1‑90 (MQ=255)
                  gCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGGCTGGTGCATACATTGATGCg                                                                   >  1:248523/1‑90 (MQ=255)
                     tAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGGCTGGTGCATACATTGATGCGGAt                                                                <  2:477056/90‑1 (MQ=255)
                          aGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGGCTGGTGCATACATTGATGCGGATGAAtc                                                           >  2:552275/1‑90 (MQ=255)
                               cggcgCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGGCTGGTGCATACATTGATGCGGATGAATCTCAcc                                                      <  2:429130/90‑1 (MQ=255)
                                              cAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGGCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGc                                       <  1:80415/90‑1 (MQ=255)
                                                     tATGTCATTTTGAGGGGAATGCACAAGGCATTCGGCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTAttt                                >  1:9802/1‑90 (MQ=255)
                                                      aTGTCATTTTGAGGGGAATGCACAAGGCATTCGGCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTAtttt                               <  1:434538/90‑1 (MQ=255)
                                                      aTGTCATTTTGAGGGGAATGCACAAGGCATTCGGCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTAtttt                               <  2:173220/90‑1 (MQ=255)
                                                           aTTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGCTGGCGGTATTTTAAAat                          >  1:296241/1‑90 (MQ=255)
                                                             tttGAGGGGAATGCACAAGGCATTCGGCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAatat                        >  1:340764/1‑90 (MQ=255)
                                                               tGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATAc                      <  1:224776/90‑1 (MQ=255)
                                                                   gggAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGt                  <  1:20768/90‑1 (MQ=255)
                                                                     gAATGCACAAGGCATTCGGCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTcc                >  1:26541/1‑90 (MQ=255)
                                                                              aGGCATTCGGCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCgtggtg       >  1:73106/1‑90 (MQ=255)
                                                                                  aTTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGgtggcgt   >  1:246239/1‑90 (MQ=255)
                                                                                   ttCGGCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCgg              >  1:564555/1‑78 (MQ=255)
                                                                                   ttCGGCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCgg              <  2:564555/78‑1 (MQ=255)
                                                                                   ttCGGCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGgtggcgtg  >  1:102403/1‑90 (MQ=255)
                                                                                   ttCGGCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGgtggcgtg  >  1:80740/1‑90 (MQ=255)
                                                                                       |                                                                                     
ACGGGACGGGGAAGAACCGCTTAACGAGCTGCGGCGCAAGATTCGTCAGGACTTATGTCATTTTGAGGGGAATGCACAAGGCATTCGCCTGGTGCATACATTGATGCGGATGAATCTCACCTGGGCACAGGTTGGCGGTATTTTAAAATATACCCGTCCGGCGTGGTGGCGTG  >  NZ_CP009273/176209‑176381

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