Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,125,301 A→T 12.2% intergenic (+30/+51) PP_3628 → / ← PP_3629 choline/carnitine/betaine transporter family protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,125,3010AT12.2% 105.5 / 6.3 41intergenic (+30/+51)PP_3628/PP_3629choline/carnitine/betaine transporter family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/17);  new base T (2/3);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGTGGC  >  NC_002947/4125156‑4125449
                                                                                                                                                 |                                                                                                                                                    
tCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAag                                                                                                                                                    <  8:107723/148‑1 (MQ=255)
   agaAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTagcagc                                                                                                                                                >  2:130979/1‑149 (MQ=255)
   agaAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAagcagc                                                                                                                                                >  5:32951/1‑149 (MQ=255)
       ggCCGTACGTTATGGCGCGGTGATGCAGGGGGATGAGTTCTACATTGTGGAGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTAGCAGCCCGt                                                                                                                                            <  3:156600/149‑1 (MQ=255)
        gCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTAGCAGCCCGtg                                                                                                                                           <  1:130979/149‑2 (MQ=255)
         ccGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTAGCAGCCCGtgg                                                                                                                                          <  1:234128/149‑3 (MQ=255)
          cGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTcg                                                                                                                                         >  7:230495/1‑149 (MQ=255)
                    ggCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATg                                                                                                                                    <  8:185116/144‑1 (MQ=255)
                    ggCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATg                                                                                                                                    >  7:185116/1‑144 (MQ=255)
                    ggCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCgaga                                                                                                                               >  4:97049/1‑149 (MQ=255)
                       gcgGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACcg                                                                                                                            >  4:165021/1‑149 (MQ=255)
                             aaGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGTGAGACCGCGCctg                                                                                                                      >  6:121372/1‑149 (MQ=255)
                               gCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCctgct                                                                                                                    >  2:58123/1‑149 (MQ=255)
                                                aCATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCtt                                                                                                   >  2:86380/1‑149 (MQ=255)
                                                 cATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGc                                                                                                     <  3:257695/146‑1 (MQ=255)
                                                        gtgCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTc                                                                                           <  3:97049/149‑1 (MQ=255)
                                                        gtgCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTc                                                                                           <  6:32951/149‑1 (MQ=255)
                                                            aGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAc                                                                                       >  1:269305/1‑149 (MQ=255)
                                                              gagaAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGccc                                                                                                                                              >  8:180560/1‑92 (MQ=255)
                                                              gagaAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGccc                                                                                                                                              <  7:180560/92‑1 (MQ=255)
                                                                gaagaaGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACAt                                                                                                             >  8:95434/1‑123 (MQ=255)
                                                                gaagaaGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACAt                                                                                                             <  7:95434/123‑1 (MQ=255)
                                                                     agGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAcgccgccgc                                                                               <  2:269305/148‑1 (MQ=255)
                                                                       gCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTa                                                                             >  6:17631/1‑148 (MQ=255)
                                                                             gCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGc                                                                      <  5:17631/149‑1 (MQ=255)
                                                                             gCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGc                                                                      <  1:58123/149‑1 (MQ=255)
                                                                                  gACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACAt                                                                 >  4:117779/1‑149 (MQ=255)
                                                                                          cAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTc                                                                                           <  7:109445/115‑1 (MQ=255)
                                                                                          cAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTc                                                                                           >  8:109445/1‑115 (MQ=255)
                                                                                          cAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTg                                                         <  8:31769/149‑1 (MQ=255)
                                                                                                 cGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCt                                                  >  4:260441/1‑149 (MQ=255)
                                                                                                 cGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCt                                                  >  4:243774/1‑149 (MQ=255)
                                                                                                         cccGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATg                                          <  3:260441/149‑1 (MQ=255)
                                                                                                          ccGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGc                                                        >  6:14776/1‑134 (MQ=255)
                                                                                                                   aCGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGTGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAac                                 <  5:121372/148‑1 (MQ=255)
                                                                                                                          cAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCtt                         <  8:230495/149‑1 (MQ=255)
                                                                                                                              ccGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGgcc                     >  2:188831/1‑149 (MQ=255)
                                                                                                                                      cACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACg             >  5:10155/1‑149 (MQ=255)
                                                                                                                                       aCGGGCTGCTTAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACgg            >  2:219319/1‑149 (MQ=255)
                                                                                                                                            ctgctAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCggg       <  3:243774/149‑1 (MQ=255)
                                                                                                                                               ctAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGt     <  3:117779/148‑1 (MQ=255)
                                                                                                                                                 aaGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGTGGc  <  3:165021/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGTGGC  >  NC_002947/4125156‑4125449

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