Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 948,353 T→G 23.2% L140R (CTT→CGT)  U712_04645 → putative lipoprotein yhcJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881948,3530TG23.2% 69.6 / 6.7 41L140R (CTT→CGT) U712_04645putative lipoprotein yhcJ
Reads supporting (aligned to +/- strand):  ref base T (9/22);  new base G (8/2);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 8.27e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTAATTTGGCACCGCTCGGGACGACATACATAACGCCCATGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGACACCTGAAATTAAAAGCTGTCCGGCAACAGGATGCAGTGAG  >  CP006881/948217‑948482
                                                                                                                                        |                                                                                                                                 
cTTAATTTGGCACCGCTCGGGACGACATACCTAACGCCCCTGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTa                                                                                                                                 <  1:303507/139‑1 (MQ=255)
   aaTTTGGCACCGCTCGGGACGACATACATAACGCCCATGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGTGGTGTAGTTTCTGTTCCGTATAAAGCGTTCGATTTTGGCAGGGAGCTTACg                                                                                                                               <  4:238722/138‑1 (MQ=255)
     tttGGCACCGCTCGGGACGACATACATAACGCCCATGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTg                                                                                                                            <  5:38257/139‑1 (MQ=255)
                        acatacatAACGCCCATGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTGTTTGGTGCCGCTTAACCGTGGGCTTTTGGGGGGGCGGGTTCGGTGCTGCCCGGAGCGGGGTTGc                                                                                                         >  1:257168/1‑139 (MQ=255)
                         catacataACGCCCATGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGACGTTTCTGGTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGc                                                                                                         >  8:146996/1‑138 (MQ=255)
                                   gCCCATGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTTTGTGTTGGGGTTAAAGCGGTGGATTTTTGGCGGGGGCGTTCGGTGGTCCAGGTAGCGGGTTTGGTAACGCTGaa                                                                                               >  7:116031/1‑138 (MQ=255)
                                   gCCCATGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTGTGGGTTGCGGCTTAATCGGTAGGTTTTGGCAGGGGGCGGACGGTGGTCCAGGTTGCGGGGTTGCGTACCCTGaa                                                                                               >  2:116802/1‑138 (MQ=255)
                                   gCCCATGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTGTATGTTGCGGATTAAAAGTTAGATTTTTGGAATGCGCGTTCCTTGCTCCCTGAAACCGGTTTTCTAAAGCAGaa                                                                                               >  4:3406/1‑138 (MQ=255)
                                      cATGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAAc                                                                                           <  1:131124/139‑1 (MQ=255)
                                       aTGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTGTGTGTTGCGGTTAAAGCGTTGGCTTTTGGCCGGGCGCTGTCGGTGGTCCAGGAAGCGGGGTTGGTAACGGGGCAAAACg                                                                                          >  8:33059/1‑139 (MQ=255)
                                       aTGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTGTGTGTTCCGGATAAAGCGGTGGCTTTTGGGAGGGCGCGTACGGTGGTCCCGGGAGCGGGTTTGCTAACGCTGTAAAACg                                                                                          >  7:67994/1‑139 (MQ=255)
                                       aTGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTGTGTGGTGCGGATAAAAAGTTGGATTTTGGCAGGGCCCGTTCGGTGCTACAGGAAACGGGGTTTGTAACGCTGaaaaaag                                                                                          >  2:264180/1‑137 (MQ=255)
                                        tGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTGTCTGGTGCGGAGAAACCCTTGGGTTTTGGCCGGGGGGGGAAGGGGCGGCCGGAGGCGGGTTTGGTTACGATGaaaaaag                                                                                          >  6:286836/1‑136 (MQ=255)
                                              tttATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTaa                                                                                   <  8:231292/139‑1 (MQ=255)
                                                tATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAAc                                                                                  <  5:340343/138‑1 (MQ=255)
                                                     aaaaCGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAAc                                                                             <  4:384350/138‑1 (MQ=255)
                                                         cGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCGTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGtt                                                                        <  5:266743/139‑1 (MQ=255)
                                                            tatGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTGTGTGTTGCGGGTAAAGCGTTGGATTTTTGGAGGGGGCGTACGGGGGGCCCGGTAGCGGGGTTTGTAACGCGGAAAAACGGCTTTTACGGTACAgggtgag                                                                     >  6:70017/1‑134 (MQ=255)
                                                                 aCGAATCCGAGATATATCGCGGGGAGCTGTCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGat                                                                <  1:264180/139‑1 (MQ=255)
                                                                  cGAATCCGAGATATATCGCGGGGAGCTGTCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGata                                                               <  2:252797/139‑1 (MQ=255)
                                                                        cGAGATATATCGCGGGGAGCTGTCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAATTGATATGAtt                                                          <  8:209387/138‑1 (MQ=255)
                                                                               tatCGCGGGGAGCTGTCGTTGGGGTTCCGGGTTAATAGGTAGATTTTTGGCCGGCGGTTAGCGTGGGCCAGGCGGCGGGGTTTGTAACGCTTAAAAAACGATTTAAAGGCACACGTTAAGTTTTTATGTTTTAAGATaa                                                  >  6:309220/1‑139 (MQ=255)
                                                                                      gggAGCTGTCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCgaga                                           <  5:309220/139‑1 (MQ=255)
                                                                                              tCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGACACCTg                                     <  1:26548/137‑1 (MQ=255)
                                                                                              tCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGACACCTg                                     <  3:7354/137‑1 (MQ=255)
                                                                                                            gATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTaa                                                                                                      >  2:39949/1‑58 (MQ=255)
                                                                                                            gATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTaa                                                                                                      <  1:39949/58‑1 (MQ=255)
                                                                                                            gATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGACACCTGAAATTAAAAGCTGTcc                     >  3:35236/1‑139 (MQ=255)
                                                                                                              tAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGACACCTGAAATTAAAAGCTGTCCgg                   <  7:146996/139‑1 (MQ=255)
                                                                                                                 aGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGACACCTGAAATTAAAAGCTGTCCGGCaa                >  1:113039/1‑139 (MQ=255)
                                                                                                                 aGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGACACCTGAAATTAAAAGCTGTCCGGCaa                >  8:316734/1‑139 (MQ=255)
                                                                                                                       cGATTTTGGCCGGGCGCGTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGCCACCTGAAATTAAAAGCTGTCCGGCAACAGGAt          <  7:360114/139‑1 (MQ=255)
                                                                                                                           tttGGCATGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGACACCTGAAATTAAAAGCTGTCCGGCAACAGGATGCa       <  3:3406/138‑1 (MQ=255)
                                                                                                                             tGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGat                                                                >  8:50565/1‑79 (MQ=255)
                                                                                                                             tGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGat                                                                <  7:50565/79‑1 (MQ=255)
                                                                                                                               gCAGGGCGCTTACGGTGCTCCATGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGACACCTGAAATTAAAAGCTGTCCGGCAACAGGATGCAGTGAg  <  7:71987/139‑1 (MQ=255)
                                                                                                                               gCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGACACCTGAAATTAAAAGCTGTCCGGCAACAGGATGCAGTGAg  <  8:79572/139‑1 (MQ=255)
                                                                                                                               gCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGACACCTGAAATTAAAAGCTGTCCGGCAACAGGATGCAGTGAg  <  5:70017/139‑1 (MQ=255)
                                                                                                                               gCAGGGCGCTTACGGTGCGCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGACACCTGAAATTAAAAGCTGTCCGGCAACAGGATGCAGTGAg  <  5:122471/139‑1 (MQ=255)
                                                                                                                                  gggCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAAGGCTGAAAAACGGa                                                                                        >  6:283253/1‑50 (MQ=255)
                                                                                                                                  gggCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGa                                                                                        <  5:283253/50‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                 
CTTAATTTGGCACCGCTCGGGACGACATACATAACGCCCATGGGCATTTATTCAAAACGATATGAACGAATCCGAGATATATCGCGGGGAGCTGTCGTTTCTGTTCCGGATAAAGCGTTCGATTTTGGCAGGGCGCTTACGGTGCTCCAGGAAGCGGGTTTGCTAACGCTGAAAAACGGATTTAACGGAACAGGTTCAGTTGATATGATTAAAGATAATCCGAGACACCTGAAATTAAAAGCTGTCCGGCAACAGGATGCAGTGAG  >  CP006881/948217‑948482

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: