Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,123,224 A→T 32.2% intergenic (+39/‑71) yceM → / → murJ putative oxidoreductase/putative peptidoglycan lipid II flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,123,2240AT32.2% 30.5 / 20.2 31intergenic (+39/‑71)yceM/murJputative oxidoreductase/putative peptidoglycan lipid II flippase
Reads supporting (aligned to +/- strand):  ref base A (12/9);  new base T (4/6);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

TTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTT  >  CP009273/1123082‑1123359
                                                                                                                                              |                                                                                                                                       
ttGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCTGGTAGa                                                                                                                                   <  2:295127/149‑1 (MQ=255)
ttGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGa                                                                                                                                   <  1:217204/149‑1 (MQ=255)
         tgCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGCCAAGATCTGGCGCGATGCGATCAGTGCCCCATCGCGGTATCATCGGGCGGTAGCCATTCACCATAATCCTGGTAGACTATTCGcc                                                                                                                          >  1:76325/1‑149 (MQ=255)
          gCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGACCTGGCGCGATGCGATGAGTGAATAAACCATGTAACATCTGGCGGGAGCAATTCACCATAATCCTGGTAGACTAGTCGcca                                                                                                                         >  2:156425/1‑148 (MQ=255)
                       gTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCTGGTAGACTATTCGCCCCTTTCAGCgcctg                                                                                                            >  2:101216/1‑149 (MQ=255)
                            gCAAGCCGCCGGCGAACAGGCCGTGCGGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTACCATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTg                                                                                                       <  2:140621/149‑1 (MQ=255)
                             cAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGc                                                                                                      >  1:249033/1‑149 (MQ=255)
                             cAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATAGCTAAGCAGATCTGGCGCGATGCGCTGAGTGAATCAACCCTGTAACATCTGGCGGAAGCAATTCACCATAATCCAGGTAGCCTATTCGGCTCTTTCAGCGCCTGCCTTGc                                                                                                      >  1:93621/1‑149 (MQ=255)
                                         gAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGTGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGcc                                                                                          <  1:290182/149‑1 (MQ=255)
                                            cAGGCCGTGCTGGCGGAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCGGTAACATCTGGCGGTAGCAATTCACCATAATCAAGGTAGACTATTCGCCTCGTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGt                                                                                       >  1:225701/1‑149 (MQ=255)
                                               gCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCTGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTggg                                                                                    >  1:351057/1‑149 (MQ=255)
                                                         cgcAACGTACCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCTGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAg                                                                          <  2:73368/149‑1 (MQ=255)
                                                            aaCGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAac                                                                       >  2:193005/1‑149 (MQ=255)
                                                                 atccgtgactacagcTGGCGCAATGTGAGGAGTGAATAAACCCTGCAACATCTGGCGGTAGCATTTCACCATAAGCCTGGTAGACTGTGCGCCTCTTTCAGCGCCTGCCTTGCTGGCGGTTTGCCCGTGGGTCTGGATTAGAACACCGa                                                                  <  1:71946/135‑1 (MQ=255)
                                                                         cAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCACTTCACCATAATCCAGGTCGACTATTCGCCTCTTTCAGCGCCTGCCTGGCAGGCGTTTTGCCCGTGGGTCGGGATTAGAACACCGATGCGTtta                                                          >  2:89973/1‑149 (MQ=255)
                                                                         cAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATtta                                                          >  1:36978/1‑149 (MQ=255)
                                                                         cAAGATCTGGCGCGATGCGATGAGGGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATtta                                                          <  1:152383/149‑1 (MQ=255)
                                                                         cAAGAGCTGGAGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACCCCGATTAATtta                                                          <  1:151151/149‑1 (MQ=255)
                                                                             gtCTGGCGCGATGCGATGAGTGTATCAACCCCGTAACATCTGGCGGTAGCAATTCACCATAATCCTGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTaa                                                      <  1:234087/148‑1 (MQ=255)
                                                                             aTCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCTGGTAGACTATTCGCCTCTTTAAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTaa                                                      <  2:140638/149‑1 (MQ=255)
                                                                                  ccgcgATGCGATGCGTCAATCGGCCCGAGAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCg                                                 <  1:250435/148‑1 (MQ=255)
                                                                                      gatgcgatgAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACCCCGATGAATTTATTAAAATCGCTgg                                             >  2:316111/1‑149 (MQ=255)
                                                                                       atgcgatgAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATACAGGTACAGTATTGGCCTCTGTCGGGGCCTGCCTTGGAGGCGTTTTGCCCGGTGGTCTGGATTAGAACCCCGCTGAAGTTATTAAGTTCGCTGgc                                            >  1:234500/1‑149 (MQ=255)
                                                                                        tgtgatgAGTGCATAAACCCTTTAACCTCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGgcc                                           <  1:179520/146‑1 (MQ=255)
                                                                                                 tGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTGGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTGTGGCCCGTGGGTCTGGATTAGAACACCGATGAGTTTATTAAAATCGCTGGCCGCCGTCAGc                                  >  2:305351/1‑149 (MQ=255)
                                                                                                 tGAATAAACCCTGTAACATCTGGCGGTAGCAATGCACCATAATCCAGGTGGCCTATTAGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGACTAGACCAGCGATGAATTTATTAAAATCTCCGGGCGCCCTCAGc                                  >  1:315359/1‑149 (MQ=255)
                                                                                                   aaTAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTc                                >  1:37232/1‑149 (MQ=255)
                                                                                                          ccTGTATCATCTGTCGGTAGCAATTCACCATAATCCAGGTAGACTATTTGCCGCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACTCCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGAcc                         <  1:252043/149‑1 (MQ=255)
                                                                                                           cTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGCCCa                        >  2:65155/1‑149 (MQ=255)
                                                                                                                              cAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTgg     <  2:37232/149‑1 (MQ=255)
                                                                                                                                 ttCACCATAATCCTGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTGTATCCGTTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCtt  <  2:351057/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                       
TTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTT  >  CP009273/1123082‑1123359

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