Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,515,586 G→T 16.5% F25L (TTC→TTA)  U712_18010 ← Uncharacterized protein ywtC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,515,5860GT16.5% 112.9 / 3.8 49F25L (TTC→TTA) U712_18010Uncharacterized protein ywtC
Reads supporting (aligned to +/- strand):  ref base G (25/14);  new base T (2/7);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 3.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAAATAAAAAGAGAGTGTCTCATCAGACACTCTCTTTTTTGAATTATTTATTGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGTTTGTCACTATGATCAATATCAAACGTCAGTTTT  >  CP006881/3515448‑3515723
                                                                                                                                          |                                                                                                                                         
ggCATATAAAAAGAGAGAGTCTCTTCAGACACTCTCTTTTTTGAATTATTTATTGGCTTTTACCGGAGCTTCCTGCTGCTGGATTTTATTCAATTCCATATCAATCTTCTGTTATTCCTGATCTGTCAGCTGATCATTg                                                                                                                                           <  6:216761/139‑1 (MQ=255)
 gCAAATAAAAAGAGAGTGTCTCATCAGACACTCTCTTTTTTGAATTATTTATTGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCTATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTg                                                                                                                                           >  7:265339/1‑138 (MQ=255)
        aaaaGAGAGTGTCTCATCAGACACTCTCTTTTTTGAATTATTTATTGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTaa                                                                                                                                    <  1:34513/138‑1 (MQ=255)
           agagagTGTCTCATCAGACACTCTCTTTTTTGAATTATTTATTGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGc                                                                                                                                <  2:113650/139‑1 (MQ=255)
               agTGTCTCATCAGACACTCACTTTTTTGAATTATTTATTGGCCTTTACAGGTTCTTCCTGCTGCTTGATTTTATTCATTTCCTTATCAAAATTCTTTTTTTCCTTATCTGTCATCTGATCAATTAATTTaaaacatga                                                                                                                             >  4:151621/1‑133 (MQ=255)
                gtgtCTCATCAGACACTCTCTTTTTTGAATTATTTATTGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAgattaa                                                                                                                            >  2:251156/1‑133 (MQ=255)
                gtgtCTCATCAGACACTCTCTTTTTTGAATTATTTATTGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACa                                                                                                                           <  1:1958/139‑1 (MQ=255)
                 tgtCTCATCAGACACTCTCTTTTTTGAATTATTTATTGTCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTACATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGATTTTAAAAGCTGACa                                                                                                                           >  1:256914/1‑138 (MQ=255)
                 tggctcATCAGACACTCTCTTTTTTGAATTATATATTGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGGCAGCTTATCATTGAATTTAAAAGCTGACa                                                                                                                           <  3:15429/135‑1 (MQ=255)
                   tctcATCAGACACTCTCTTTTTTGAATTATTTATTGGCGTTTACCGGGTCTTCCTGCTGCTGGAGTTTATTAATTTCCATATCCATCTTCTTTTTTTCCTTATCTGACAGCTGATCATTGAAATTAAAAGCAGACATa                                                                                                                         >  4:7750/1‑138 (MQ=255)
                         caGACACTCTCTTTTTTGAATTATTTATTGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATg                                                                                                                   >  2:224622/1‑138 (MQ=255)
                           gACACTCTCTTTTTTGAATTATTTATTGGCTTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATAAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGTATTTAAAAGATGACATAAACATGAac                                                                                                                >  6:259397/1‑139 (MQ=255)
                              acTCTCTTTTTTGAATTATTTATTGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACg                                                                                                             <  6:260590/139‑1 (MQ=255)
                              acTCTCTTTTTTGAATTATTTATTGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACg                                                                                                             <  5:207640/139‑1 (MQ=255)
                               ctctctTTTTTGAATTATTTATTGGCGTATACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATTTCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTTAATTTAAAAGCTGACATAAACATGAACACGa                                                                                                            >  6:69535/1‑139 (MQ=255)
                                  tctTTTTTGAATTATTTATTGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAATCTGACATAAACATGAACACGATgg                                                                                                         >  4:196445/1‑139 (MQ=255)
                                          taaTTATTTATTGGCGTTTACCGGGTCTTCCTGCTTCTGGATTTTATTCATTTCCATATCAATCGTCTTTTTTTCCTGATCTGTCACCTGATCATTGATTTTAAGAGCTGCCATTAAGTGGACAACGATGGCTACTGct                                                                                                 >  2:289332/2‑138 (MQ=255)
                                            attattTATTGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGcaa                                                                                                >  5:7044/1‑138 (MQ=255)
                                                tttattGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTTAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACaaa                                                                                            <  1:224622/138‑1 (MQ=255)
                                                       gCGTTTACCGGTTCTTCCTTCTGCTGTATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGGCAGCTGATCATTGAATTTAAAAGCTGACATAAACATTAACACGATGGCTACGGCAACAAACGGCCAGa                                                                                    <  3:195816/139‑1 (MQ=255)
                                                       gCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGc                                                                                        >  6:208512/1‑135 (MQ=255)
                                                            tACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCAATCAATTTAAAATCTCACATAAACATCAACACGATGGCTACGCCAACAAACTGCCAGATGGca                                                                               >  1:212541/1‑138 (MQ=255)
                                                                      tcCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAAtt                                                                     >  8:238546/1‑139 (MQ=255)
                                                                      tcCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAAtt                                                                     >  2:209496/1‑139 (MQ=255)
                                                                       ccggctgctgGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGAGCTGTCAGCTGATCATTGAATTAAAAAGCTGACATAAACATCAACCCGAGGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAAtt                                                                     <  4:241135/135‑1 (MQ=255)
                                                                       ccTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATa                                                                                                                         <  3:25035/86‑1 (MQ=255)
                                                                       ccTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATa                                                                                                                         >  4:25035/1‑86 (MQ=255)
                                                                       ccTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAAtt                                                                     >  2:124492/1‑138 (MQ=255)
                                                                           ctgctgGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATc                                                                >  4:267688/1‑139 (MQ=255)
                                                                           ctgcggGATTTTAATCATTTACATATCAATCTTCTGTTTTTCCTGATCTGACAGCTTATCATTTAATTTAAAAGCGAACATAACCATCAACACGATGGCTACGCCAACAAACGGCCAGATGGCTTTGACAAATTTCAt                                                                 <  5:259397/138‑1 (MQ=255)
                                                                              ctgGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATa                                                                                                                         >  2:49999/1‑79 (MQ=255)
                                                                              ctgGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATa                                                                                                                         <  1:49999/79‑1 (MQ=255)
                                                                                ggATTTTATTGATTTCCATATCAATTTTCTGTTTTTCCTGATCTGTCAGCTGTTCATTGAATTTAAAAGCTGAGATAAACATGACCAGGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTGATCTTTAt                                                           <  4:236188/139‑1 (MQ=255)
                                                                                 gATTTTATTCTTTTCAATATAATTCTTGTTTTTTTAGTGTTCTGTAAGTTTTTATTTTAATTAAAAAGCTGACATAAACTTGAACACTTTGGCTACGGCAACAAACGGCCAGATGGCTTTGAGAAATTTCATCTTTAt                                                           <  2:212541/138‑1 (MQ=255)
                                                                                    ttttttCTTTTCCATTTCATTGTTATGTTTTTCTTTATCTGTCAGCTTATCATTTAATTTAAAAGATTACATAAACATTAACACTATGGCTACGGCAACACACGGCCAGATGGCTTTGACAAATTTCATCTTTATCAc                                                        <  6:7044/138‑1 (MQ=255)
                                                                                         tcttttCCTTATAAATCTTCTTTTTTTCGTGTGGTGGGAGCTTATCTTTTAATTTAAATGCTGACATAAACATTAAAACGATGGCTACGGCAACAAACGGCAAGATGGCTTTGACAAATTTCATCTTTATCACTCCGa                                                   <  3:151621/135‑1 (MQ=255)
                                                                                             ttCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGattatt                                              >  5:129645/1‑139 (MQ=255)
                                                                                             ttCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGattatt                                              >  5:3635/1‑139 (MQ=255)
                                                                                                  tatCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGAt                                         >  6:252894/1‑139 (MQ=255)
                                                                                                   atCAATCTTCTGTTTTTATTTACATTTCAGATGATAATTTATTTAAAAAGCTTAAATAAACATCACCACGATGGCTAGGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGAt                                         <  3:7750/138‑1 (MQ=255)
                                                                                                       aTCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAg                                    >  3:127161/1‑139 (MQ=255)
                                                                                                             tGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGTTTGt                               >  5:268651/1‑138 (MQ=255)
                                                                                                                    ccTGATCTGTCATCTGATCATTGAATTTAAAAGATGACATAAACATGAACAAGATGTCTACGGCAACAAACTGCCAGATGCCTTTGACAAATATCATCTTTATAACTCCGATTATTAAGATTTAAGTTTGTCACTATGa                       >  2:78595/1‑139 (MQ=255)
                                                                                                                    ccTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGTTTGTCACTATGa                       >  7:175802/1‑139 (MQ=255)
                                                                                                                     cTGATCTGTCAGCTGATCATTTAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGTTTGTCACTATGAt                      <  6:3635/139‑1 (MQ=255)
                                                                                                                              cAGATGATCATTGATTTTAAAAGCTCAAATAAACATCAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGTTTGTCACTATGATCAATATCa              <  1:251156/138‑1 (MQ=255)
                                                                                                                                    aTCATTGAATTTAAAAGCTGACATAAACATGAACACGATTGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCTATTATTTAGATTTTAGTTTGTCACTATGATCAATATCAAACGt         >  1:85603/1‑137 (MQ=255)
                                                                                                                                    aTCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGTTTGTCACTATGATCAATATCAAACGt         <  6:268651/137‑1 (MQ=255)
                                                                                                                                    aTCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGTTTGTCACTATGATCAATATCAAACGt         >  1:268312/1‑137 (MQ=255)
                                                                                                                                    aTCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGGTTGTCCCTATGATCAATATCAACCGt         >  2:110715/1‑137 (MQ=255)
                                                                                                                                       aTTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGTTTGTCACTATGATCAATATCAAACGTCAGt     <  8:175802/138‑1 (MQ=255)
                                                                                                                                          gAATTTAAAAGCTGACATAATCATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGTTTGTCACTATGATCAATATCAAACGTCAGtttt  <  5:138543/138‑1 (MQ=255)
                                                                                                                                          gAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGTTTGTCACTATGATCAATATCAAACGTCAGtttt  <  1:209496/138‑1 (MQ=255)
                                                                                                                                          gAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGTTTGTCACTATGATCAATATCAAACGTCAGtttt  <  7:238546/138‑1 (MQ=255)
                                                                                                                                          gAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGTTTGTCACTATGATCAATATCAAACGTCAGtttt  <  1:124492/138‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
GGCAAATAAAAAGAGAGTGTCTCATCAGACACTCTCTTTTTTGAATTATTTATTGGCGTTTACCGGTTCTTCCTGCTGCTGGATTTTATTCATTTCCATATCAATCTTCTGTTTTTCCTGATCTGTCAGCTGATCATTGAATTTAAAAGCTGACATAAACATGAACACGATGGCTACGGCAACAAACGGCCAGATGGCTTTGACAAATTTCATCTTTATCACTCCGATTATTTAGATTTTAGTTTGTCACTATGATCAATATCAAACGTCAGTTTT  >  CP006881/3515448‑3515723

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: ‑