Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 790,348 T→G 19.2% intergenic (‑103/‑141) U712_03825 ← / → U712_03830 Uncharacterized protein yfmB/General stress protein 17M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881790,3480TG19.2% 82.0 / 2.6 48intergenic (‑103/‑141)U712_03825/U712_03830Uncharacterized protein yfmB/General stress protein 17M
Reads supporting (aligned to +/- strand):  ref base T (19/18);  new base G (8/2);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTTCTGGTGAAAAGTATTGCATTAGACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGTTTCAGGTACAGACGATCGGGTATGAAAGAAATATAGAAAACATGAAGGAGGAATAT  >  CP006881/790222‑790484
                                                                                                                              |                                                                                                                                        
cTGTTCTGGTGAAAAGTATTGCATTAGACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCGCCTTAATTTTTATTTTTGTGAAAATAGGTCCCGTTCCTATAAcgt                                                                                                                               >  6:184534/1‑136 (MQ=255)
cTGTTCTGGTGAAAAGTATTGCATTAGACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTACGTGCAGCGGCTTGTACTTTATTTTGATTTTTATGAAAATATGTCCCGGTACATTCACttt                                                                                                                              >  5:106377/1‑139 (MQ=255)
   ttCTGGTGAAAAGTATTGCATTAGACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGGGCTCCGGGCCGCGGCGTCTCCGGAAATTTTAATTTTATGAAAAAAGGGCCCCGGCCTATCACCTTaac                                                                                                                           >  2:42653/1‑137 (MQ=255)
     cTGGTGAAAAGTATTGCATTAGACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTACCGGTCCTATCACTTTACCAt                                                                                                                         >  8:30734/1‑139 (MQ=255)
     cTGGTGAAAAGTATTGCATTAGACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCGTAATTTTGTTTTTTCTGAAAATAGGTCACGTTCCTATCACTTTAACAt                                                                                                                         >  8:308449/1‑139 (MQ=255)
        gTGAAAAGTATTGCATTAGACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGAGGCTTCGCCATTATTTTTCTTTTTTTTTAAAAATGTCCCGGTACGATCAAATTACCAAGGa                                                                                                                      >  1:5045/1‑139 (MQ=255)
             aaGTATTGCATTAGACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTGCTTAATTTTGATTTTTCTGAAAATAGGGCACGTTGCTATAACATTACCATgtccggg                                                                                                                  >  8:253235/1‑132 (MQ=255)
                  ttGCATTAGACATACCTCCTGCTCGTACGGATACAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAACATAGGTACCGTTCCTATCACTTTACCATGGACGGAAAACac                                                                                                            <  6:129327/139‑2 (MQ=255)
                  ttGCATTAGACATACCTCCAGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACaa                                                                                                            <  8:286676/139‑1 (MQ=255)
                    gCATTAGACCTCCCTCCTGCTCGTACGGAAAACCGCAGCGGTTCATGGGCGTGTGCTCCGTGCAGCGGCTTCTCCTTACTTTTGATGTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACaaa                                                                                                           <  1:117449/138‑1 (MQ=255)
                       ttAGACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGc                                                                                                       <  5:184534/139‑1 (MQ=255)
                        tAGACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTTATTTTTATGAAAATAGGTCCCGTTACTATCACTTTAACATTGACGGAAAACAAAAAGc                                                                                                       >  6:217249/1‑138 (MQ=255)
                         aGACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGc                                                                                                       <  8:37478/137‑1 (MQ=255)
                          gACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGcta                                                                                                     <  1:42653/138‑1 (MQ=255)
                          gACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTAACATGGACGGAAAACAAAAATctaa                                                                                                    >  3:133664/1‑138 (MQ=255)
                               aCCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGGTCCGTGCCGCGGCTTGTCCGTAATTTTTTTTTTTGTGAAAATAGGTCCCGTTGCGATCACTTTACCATGGGCGGAAAACAAATATCTACTAcc                                                                                                >  5:233955/1‑138 (MQ=255)
                               aCCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGGAGCGGCGTGGGCGTTATTTTTTTTTTTGTGAAAATAGGTGCCGTTCCGATAACGTTAACAAGGGCGGGAAACAAATATctgctcac                                                                                                >  2:220906/1‑135 (MQ=255)
                               aCCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTATGAAAATAGGTCCCGTTCCTATCACTTTTCCATGGACGGAAAACAAATAGCTACGACCa                                                                                               >  5:123102/1‑139 (MQ=255)
                               aCCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGGGCAGCGGCTTCGCCGTAAATTTGATTTTTATGAAAATAGGTCCCGGTGCTATCACATTTACATGGCCGGAAAAAAAAAAACtaataca                                                                                                >  4:99195/1‑137 (MQ=255)
                               aCCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCGCCGTGGAGCGGCTTGTCCAATATTTTGTATTTTGTGAAAAAAAGTACCCGTGCGATCACTTTTACCTGGACGGAAAAAAAAAATCTAATAcc                                                                                                >  5:59431/1‑138 (MQ=255)
                                       cTCGTACGGATAAAGGCAGCGGTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATtcctc                                                                                         <  7:253235/137‑1 (MQ=255)
                                             cGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGttttt                                                                                 <  5:217249/139‑1 (MQ=255)
                                                      gCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAAt                                                                        <  2:88858/139‑1 (MQ=255)
                                                          cGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCGATCCCTTTAACATGGTCGGGAAACAAATAGCTACTACCCTTGCTCCTGTTTTTGTGTTAAATGTTc                                                                    >  7:6827/1‑139 (MQ=255)
                                                          cGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCGTTATTTTTATTTTTCTGAAAATAGGTCCCGGTCCTTTAACTTTAACATGGTCGGAAAAAAAATAGCTACTAACAATCCGCCGGGTTTTGTGTTTAAAGttg                                                                    >  3:188531/1‑138 (MQ=255)
                                                                aTGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGaa                                                               <  4:188531/138‑1 (MQ=255)
                                                                aTGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAAt                                                              <  2:169087/139‑1 (MQ=255)
                                                                aTGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATGTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAAt                                                              <  2:5045/139‑1 (MQ=255)
                                                                      gtgtGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAACTAGGTCCCGGTCCTATCACTTTCCCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGttt                                                        <  5:64182/139‑1 (MQ=255)
                                                                      gtgtGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGttt                                                        <  1:56713/139‑1 (MQ=255)
                                                                      gtgtGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCACGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTACTGTTTTTCTCTTCAATGTTCTGGAATCTGttt                                                        <  5:170988/139‑1 (MQ=255)
                                                                              cGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTATCTTCAATGTTCTGGAATCTGTTTCAGGTa                                                  >  7:252106/1‑137 (MQ=255)
                                                                                tGCAGCGGCTTCTGCTTAATTTTGATTTTTATGAAAATAGGTCCCGGTCCTATCACTTTACCATGGACGGAAAAAAAAAAGCTACGACCAATACGCCTGTTTTTATCTTCAATGTTCTGGAATCTGTTTCAGGTAAAGa                                              >  8:48135/1‑139 (MQ=255)
                                                                                tGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGTTTCAGGTACAGa                                              >  3:305845/1‑139 (MQ=255)
                                                                                tGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCAATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTACTCCTGTTTTTCTCTTCAATGTTATGGAATCTGTTTCATGTTCAg                                               >  3:245275/1‑138 (MQ=255)
                                                                                     cGGCTTCTCCTTAATTTTTATTTTTCTAAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGTTTCAGGTACAGACGATc                                         >  1:76808/1‑139 (MQ=255)
                                                                                          tctcCGTAATTTTTATTTTTATGAAAATATGTCCCGTTCCGATAACGTTTCCATGGGCGGAAAACAAATAGCTACATCCATTCCTGCTGTTTTTGTGTTAAATGTTGTGGAATATGTTTGAGGTTAAGGCGGTAGGGTa                                    >  8:17901/1‑139 (MQ=255)
                                                                                                 aaTTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGATTCAGGTACAGACTATCGGGTATgaaag                              >  1:149995/1‑138 (MQ=255)
                                                                                                 aaTTTTGATTTTTCTGAAAATAGGTCCCGGTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGTTTCAGGTACAGACGATCGGGTATgaaa                               <  6:123102/137‑1 (MQ=255)
                                                                                                      tGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGTTTCAGGTACAGACGATCGGGTATGAAAGAAata                         >  3:91211/1‑138 (MQ=255)
                                                                                                      tGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGTTTCAGGTACAGACGATCGGGTATGAAAGAAata                         <  7:48135/138‑1 (MQ=255)
                                                                                                        aTTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGTTTCAGGTACAGACGATCGGGTATGAAAGAAata                         >  6:246708/1‑136 (MQ=255)
                                                                                                        aTTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCATCCTGTTTTTCTCTTCAATGTACTGGAATCTGTTTCAGGTACAGACGATCGGGTATGAAAGAAata                         >  4:257708/1‑136 (MQ=255)
                                                                                                                  aaaTAGGTCCCGTTCCTATCTCTTTACCATGTAAGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGTTTCAGGTACAGACGATCGGGTATGAAAGAAATATAGAAAACATGa             <  3:12387/138‑1 (MQ=255)
                                                                                                                    aTAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGTTTCAGGTACAGACGATCGGGTATGAAAGAAATATAGAAAACATGAagg          <  7:308449/139‑1 (MQ=255)
                                                                                                                    aTAGGTCCCGATCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCATTGTTCTGGATTCTGTTTGAGGTACAGACGATCGGGTATGAAATAAATATAGAAAACTTGAagg          <  2:226801/139‑1 (MQ=255)
                                                                                                                           ccGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGTTTCAGGTACAGACGATCGGGTATGAAAGAAATATAGAAAACATGAAGGAGGAata   >  4:296329/1‑139 (MQ=255)
                                                                                                                            cGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGTTTCAGGTACAGACGATCGGGTATGAAAGAAATATAGAAAACATGAAGGAGGAatat  <  2:76808/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                        
CTGTTCTGGTGAAAAGTATTGCATTAGACATACCTCCTGCTCGTACGGATAAAGGCAGCGTTTCATGGTCGTGTGCTCCGTGCAGCGGCTTCTCCTTAATTTTGATTTTTCTGAAAATAGGTCCCGTTCCTATCACTTTACCATGGACGGAAAACAAATAGCTACTACCATTCCTCCTGTTTTTCTCTTCAATGTTCTGGAATCTGTTTCAGGTACAGACGATCGGGTATGAAAGAAATATAGAAAACATGAAGGAGGAATAT  >  CP006881/790222‑790484

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: