Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,786,467 T→C I193V (ATT→GTT)  yibB ← protein YibB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,786,4670TC90.5% 58.2 / ‑0.2 21I193V (ATT→GTT) yibBprotein YibB
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (13/6);  total (14/7)
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.52e-01

TTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGG  >  NZ_CP009273/3786384‑3786551
                                                                                   |                                                                                     
ttCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGaa                                                                                 >  1:354392/1‑90 (MQ=255)
       aCTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCaaaa                                                                          <  2:276169/90‑1 (MQ=255)
           aTTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAATTCAAAAACCt                                                                      <  2:323290/90‑1 (MQ=255)
                   aaTTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTga                                                              <  1:81099/90‑1 (MQ=255)
                                     tgtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCt                                           <  2:5346/90‑1 (MQ=255)
                                      gtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCtt                                          >  2:346061/1‑90 (MQ=255)
                                                   tGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGaa                             <  2:262194/90‑1 (MQ=255)
                                                                   tataCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACtt             >  1:139696/1‑90 (MQ=255)
                                                                   tataCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACtt             >  2:272146/1‑90 (MQ=255)
                                                                   tataCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACtt             >  1:2912/1‑90 (MQ=255)
                                                                      aCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtca          >  2:363522/1‑90 (MQ=255)
                                                                       cATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcat         >  1:62701/1‑90 (MQ=255)
                                                                         tGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcatca       >  2:291686/1‑90 (MQ=255)
                                                                          gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAg                                                <  2:14542/49‑1 (MQ=255)
                                                                          gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAg                                                >  1:14542/1‑49 (MQ=255)
                                                                          gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTTATAGTGAAAAGATGCATCTTACTTtcatca       >  2:209464/1‑90 (MQ=255)
                                                                           aTTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaaa     >  1:38952/1‑90 (MQ=255)
                                                                            ttGCCAATCATGAAATCAAAAACCTGTTGTTgaga                                                            <  1:330081/35‑1 (MQ=255)
                                                                            ttGCCAATCATGAAATCAAAAACCTGTTGTTgaga                                                            >  2:330081/1‑35 (MQ=255)
                                                                            ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt    >  1:107300/1‑90 (MQ=255)
                                                                              gCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATgg  >  1:302035/1‑90 (MQ=255)
                                                                                   |                                                                                     
TTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGG  >  NZ_CP009273/3786384‑3786551

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