Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,803,201 A→G V97A (GTC→GCC)  yniB ← YniB family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,803,2010AG78.3% 35.2 / 9.7 23V97A (GTC→GCC) yniBYniB family protein
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base G (8/10);  total (9/14)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

ATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGAATAGAT  >  NZ_CP009273/1803120‑1803285
                                                                                 |                                                                                    
aTGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATc                                                                              <  2:598508/90‑1 (MQ=255)
           ttCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAAt                                                                   <  1:47771/90‑1 (MQ=255)
           ttCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAAt                                                                   <  2:595613/90‑1 (MQ=255)
            tCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATc                                                                  >  2:734063/1‑90 (MQ=255)
                cacCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCAccc                                                              <  2:110263/90‑1 (MQ=255)
                    ttCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGaa                                                          >  2:574189/1‑90 (MQ=255)
                    ttCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGAAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGaa                                                          >  2:800218/1‑90 (MQ=255)
                     tCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAAc                                                         >  2:303115/1‑90 (MQ=255)
                                   cgcCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGcc                                           <  2:71500/90‑1 (MQ=255)
                                    gcCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGcca                                          <  2:94415/90‑1 (MQ=255)
                                    gcCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGcca                                          <  2:836384/90‑1 (MQ=255)
                                                 tCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTgg                             >  1:620052/1‑90 (MQ=255)
                                                  ccGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTggg                            <  1:765995/90‑1 (MQ=255)
                                                         ccGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCggg                     <  2:290195/90‑1 (MQ=255)
                                                         ccGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCggg                     <  2:59430/90‑1 (MQ=255)
                                                           ggAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCggg                     <  1:622659/88‑1 (MQ=255)
                                                           ggAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCggg                     >  2:622659/1‑88 (MQ=255)
                                                                  tGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATAc            <  2:306160/90‑1 (MQ=255)
                                                                   gCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGtt                >  1:183654/1‑85 (MQ=255)
                                                                   gCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGtt                <  2:183654/85‑1 (MQ=255)
                                                                          cAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCggg                     <  2:322602/73‑1 (MQ=255)
                                                                          cAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCggg                     >  1:322602/1‑73 (MQ=255)
                                                                            gTCCGGCAATAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGAATAGAt  >  2:261545/1‑90 (MQ=255)
                                                                                 |                                                                                    
ATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGAATAGAT  >  NZ_CP009273/1803120‑1803285

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