Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,183,476 C→A 20.6% G305V (GGA→GTA)  yehB ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,183,4760CA20.6% 60.0 / 2.2 24G305V (GGA→GTA) yehBputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (9/10);  new base A (3/2);  total (12/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 = 9.78e-01

TTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAACCTCTTTCTGATAAACCACAAAACCATTCTGTTCAATAGTTACCAGCGCGTTACTCTGAGCAATCCCCTGCACCCGTGGCGTAAAATTTTGTTTCGAGTTAGGCAACATCTGCATATCACGAAACAACCGCACACCTCTGAAGCGAACAGAATCAAAAATATCGCTTGATG  >  CP009273/2183284‑2183695
                                                                                                                                                                                                |                                                                                                                                                                                                                           
ttattaAAACCATACTGATAACCAGCATGGACAAAATCACTTTGTTTGCTCGCCCCATCAATATGGCTACGACCAGACGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCGTTCACGCGGACATCAAGATCTGCTACACCACCGGCCAACTGCAAAACTGCAATCGAAAACGAACAAGGAGGAACCGATAGTCGAGAAACCACAAAAANATTCTgccc                                                                                                                                           >  2:262509/1‑272 (MQ=255)
                 aTAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCaccacc                                                                                                                                                                                                                      >  2:135890/1‑183 (MQ=255)
                     cccGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGcc                                                                                                                                                                                                                  <  1:135891/183‑1 (MQ=255)
                                                                  cTACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATcg                                                                                                                                                                                               >  1:397211/1‑157 (MQ=255)
                                                                      gACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGaa                                                                                                                                                                                           <  2:397206/157‑1 (MQ=255)
                                                                        cccGCCGCTAAATGATATTTCGACACGCCGGGTTGCAGAATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTGAACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAACCTCTTTCTGATAAACCACAAAACCATTCTGTTCAATAGTTACCAGCGCGTTACTCTGAGCAATCCCCTGCACCCGTGGCGTAAAATTTTGTTTCGAGTTAGGCaa                                                                   <  1:383718/275‑1 (MQ=255)
                                                                                       tatTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATcg                                                                                                                                                                                               >  2:177152/1‑136 (MQ=255)
                                                                                           tCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGaa                                                                                                                                                                                           <  1:177156/136‑1 (MQ=255)
                                                                                                   acGGGTTGCANAATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAACATCTTTCTGATAAACCACAAAACCATTCTGTTCAATAGTTACCAGCGCGTTACTCTGAGCAATCCCCTGCACCCGTGGCGTAAAATTTTGTTTCGAGTTAGGCAACATCTGCATATCACGAAACAACCGca                                         <  2:52141/273‑1 (MQ=255)
                                                                                                            agcaTATTTGGCACCGCTGCATAATGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAACCTCTTTCTGATAAACCACAAAACCATTCTGTTCAATAGTTACCAGCGCGTTACTCTGAGCAATCCCCTGCACCCGTGGCGTAAAATTTTGTTTCGAGTTAGGCAACATCTGCATATCACGAaacaac                                             >  1:52141/1‑261 (MQ=255)
                                                                                                                                      cACCAGATGGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGCGCCAGGAGGAACCTCTTTCTGATAAACCACAAAACCATTCTGTTCAATAGTTACCAGCGCGTTACTCTGAGCAATCCCCTGCACCCGTGGCGTAAAATTTTGTTTCGATTTAGGCAACATCTGCATATCACGAAACAACCGCACACCTCTGAAGCGAACAGAATCAAAAATATCGCTTg     <  1:262513/275‑1 (MQ=255)
                                                                                                                                       aCCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTACACCACCGGCCAACTGCAAATCTGTAATcg                                                                                                                                                                                               >  1:483465/1‑88 (MQ=255)
                                                                                                                                         cAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAACCTCTTTCTGATAAACCACAAAACCATTCTGTTCAATAGTTACCAGCGCGTTACTCTGAGCAATCCCCTGCACCCGTGGCGTAAAATTTTGTTTCGAGTTAGGCAACATCTGCATATCACGAAACAACAGCAAACCTCTGAAGCGAACAGAACCAAAAAAAACGCTTAAtg  >  1:182356/1‑275 (MQ=255)
                                                                                                                                           gctAGGTGGTTACCGCGCCGTCCGCCTCTTTCACGCTGACATCACGATCTGCTACACCCCCGGCCAACTGCAAATCTGTAATCGCGaa                                                                                                                                                                                           <  2:483457/86‑1 (MQ=255)
                                                                                                                                              aGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAACCTCTTTCTGATAAACCACAAAACCATTCTGTTCAATAGTTACCAg                                                                                                                               >  1:232206/1‑145 (MQ=255)
                                                                                                                                                gtggtTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAAc                                                                                                                                                                            >  2:245736/1‑98 (MQ=255)
                                                                                                                                                  ggtTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAACCTCTTTCTGATAAACCACAAAACCATTCTGTTCAATAGTTACCAGcgcg                                                                                                                           <  2:232202/145‑1 (MQ=255)
                                                                                                                                                    ttACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAACctct                                                                                                                                                                        <  1:245740/98‑1 (MQ=255)
                                                                                                                                                          aGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTACACCACCAGCCAACTGCAAATCTGTAATcg                                                                                                                                                                                               >  1:470916/1‑69 (MQ=255)
                                                                                                                                                            ccGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCTAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAACCTCTTTCTGATAAACCACAAAACCATTCTGTTc                                                                                                                                           >  1:620498/1‑119 (MQ=255)
                                                                                                                                                              gTCCGCCTCTTTCACGCTGACATCAAGATCTGCTACACCACCAGCCAACTGCAAATCTGTAATCGCGaa                                                                                                                                                                                           <  2:470908/69‑1 (MQ=255)
                                                                                                                                                                ccGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAACCTCTTTCTGATAAACCACAAAACCATTCTGTTCAATa                                                                                                                                       <  2:620487/119‑1 (MQ=255)
                                                                                                                                                                              cTGACAGCAAGATCTGCTCCACCACCGGTCAACTGCAAATCTGTAATcgc                                                                                                                                                                                              >  1:503410/1‑50 (MQ=255)
                                                                                                                                                                                  cATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAAc                                                                                                                                                                                          <  2:503400/50‑1 (MQ=255)
                                                                                                                                                                                                |                                                                                                                                                                                                                           
TTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTTGCTCGCCCCTTCAATATGGCTACGACCCGCCGCTAAATCATATTTCGACACGCCGGGTTGCAGCATATTTGGCACCGCTGCATAAGGCACCAGATAGGTGGTTACCGAGCCGTCCGCCTCTTTCACGCTGACATCAAGATCTGCTCCACCACCGGCCAACTGCAAATCTGTAATCGCGAACGGGCCAGGAGGAACCTCTTTCTGATAAACCACAAAACCATTCTGTTCAATAGTTACCAGCGCGTTACTCTGAGCAATCCCCTGCACCCGTGGCGTAAAATTTTGTTTCGAGTTAGGCAACATCTGCATATCACGAAACAACCGCACACCTCTGAAGCGAACAGAATCAAAAATATCGCTTGATG  >  CP009273/2183284‑2183695

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