Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,627,056 C→T 100% G100D (GGC→GAC)  gacA ← Response regulator GacA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,627,0560CT97.9% 167.1 / ‑4.3 47G100D (GGC→GAC) gacAResponse regulator GacA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (22/24);  total (22/25)
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.99e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGACCCCTGCGGCTGGAACGGTTTCAGCGCCAACTGCTGGGCGATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGG  >  ACB143/4626926‑4627182
                                                                                                                                  |                                                                                                                              
agaCCCCTGCGGCTGGAACGGTTTCAGCGCCAACTGCTGGGCGATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCa                                                                                                                           <  1:293906/136‑1 (MQ=255)
   ccccTGCGGCTGGAACGGTTTCAGCGCCAACTGCTGGGCGATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAgcg                                                                                                                        >  4:272617/1‑136 (MQ=255)
      cTGCGGCTGGAACGGTTTCAGCGCCAACTGCTGGGCGATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCg                                                                                                                     <  5:320521/136‑1 (MQ=255)
      cTGCGGCGGGAACGGTTTCAGCGCCAACTGCTGGGCGATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCg                                                                                                                     <  2:73/136‑1 (MQ=255)
          ggCTGGAACGGTTTCAGCGCCAACTGCTGGGCGATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCt                                                                                                                 >  4:114940/1‑136 (MQ=255)
           gCTGGAACGGTTTCAGCGCCAACGGCTGGGCGATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTg                                                                                                                <  7:14695/136‑1 (MQ=255)
                 aCGGTTTCAGCGCCAACTGCTGGGCGATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGcaac                                                                                                            >  3:240761/1‑134 (MQ=255)
                               aaCTGCTGGGCGATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACg                                                                                            <  7:75969/136‑1 (MQ=255)
                               aaCTGCTGGGCGATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACg                                                                                            <  7:219820/136‑1 (MQ=255)
                                  tgctgGGCGATCTGCGGGCTGATATAGCGCTGGTCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGt                                                                                         >  2:338481/1‑136 (MQ=255)
                                  tgctgGGCGATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGt                                                                                         >  7:188818/1‑136 (MQ=255)
                                    ctgGGCGATCTGTGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTc                                                                                        <  8:46563/135‑1 (MQ=255)
                                    ctgGGCGATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTc                                                                                        <  2:28055/135‑1 (MQ=255)
                                               gCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCacac                                                                            >  3:353360/1‑136 (MQ=255)
                                                        tataGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCg                                                                    <  3:112492/135‑1 (MQ=255)
                                                                tGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCCCCGCCCCTACCt                                                           >  5:262819/1‑136 (MQ=255)
                                                                tGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTAccg                                                           >  6:246745/1‑135 (MQ=255)
                                                                tGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCt                                                           >  3:255660/1‑136 (MQ=255)
                                                                tGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCt                                                           >  3:319133/1‑136 (MQ=255)
                                                                 ggCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCtt                                                          <  6:108097/136‑1 (MQ=255)
                                                                    cGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGAt                                                       >  6:260110/1‑136 (MQ=255)
                                                                     ggCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATg                                                      >  1:349026/1‑136 (MQ=255)
                                                                           ggCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAggg                                                <  8:379762/136‑1 (MQ=255)
                                                                                   gAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTgcg                                        <  1:336354/136‑1 (MQ=255)
                                                                                   gAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTgcg                                        <  8:307203/136‑1 (MQ=255)
                                                                                   gAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTgcg                                        <  5:41523/136‑1 (MQ=255)
                                                                                   gAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTgcg                                        <  3:19510/136‑1 (MQ=255)
                                                                                         ccTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGcagcag                                  >  1:364013/1‑136 (MQ=255)
                                                                                             gACCATTTCGTCAATGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTg                              <  1:395497/136‑1 (MQ=255)
                                                                                             gACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGtt                                <  2:9278/134‑1 (MQ=255)
                                                                                              aCCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACATGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGc                             <  5:183282/136‑1 (MQ=255)
                                                                                              aCCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGc                             >  4:360136/1‑136 (MQ=255)
                                                                                               ccATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCg                            <  8:261255/136‑1 (MQ=255)
                                                                                               ccATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCg                            <  8:141929/136‑1 (MQ=255)
                                                                                                cATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCgg                           >  5:354505/1‑136 (MQ=255)
                                                                                                     cGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGa                                                        >  7:109550/1‑102 (MQ=255)
                                                                                                     cGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGa                                                        <  8:109550/102‑1 (MQ=255)
                                                                                                     cGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGc                      >  6:284616/1‑136 (MQ=255)
                                                                                                     cGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGc                      >  5:331984/1‑136 (MQ=255)
                                                                                                     cGTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGc                      >  8:362384/1‑136 (MQ=255)
                                                                                                      gTCAAGGCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGcc                     >  7:16782/1‑136 (MQ=255)
                                                                                                            gCCTGCACCCTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAgg               >  6:214255/1‑136 (MQ=255)
                                                                                                                   cccTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGa        <  4:255660/136‑1 (MQ=255)
                                                                                                                   cccTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGa        <  2:69282/136‑1 (MQ=255)
                                                                                                                   cccTTGGTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGa        <  4:319133/136‑1 (MQ=255)
                                                                                                                         gTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCgg  >  6:291391/1‑136 (MQ=255)
                                                                                                                         gTCAGGTAGTCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCgg  <  2:7607/136‑1 (MQ=255)
                                                                                                                                  |                                                                                                                              
AGACCCCTGCGGCTGGAACGGTTTCAGCGCCAACTGCTGGGCGATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGG  >  ACB143/4626926‑4627182

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