Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 804,403 T→G 30.1% C238W (TGT→TGG)  U712_03910 → PTS system N‑acetylglucosamine‑specific EIICB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881804,4030TG30.1% 36.3 / 5.5 41C238W (TGT→TGG) U712_03910PTS system N‑acetylglucosamine‑specific EIICB component
Reads supporting (aligned to +/- strand):  ref base T (9/19);  new base G (6/6);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATATTTTCTGGTTCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCATTTATGTTTTTATCTCCGCT  >  CP006881/804265‑804537
                                                                                                                                          |                                                                                                                                        
cAATATTTTCTGGTTCCATTTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCgtgt                                                                                                                                          <  8:86264/139‑1 (MQ=255)
   tatTTTCTGGTTCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTg                                                                                                                                       >  4:96702/1‑139 (MQ=255)
   tatTTTCTGGTTCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGGCTg                                                                                                                                       >  8:124962/1‑139 (MQ=255)
       ttCTGGTTCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCt                                                                                                                                    <  8:219528/138‑1 (MQ=255)
             ttCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCgtgt                                                                                                                                          <  6:169469/126‑1 (MQ=255)
             ttCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCgtgt                                                                                                                                          >  5:169469/1‑126 (MQ=255)
             ttCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTtg                                                                                                                             >  7:22774/1‑139 (MQ=255)
              tCCAATTCGGAGCATATAACGGAGGGACCGGGGCTTTGGCGCGCGTTTTCGCGAAAGATCCGACTGCGGGCACATATAAGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTtgt                                                                                                                            <  5:186714/139‑1 (MQ=255)
               ccAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGGCTTCCGGCGGCGTGTCTGGCTATGGTtgt                                                                                                                            <  6:168265/138‑1 (MQ=255)
                                    aGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACgg                                                                                                     >  3:104671/1‑139 (MQ=255)
                                      tGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGGGGGGGTTTTGTTGTGTCTGCGAAGCCCT‑CTAAACGGa                                                                                                    >  1:106813/1‑138 (MQ=255)
                                            gggATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGCCCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGGCGGGGTAATGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCggc                                                                                             >  3:123527/1‑137 (MQ=255)
                                                    gcgcgcTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGa                                                                                     <  3:96702/138‑1 (MQ=255)
                                                        gcTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGGGGGGGGTTGGTTGTGAGTTCGACGCCGT‑CGAAAACGGAAGCG‑CCGGCGGGGattga                                                                                >  2:283929/1‑136 (MQ=255)
                                                        gcTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGGGTCTGGCGATGGGTGGGACTGCGAAGCCGT‑GTAACCGGCAAGCG‑CCCGGTGGAatgat                                                                                >  4:234885/1‑139 (MQ=255)
                                                            tttCGCTAAAGATCGGACTGCGGGCACTTAAATCACTGTATTTTACCCAAAAATCTGGTGCGGGCTTCCGGCGGGGTGGCTGGATATGTTTTTCACTCCCACGCCGA‑CAAAGCGCAACGCG‑ACGGCTGGAAtgatga                                                                              <  2:278698/137‑1 (MQ=255)
                                                                 cgAAATATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATATTGTTGTTCGGTCTTCCGGGGGCGTGGCTGGCTATGGTTGTGACTGCGAAGCCGT‑ATAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGAtt                                                                       <  3:9124/139‑1 (MQ=255)
                                                                 cgAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGAtt                                                                       <  8:79542/139‑1 (MQ=255)
                                                                 cgAAAGAGCCGACTGCGGGCACATATCTGACCGGAGTTTGCCCAATCCTGATGTTCGGTCTGCCGGCGGCGTGTCTGGCGAGGGTTGTGACTGCGAAGCCGT‑CTAAACGTAAAGCG‑ACGGCTGGAATGCTGATGGGAtt                                                                       <  2:46379/139‑1 (MQ=255)
                                                                          cGACTGCGGGCACATTTATGAGCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGGCACTGCGAAGCCGT‑CTAAACGGAAAGCG‑TCGGCTGGAATGATGATTGGATTTGCTTTaa                                                               <  5:140493/138‑1 (MQ=255)
                                                                           gACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGCCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACg                                                             <  5:309282/139‑1 (MQ=255)
                                                                              tGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATTGGTTTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCGGGAATGAGGATTGGAGTTGCTTTAACGGCg                                                          >  1:91993/1‑139 (MQ=255)
                                                                                cGGGCACATATATGACCGGATTTTTCCCAATCTTTTTGTTCGGGGTTCCGGGGGCGTGTATGGCTTGGGTTGGGACTGCGACGCCGT‑CTAACCGTATAGCG‑ACGGATGGAATGATGATTGGATTTGCTTTAACGGCGt                                                         <  1:38275/138‑1 (MQ=255)
                                                                                cGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTGTGGCTATGGTTGTGGCTGCGAAGCCCTGATAACCGAAACGCG‑ACGGGTGAAATGTTGATTGAATTTGGTATAACGGCGt                                                         >  2:205536/1‑139 (MQ=255)
                                                                                  ggCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTGTGGGGTTGGTTGTGACTGCGAAGCCGT‑CTAAACGGTAAGCG‑ACGGGGGGTATGATGATTGGATTTGGTTTTACGGCGgttt                                                      >  1:238451/1‑138 (MQ=255)
                                                                                      caTATATGACCGGATTTTTCCCAATAATGATGTGCGGTGTGCGGGGGGCGGGGCGGGGTATGGTTGGCTGTGGGAAGCGGT‑CAAAAGGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATcac                                                  <  6:47272/139‑1 (MQ=255)
                                                                                      caTATATGACCGGATTTTTCCCAAGCATGATGTTCGGTCTTCCGGCGGCGTGGCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATcac                                                  <  8:297559/139‑1 (MQ=255)
                                                                                                 ggATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACGGCGAAGCCGT‑CTAAACGGAACGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACgg                                       <  1:283929/139‑1 (MQ=255)
                                                                                                 ggATTTTTCCCAATCATGATGTTCGGGCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACgg                                       <  3:234885/139‑1 (MQ=255)
                                                                                                  gATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTGTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATTATGATTGGATTTTCTTTAACGGCGTTTATCACAGGGATTaag                                        >  8:152498/1‑135 (MQ=255)
                                                                                                  gATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACg                                        >  3:81701/1‑137 (MQ=255)
                                                                                                  gATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGGGGGGCGTTGGGTGTGTCCGCGAAGCCCG‑CGAAACCGAAAGCGCACGCCGGGAA‑GATGATTCGGTTTGCTTTTACGGCGGTGTTACCAGGGGTTACggc                                      >  6:4169/1‑138 (MQ=255)
                                                                                                        tCCCAATCATGATGTTCGGTCTTCGGGCGGCGTGTCTGGCTATGGTGGTGACTGGGATGCGGT‑CTAAAGGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGa                                 <  5:167255/138‑1 (MQ=255)
                                                                                                         cccAATCATGATGTTCGGTCTTCCGGCGGCGTGGCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGAt                                <  5:236411/138‑1 (MQ=255)
                                                                                                             aTCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGt                           <  7:152498/139‑1 (MQ=255)
                                                                                                             aTCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGt                           <  6:126975/139‑1 (MQ=255)
                                                                                                                     gtgcGGGCTTCCGGCGGCGTGTGTGGCTATGGGTGTGACGGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCAttt                   <  4:202102/136‑1 (MQ=255)
                                                                                                                      ttCGGTCTTGCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCAttt                   <  1:205536/138‑1 (MQ=255)
                                                                                                                      ttCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCAttt                   <  8:22774/138‑1 (MQ=255)
                                                                                                                      ttCGGTCTTCCGGCGGCGTGGGTGGGGATGGTTGTGACTGCGAAGCCCT‑CTAAACGGAAATCG‑ACGGCGGGGATGATGTTTTGTTTTTGTTTAACGGCGTTTTTAACAGGGATTACCGAGCCGATTGTGTTTggaattt                  >  6:263591/1‑133 (MQ=255)
                                                                                                                      ttCGGTCTTCCGGCGGCGGGGCTGGCTATGGTTGGGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCAttt                   <  4:104671/138‑1 (MQ=255)
                                                                                                                                  gcggcgTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGTG‑ACGGCTGGAATGATGATTGGATTTGCTGTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCATTTATGTTTTTAtctc      <  6:254179/139‑1 (MQ=255)
                                                                                                                                    ggcgTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCATTTATGTTTTTATCTCCg    <  8:149048/139‑1 (MQ=255)
                                                                                                                                      cgTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCATTTATGTTTTTATCTCCGCt  <  2:34967/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                        
CAATATTTTCTGGTTCCAATTCGGAGAATATAACGGAGTGACAGGGGATTTGGCGCGCTTTTTCGCGAAAGATCCGACTGCGGGCACATATATGACCGGATTTTTCCCAATCATGATGTTCGGTCTTCCGGCGGCGTGTCTGGCTATGGTTGTGACTGCGAAGCCGT‑CTAAACGGAAAGCG‑ACGGCTGGAATGATGATTGGATTTGCTTTAACGGCGTTTATCACAGGGATTACGGAGCCGATTGAGTTTGCATTTATGTTTTTATCTCCGCT  >  CP006881/804265‑804537

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