Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,324,255 G→A 8.0% H34H (CAC→CAT)  sfsA ← Sugar fermentation stimulation protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,324,2550GA 8.0% 156.6 / 3.9 50H34H (CAC→CAT) sfsASugar fermentation stimulation protein A
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base A (2/2);  total (25/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 = 8.43e-01

AGGCCAACCGCCCCTGGGGGGTTTCGCTGATTTCCCAGGTCCCTGGCAGCTTGCGCTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGCCAGCTCTGCAGGCCACGGGCA  >  ACB122/5324122‑5324386
                                                                                                                                     |                                                                                                                                   
aGGCCAACCGCCCCTGGGGGGTTTCGGTGATTTCCCTGGTCCCTGGCGGCTTGGGCTTGGGGTCGTTGGAGCGGCTGAACCACCCAGGCCCACCTTCACGCATGCAGTTGTGCATGGGGCCGGTTTTCGGGAAGt                                                                                                                                    >  2:3366/1‑135 (MQ=255)
             cTGGGGGGTTTCGCTGATTTCCCAGGTCCCTGGCAGCTTGCGCTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTg                                                                                                                       >  2:248308/1‑135 (MQ=255)
               ggggggTTTCGCTGATTTCCCAGGTCCCTGGCAGCTTGCGCTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCt                                                                                                                     >  1:175631/1‑135 (MQ=255)
                 ggggTTTCGCTGATTTCCCAGGTCCCTGGCAGCTTGCGCTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTcg                                                                                                                   >  8:202105/1‑135 (MQ=255)
                 ggggTTTCGCTGATTTCCCAGGTCCCTGGCAGCTTGCGCTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTcg                                                                                                                   <  5:214056/135‑1 (MQ=255)
                 ggggTTTCGCTGATTTCCCAGGTCCCTGGCAGCTTGCGCTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTcg                                                                                                                   <  3:203083/135‑1 (MQ=255)
                      ttCGCTGATTTCCCAGGTCCCTGGCAGCTTGCGCTTGGGGTCGTGGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCt                                                                                                              <  7:47363/135‑1 (MQ=255)
                        cGCTGATTTCCCAGGTCCCTGGCAGCTTGCGCTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTgg                                                                                                            >  6:74937/1‑135 (MQ=255)
                                  ccAGGTCCCTGGCAGCTTGCGCTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGa                                                                                                   <  1:15174/134‑1 (MQ=255)
                                          cTGGCAGCTTGCGCTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCa                                                                                          <  3:126233/135‑1 (MQ=255)
                                            ggCAGCTTGCGCTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGt                                                                                                                                    >  8:34340/1‑91 (MQ=255)
                                            ggCAGCTTGCGCTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGt                                                                                                                                    <  7:34340/91‑1 (MQ=255)
                                              cAGCTTGCGCTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGaaa                                                                                      >  2:50026/1‑135 (MQ=255)
                                               aGCTTGCGCTTGGGGTCCTTGGAGCGGCTGAACCACCCCAGCCCACCTTCACGCAAGCAGATGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAg                                                                                     >  8:227067/1‑135 (MQ=255)
                                                    gcgcTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTgg                                                                                                                                  <  6:1729/85‑1 (MQ=255)
                                                    gcgcTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTgg                                                                                                                                  >  5:1729/1‑85 (MQ=255)
                                                    gcgcTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCtt                                                                                >  7:297874/1‑135 (MQ=255)
                                                      gcTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCt                                                                                                                        >  3:36282/1‑93 (MQ=255)
                                                      gcTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCt                                                                                                                        <  4:36282/93‑1 (MQ=255)
                                                        ttGGGGTCGGTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAg                                                                            <  4:8834/135‑1 (MQ=255)
                                                        gtGGGGGCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCCTGCTGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATGGGCCAGAAAGCGCTTGTAg                                                                            <  6:257042/134‑1 (MQ=255)
                                                                  tGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACa                                                                  <  7:202105/135‑1 (MQ=255)
                                                                         gCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGcc                                                           >  5:235793/1‑135 (MQ=255)
                                                                         gCTGAACCACACCTGCCCACCTTCACGCATCCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGcc                                                           >  8:102749/1‑135 (MQ=255)
                                                                          cggAACCACACCTGCCCACCTTCACGCATGCAGTGGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTGCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGccc                                                          <  2:285720/133‑1 (MQ=255)
                                                                                      tGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGaa                                              >  6:137958/1‑135 (MQ=255)
                                                                                           aCCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGgaga                                         <  2:173410/135‑1 (MQ=255)
                                                                                               tCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCGTGTAGCGACGCAACAGGCGGCCCTGCTGGAGTGAAGGAGAGaaa                                     >  4:178703/1‑135 (MQ=255)
                                                                                                cACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAAc                                    >  2:42966/1‑135 (MQ=255)
                                                                                                 aCGCATGCAGTTGAGCATGGAGGCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAAc                                    <  1:50026/134‑1 (MQ=255)
                                                                                                  cGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGc                                  >  1:138777/1‑135 (MQ=255)
                                                                                                     aTGCAGGTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATc                               <  4:248749/135‑1 (MQ=255)
                                                                                                          gTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCt                          >  4:120190/1‑135 (MQ=255)
                                                                                                          gTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCt                          >  2:275469/1‑135 (MQ=255)
                                                                                                              aGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGc                       <  1:248308/134‑1 (MQ=255)
                                                                                                              aGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGc                       <  8:155296/134‑1 (MQ=255)
                                                                                                                  tGGAGCCGGTATTCGGGCAGTGGATAGTCCTCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGATAGAAACGCATCAGCCTTGCCAGct                  <  2:29564/135‑1 (MQ=255)
                                                                                                                   ggAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTc                                                     <  8:239881/99‑1 (MQ=255)
                                                                                                                   ggAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTc                                                     >  7:239881/1‑99 (MQ=255)
                                                                                                                      gCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGCCAGctct                >  5:28180/1‑133 (MQ=255)
                                                                                                                      gCCGGTATTCGGGCAATGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGCCAGctct                >  4:208097/1‑133 (MQ=255)
                                                                                                                      gCCGGTATTCGGGCAATGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGCCAGctct                >  4:211975/1‑133 (MQ=255)
                                                                                                                        cGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGCCAGCTCTGCAg            <  7:227067/135‑1 (MQ=255)
                                                                                                                        cGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGCCAGCTCTGCAg            >  2:123624/1‑135 (MQ=255)
                                                                                                                          gTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGCCAGCTCTGCAGGc          <  6:235793/135‑1 (MQ=255)
                                                                                                                          gTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGCCAGCTCTGCAGGc          <  4:40439/135‑1 (MQ=255)
                                                                                                                          gTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGACAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTTCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGCCAGCTCTGCAGGc          <  7:140496/135‑1 (MQ=255)
                                                                                                                          gTATTCGGGCAATGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGCCAGCTCTGCAGGc          <  3:208097/135‑1 (MQ=255)
                                                                                                                          gTATTCGGGCAATGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGCCAGCTCTGCAGGc          <  3:211975/135‑1 (MQ=255)
                                                                                                                                  gCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGCCAGCTCTGCAGGCCACGGGCa  >  7:140666/1‑135 (MQ=255)
                                                                                                                                     |                                                                                                                                   
AGGCCAACCGCCCCTGGGGGGTTTCGCTGATTTCCCAGGTCCCTGGCAGCTTGCGCTTGGGGTCGTTGGAGCGGCTGAACCACACCTGCCCACCTTCACGCATGCAGTTGAGCATGGAGCCGGTATTCGGGCAGTGGATAGTCATCTGCTCGCCGCTGGCCAGTTCGATATCGGCCAGAAAGCGCTTGTAGCGACGCAACAGGCGGCCCTGCTCGAGTGAAGGAGAGAAACGCATCAGCCTTGCCAGCTCTGCAGGCCACGGGCA  >  ACB122/5324122‑5324386

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