Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,521,817 C→T 9.1% intergenic (‑32/‑35) lip ← / → PP_4855 lipase/osmotically‑inducible lipoprotein OsmE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,521,8170CT 9.1% 124.1 / 2.5 44intergenic (‑32/‑35)lip/PP_4855lipase/osmotically‑inducible lipoprotein OsmE
Reads supporting (aligned to +/- strand):  ref base C (21/19);  new base T (2/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAGCATGCCTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTC  >  NC_002947/5521673‑5521962
                                                                                                                                                |                                                                                                                                                 
ccGAGCATGCCTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTc                                                                                                                                                   <  5:262433/145‑1 (MQ=255)
ccGAGCATGCCTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCaaaa                                                                                                                                               >  5:9842/1‑149 (MQ=255)
 cGAGCATGCCTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAAc                                                                                                                                              <  8:4907/149‑1 (MQ=255)
  gAGCATGCCTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAAc                                                                                                                                              >  4:178766/1‑148 (MQ=255)
   aGCATGCCTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACca                                                                                                                                            <  4:102988/149‑1 (MQ=255)
    gCATGCCTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACcac                                                                                                                                           <  7:310664/149‑1 (MQ=255)
      aTGCCTGGCACCAGCACCAGCGGGTATCCCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACcaca                                                                                                                                          >  1:207222/1‑148 (MQ=255)
          cTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGAc                                                                                                                                     >  5:235230/1‑149 (MQ=255)
               accagcagcagcGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTGCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATGATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCttt                                                                                                                                <  3:251720/149‑1 (MQ=255)
               accagcaccagcGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCttt                                                                                                                                <  5:368795/149‑1 (MQ=255)
                    caccagcGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAg                                                                                                                           <  2:207222/149‑1 (MQ=255)
                         gcGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGACCTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCg                                                                                                                      <  5:373468/149‑1 (MQ=255)
                         gcGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCg                                                                                                                      <  6:9842/149‑1 (MQ=255)
                                tCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCAt                                                                                                               <  8:3199/149‑1 (MQ=255)
                                  gcgTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGCCGGGCCAACACTAAAGCGGCCATCATGACCATCGCATTCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATg                                                                                                              >  3:80235/1‑148 (MQ=255)
                                       gCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACaa                                                                                                         <  3:317276/148‑1 (MQ=255)
                                         cAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAg                                                                                                        >  3:120838/1‑147 (MQ=255)
                                          aaaTCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAg                                                                                                        >  1:118984/1‑146 (MQ=255)
                                          aaaTCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAg                                                                                                     >  4:207145/1‑149 (MQ=255)
                                           aaTCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTTTGAACcacagcgcg                                                                                                                                     <  7:76271/116‑6 (MQ=255)
                                                tgctgcATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTg                                                                                               >  5:285065/1‑149 (MQ=255)
                                                     cATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGt                                                                                                             >  2:139989/1‑130 (MQ=255)
                                                          ccGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTc                                                                                     <  3:178766/149‑1 (MQ=255)
                                                                   cTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAgc                                                                             <  3:207145/148‑1 (MQ=255)
                                                                    tttGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCaa                                                                               >  7:205462/1‑145 (MQ=255)
                                                                    tttGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCaa                                                                               <  8:205462/145‑1 (MQ=255)
                                                                      tGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGcc                                                                          <  2:118984/148‑1 (MQ=255)
                                                                       gCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCa                                                                         <  6:235230/148‑1 (MQ=255)
                                                                        cGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCAcc                                                                       <  6:285065/149‑1 (MQ=255)
                                                                                     ccAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCg                                                          >  2:3559/1‑149 (MQ=255)
                                                                                            tACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAgccg                                                   >  1:320930/1‑149 (MQ=255)
                                                                                                gCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAgc                                                                             <  6:83253/119‑1 (MQ=255)
                                                                                                            aCCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGt                                    <  1:3559/148‑1 (MQ=255)
                                                                                                              cATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGa                                 >  7:37681/1‑149 (MQ=255)
                                                                                                                tCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGAtt                               <  2:320930/149‑1 (MQ=255)
                                                                                                                   cAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACg                            >  5:71864/1‑149 (MQ=255)
                                                                                                                     aTCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGt                           >  3:272606/1‑148 (MQ=255)
                                                                                                                              cTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACcgc                 >  2:42382/1‑149 (MQ=255)
                                                                                                                               tGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACcgc                 >  7:261747/1‑148 (MQ=255)
                                                                                                                                  cgcgcgcTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGAt              >  8:123791/1‑148 (MQ=255)
                                                                                                                                   gcgcgcTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGAt              >  6:344357/1‑147 (MQ=255)
                                                                                                                                         tgtgGTTTTGAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGt                           <  2:72504/128‑1 (MQ=255)
                                                                                                                                         tgtgGTTTTGAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGt                           >  1:72504/1‑128 (MQ=255)
                                                                                                                                         tgtgGTTTTGAACCACACAGAACCTTTCCAAGGAGCTGGACCATGTACAACCAGACACTGGCAATCCTTCTGGCAAGCGCCACACTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCt      >  4:124985/1‑149 (MQ=255)
                                                                                                                                              ttCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTc  >  2:351426/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CCGAGCATGCCTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTC  >  NC_002947/5521673‑5521962

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