Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,786,991 T→C D18G (GAT→GGT)  yibB ← protein YibB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,786,9910TC92.3% 31.5 / ‑3.5 13D18G (GAT→GGT) yibBprotein YibB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (3/9);  total (3/10)
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.50e-01

TGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAA  >  NZ_CP009273/3786909‑3787070
                                                                                  |                                                                               
tGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctct                                                                          >  2:255282/1‑90 (MQ=255)
   ggCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTAcc                                                                       <  2:61628/90‑1 (MQ=255)
    gCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCa                                                                      <  2:495485/90‑1 (MQ=255)
     ccAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTGGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCaa                                                                     <  1:523251/90‑1 (MQ=255)
     ccAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCaa                                                                     <  2:174139/90‑1 (MQ=255)
          cGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCa                                                                <  1:112860/90‑1 (MQ=255)
                  gTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGc                                                        <  1:164451/90‑1 (MQ=255)
                   tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                       <  2:101234/90‑1 (MQ=255)
                    aaCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAg                                                      >  1:47497/1‑90 (MQ=255)
                      cTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                    <  2:54701/90‑1 (MQ=255)
                              aaCATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATtgt                                            <  2:189122/90‑1 (MQ=255)
                                               cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                           <  1:242269/90‑1 (MQ=255)
                                                                        aGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATaaaaaa  >  1:570922/1‑90 (MQ=255)
                                                                                  |                                                                               
TGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAA  >  NZ_CP009273/3786909‑3787070

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