Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 238,253 G→T 20.2% F27L (TTC→TTA U712_01105 ← Uncharacterized protein ybfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881238,2530GT20.2% 60.9 / 4.0 37F27L (TTC→TTAU712_01105Uncharacterized protein ybfE
Reads supporting (aligned to +/- strand):  ref base G (11/18);  new base T (6/2);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGACAAGTGATGAGGACAAAGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAG  >  CP006881/238123‑238382
                                                                                                                                  |                                                                                                                                 
gTAGACAAGTGATGAGGACAAAGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTAGGATTTTTTTTTGAAGAGtt                                                                                                                            >  2:159451/1‑138 (MQ=255)
gTAGACAAGTGATGAGGACAAAGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGAGTTTTCGGATTTTTGTTTGAAGCGtt                                                                                                                            >  7:72311/1‑138 (MQ=255)
gTAGACAAGTGATGAGGACAAAGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGAGGGCGTTTTCTGATTTTTGTTTGTAGCGtt                                                                                                                            >  5:66818/1‑138 (MQ=255)
        gTGATGAGGACAAAGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAAttt                                                                                                                   >  8:74176/1‑139 (MQ=255)
        gTGATGAGGACAAAGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTATGATTTTTTTTTTAATCCTTCACAAAttt                                                                                                                   >  4:130655/1‑139 (MQ=255)
           atgaGGACAAAGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTGGGATTTTTGTTTGAATCGTTCAGAAATTTGa                                                                                                                 >  1:248872/1‑138 (MQ=255)
                   aaaGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCtt                                                                                                        >  7:224522/1‑139 (MQ=255)
                          aaGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTGCGATTTTTGTTTTAATCGTTCAGAAATATTAACTTTTCTTTAAAAtt                                                                                                 >  8:154079/1‑138 (MQ=255)
                          aaGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGAGTTTTCGGATTTTTTTTTTTAGCGTACAGAAATTTTAACATTTCATGaaaaag                                                                                                 >  4:238414/1‑137 (MQ=255)
                           aGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTGTTTGGATTTTTTTTTTTAGGGGTGTGAAAATTTTACTTTTGTTGTAAATgg                                                                                                >  7:105837/1‑139 (MQ=255)
                                   tGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGaaa                                                                                           <  8:72311/136‑1 (MQ=255)
                                       ttCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTGTTGGGATTTTTTTTTTAAGAGTTAAGAAATTTTAAATTTTATTTTAAATTGGGAAACAaaaa                                                                                     >  2:205947/1‑138 (MQ=255)
                                        tCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTGGGGTTTTTTTTTGTTTGGTTGAGAAATTTTTTATTTTATTTAAAATTGGGAAAAAATaaa                                                                                    >  2:269217/1‑138 (MQ=255)
                                          gagaTTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAAc                                                                                 <  1:205947/139‑1 (MQ=255)
                                                    agagATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAgg                                                                       <  2:289525/139‑1 (MQ=255)
                                                           cTTTTTTTGTATTTTTGTTCTAATACGAAGTCGCTCATACCGACTATGCGGGCGTGGTGGGTTGTTTGTGTGTATCGGTCAGAAATTTGAACGTTTGTTTCTAATGGGGGAACAATAAAAACCCAAAAAAAGCCGagg                                                                 >  8:54688/1‑136 (MQ=255)
                                                                tttGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTGTTGGGATTTTTGTTTGAAGCGTTCAGAAATTTTAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGc                                                            >  5:211662/1‑138 (MQ=255)
                                                                    tATTTTTGTTCTAATACGACGTCGATCATACAGACTATGAGGGCGGTTTCTTATTTTTTTTTTAAGAGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAATAAACCAATAAAGGCCGATCAAGGCaaaa                                                        <  7:274325/138‑1 (MQ=255)
                                                                    tATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCaaaa                                                        <  5:27840/138‑1 (MQ=255)
                                                                            ttCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGt                                               <  2:21166/139‑1 (MQ=255)
                                                                            ttCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGt                                               <  4:43885/139‑1 (MQ=255)
                                                                                    cGAAGTCGATCATACCTAGTATGCGGGCGTTTTGGGATTTTTGTTTTAAGGTTTAATAAATTTGAACTTTTATTTAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTaa                                       <  2:269906/139‑1 (MQ=255)
                                                                                    cGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTaa                                       <  8:173126/139‑1 (MQ=255)
                                                                                    cGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTa                                        <  5:106129/138‑1 (MQ=255)
                                                                                       aGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTc                                    <  7:74176/139‑1 (MQ=255)
                                                                                            aTCATACCGACTATGCGGGCGTGTTCGGTATTTTTTTTTAAGCGGTCAGAAATTTGAAGTTTTGTTTTTAATGGTGAAAAAAAAAAAAACCCAAAAATACAGAAGAAAGCAAAAAACATTTGTTCAGCTAATTATAtt                                >  6:90615/1‑138 (MQ=255)
                                                                                             tCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTaa                              <  3:112931/139‑1 (MQ=255)
                                                                                               aTACCGACTATGCGGGCTTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAAt                             <  6:66818/138‑1 (MQ=255)
                                                                                                 aCCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGGTCAGAAATTTTAACTTTTGTTGTAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTTTACAGCTAATTATATTTATTCa                          >  5:104364/1‑139 (MQ=255)
                                                                                                                 ttttCGGATTTTTTTTTGAAGCCTTAAGAAATTTTAACTTTTATTGAAAATGGGGAAACAATAAAAAACCAAAAAAGGCCCATAAAGGCAAAAAAAAATTTTACAGCTAATTCTATTAATTAATTAAAAAGAGTGCTg           >  3:94103/1‑138 (MQ=255)
                                                                                                                 ttttCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGa          <  3:130655/139‑1 (MQ=255)
                                                                                                                      ggATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGaaaaaa     <  2:248872/139‑1 (MQ=255)
                                                                                                                         tttttGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAg  <  1:159451/139‑1 (MQ=255)
                                                                                                                         tttttGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAg  <  3:200567/139‑1 (MQ=255)
                                                                                                                          ttttttttGATGCGTTCATAAATTTGAAATTTTCTTTAAATTGGGGAAACAAAAAAAACCAAAAAAAGGCCGATAAAGGCAAAAAACATTTGTAAAGCTAATTCTATAAATTCATTAAAACGAGTGCTGAAAAAATAg  <  2:140845/138‑1 (MQ=255)
                                                                                                                          ttttGTTTGAAGCGTTCATAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAg  <  5:46258/138‑1 (MQ=255)
                                                                                                                          ttttGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAg  <  7:154079/138‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                 
GTAGACAAGTGATGAGGACAAAGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAG  >  CP006881/238123‑238382

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: