Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,629,517 C→T V203M (GTG→ATG)  arcA ← response regulator in two‑component regulatory system with ArcB or CpxA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,629,5170CT94.6% 114.5 / ‑2.7 37V203M (GTG→ATG) arcAresponse regulator in two‑component regulatory system with ArcB or CpxA
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (15/20);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGC  >  CP009273/4629451‑4629577
                                                                  |                                                            
tgggggCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACg                                                             <  1:866830/68‑1 (MQ=255)
  gtggCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTc                                                           >  1:187126/1‑68 (MQ=255)
   tggCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCt                                                          <  1:1225661/68‑1 (MQ=255)
      cGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCTACa                                                       <  1:3244148/68‑1 (MQ=255)
      cGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCTACa                                                       <  1:1954090/68‑1 (MQ=255)
      cGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACa                                                       <  1:1226943/68‑1 (MQ=255)
      cGATGATTTCCGGCGTATCCGGCGTAGATTCGAAAAGTTTACGAATACGGCGGATCGTCATGTCTACa                                                       <  1:3735084/68‑1 (MQ=255)
             ttCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGt                                                <  1:443122/68‑1 (MQ=255)
               ccggcgtatccggcgtaGATTCGAAATGTTTAAGAATACGGCGGATCGTCATGTCTACAGTACGGTCg                                              >  1:1556597/1‑68 (MQ=255)
                 ggcgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTg                                            >  1:3880257/1‑68 (MQ=255)
                   cgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCg                                          >  1:1456707/1‑68 (MQ=255)
                   cgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCg                                          >  1:3943621/1‑68 (MQ=255)
                     tatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGt                                        <  1:1246029/68‑1 (MQ=255)
                     tatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGt                                        <  1:104825/68‑1 (MQ=255)
                     tagccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGt                                        <  1:4092760/65‑1 (MQ=255)
                        ccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTc                                     >  1:500249/1‑68 (MQ=255)
                         cggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCa                                    >  1:1474755/1‑68 (MQ=255)
                           gcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGc                                  <  1:4552275/68‑1 (MQ=255)
                             gtaGATTCGAAATGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTc                                >  1:511082/1‑68 (MQ=255)
                                    cGAAATGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCg                         <  1:3042193/68‑1 (MQ=255)
                                      aaaTGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGt                       <  1:2359319/68‑1 (MQ=255)
                                      aaaTGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGt                       <  1:1948460/68‑1 (MQ=255)
                                      aaaTGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGt                       >  1:4463544/1‑68 (MQ=255)
                                         tGTTTACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCAt                    <  1:3517991/68‑1 (MQ=255)
                                            ttACGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGATTCAGCTCACGGCCGGTCAtttt                 <  1:1348896/68‑1 (MQ=255)
                                               cGAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCGGCTGACGGCCGGTCATTttctt              >  1:946183/1‑68 (MQ=255)
                                                gAATACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTttcttc             >  1:4707365/1‑68 (MQ=255)
                                                 aaTACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTca            >  1:99981/1‑68 (MQ=255)
                                                 aaTACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTca            <  1:3252797/68‑1 (MQ=255)
                                                 aaTACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTca            <  1:2570882/68‑1 (MQ=255)
                                                 aaTACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTca            >  1:1947808/1‑68 (MQ=255)
                                                 aaTACGGCGGATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTca            >  1:1675816/1‑68 (MQ=255)
                                                      ggcggATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGt       <  1:3610296/68‑1 (MQ=255)
                                                      ggcggATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGt       <  1:1926786/68‑1 (MQ=255)
                                                         ggATCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCa    <  1:1912943/68‑1 (MQ=255)
                                                         ggATCGTCATGTCTACAGAACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCa    >  1:4131253/1‑68 (MQ=255)
                                                           aTCGTCATGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCTGCAGTTCAGc  <  1:4251290/68‑1 (MQ=255)
                                                                  |                                                            
TGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGC  >  CP009273/4629451‑4629577

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