Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,143,377 C→T S283L (TCA→TTA)  frwC → putative PTS enzyme IIC component FrwC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,143,3770CT93.9% 96.7 / ‑0.8 33S283L (TCA→TTA) frwCputative PTS enzyme IIC component FrwC
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (12/19);  total (13/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-01

TGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTT  >  NC_000913/4143148‑4143597
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                            
tGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTgc                                                                                                                                                                                                         >  1:295091/1‑251 (MQ=255)
  cGTTGCTCAGGGAGTTTATACCGTGGTGGCCATCACTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGCCGCTGCTAGGTCGAAAAAATTTCTCCGCAGAAGAGAGCGCAACAGGCCAAGCGGCACTGGTAATG‑GGATGTGTCGGCGTTACTGAAGGGGCGATTACTTTCGCCGCTGCCCATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGATTAGTTTGTGGTGCAGTAACTgcgg                                                                                                                                                                                                       >  1:319425/1‑251 (MQ=255)
                         ggggggCCATCGCGGCTGGGGGTGTGTGTAGCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTTCTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTAt                                                                                                                                                                                <  2:163831/247‑1 (MQ=255)
                               ccatggcggcggggggtagtggTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGAGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCGGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGc                                                                                                                                                                          <  2:326336/230‑1 (MQ=255)
                                   cGCTGCCTTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCAGCATGGTCGGTTTAGTGTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTggg                                                                                                                                                                      <  1:165919/251‑1 (MQ=255)
                                   cGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGCTG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTACTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGTTGCGCAGTGCTATGCAGGCTggg                                                                                                                                                                      >  1:286248/1‑251 (MQ=255)
                                          gTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCt                                                                                                                                                               >  1:117909/1‑251 (MQ=255)
                                                                  gggATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGGCTGATTGTGCTGCCGGTGGTTGAAgg                                                                                                                                       >  1:45505/1‑251 (MQ=255)
                                                                                  cGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAACCCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGACGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTCGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCCGTGCTATGCAGGCTGGGGTGGTCTGATTGTTCTGCCGGGGGTTGGAGGCCAGCTGGGTTTTATc                                                                                                                       >  2:23787/1‑251 (MQ=255)
                                                                                                     gctctcCGCATAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGTCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAg                                                                                                    <  1:236459/249‑1 (MQ=255)
                                                                                                      ttccccGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAgg                                                                                                                                       <  2:303741/211‑1 (MQ=255)
                                                                                                      ttCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAgg                                                                                                                                       >  1:303741/1‑215 (MQ=255)
                                                                                                          ccGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCa                                                                                                  <  1:215551/248‑1 (MQ=255)
                                                                                                          ccGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAgt                                                                                                <  1:172960/250‑1 (MQ=255)
                                                                                                          ccGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCa                                                                                                  <  2:94467/248‑1 (MQ=255)
                                                                                                                aagaGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTAGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCAGCTGCGTGTCATTCCTTCCATCATGGACGGTTTAGATTGTGGTGCAGTAACTGCGGCGCTGGACGGTGCGCAGTGCCATGAAGGCTGGGGTGGACTGATTGTGCTGCCGGTGGGTGAAGGCAAGCTGGGTTATAGCGCAGCAGTGGCTGGCGGAGCAGTGGTGACg                                                                                         >  1:200689/1‑251 (MQ=255)
                                                                                                                   ggcggggaaCCGGCAAAGCGGACCTGGTGATG‑GGCTGTTTTGGCGTGACTGAAGTGGCGATTCCTTTCGCCTCTGCCGATCCGCTGCGTGGCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCTGTGCGCAGTGCTATGCAGGTGGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCt                                                                                      <  2:289585/244‑1 (MQ=255)
                                                                                                                       ggaagccgggtaaGCGGCACTGTTG‑TGTGGCTGTGTGGGCGTGACTGAAGTGGCGATTCCTTTCGCCGCTGTCGATCCGCGTCGTGTCATTCCTTCCAGCATGGTCGGTTTAGTTTGTGGTGCAGGAGCTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGTTGGGCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCGGCAGTGGCTGTCGGAGAAGTGGTGACGGCTGttt                                                                                  <  2:269701/240‑1 (MQ=255)
                                                                                                                         aaaCCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGt                                                                                    >  1:9488/1‑247 (MQ=255)
                                                                                                                         aaaCCGGCAAAGCGGCACTGGTGATG‑GGCGGTGTGGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGt                                                                                    <  2:9488/247‑1 (MQ=255)
                                                                                                                                                     gCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGcgcgg                                                     >  1:121207/1‑250 (MQ=255)
                                                                                                                                                                   aCTGAAGGTGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCtt                                        <  1:93625/250‑1 (MQ=255)
                                                                                                                                                                      acaggggcttttCCTTTCTCCGCTGACGATCCTCTGCGTGTCAGTCCTGCCAGCATTGTCGGTTTAGGTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTTCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTc                                       <  2:304003/238‑1 (MQ=255)
                                                                                                                                                                              gATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGaa                             >  1:75355/1‑250 (MQ=255)
                                                                                                                                                                                  ccTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCGGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAaag                         <  1:23787/250‑1 (MQ=255)
                                                                                                                                                                                     ttCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCggt                                                                                                                                                <  2:165169/128‑1 (MQ=255)
                                                                                                                                                                                     ttCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCggt                                                                                                                                                >  1:165169/1‑128 (MQ=255)
                                                                                                                                                                                      tCGCCGCTGCCGATCCGCTGCGTGTCATTCCGTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGGCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAagaaga                     <  2:106658/250‑1 (MQ=255)
                                                                                                                                                                                       cgccgcTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAgacg                   <  1:318349/251‑1 (MQ=255)
                                                                                                                                                                                              gCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGc                                                                      <  1:300373/193‑1 (MQ=255)
                                                                                                                                                                                                  aTCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCtggatttg        >  1:172703/1‑251 (MQ=255)
                                                                                                                                                                                                   tCCGCTGCGTGTCATTCCTTCCATCAGGGGCGGTTTAGGTTGTGGTGCAGTAACTGCGGCGCTGGGCAGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCtggatttg        <  2:284172/250‑1 (MQ=255)
                                                                                                                                                                                                         gtgtgttattccgctcatactGGTCGGTTTGGTGTGTGGGGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGGCTGAGTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGTATTTGGAtttt  <  2:275187/230‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                            
TGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATG‑GGCTGTGTTGGCGTTACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTT  >  NC_000913/4143148‑4143597

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