Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,555,648 T→A 10.7% intergenic (‑124/+47) PP_4887 ← / ← PP_4888 hypothetical protein/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,555,6480TA10.7% 111.2 / 5.2 47intergenic (‑124/+47)PP_4887/PP_4888hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (15/26);  new base A (2/3);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

GGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCG‑‑‑CCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCGTGGTCTGCATCACCACCTCAGCCGAGCGTT  >  NC_002947/5555513‑5555795
                                                                                                                                          |                                                                                                                                                    
gacatgagcaactccttttgtcatgtgtcgccgacagcatagcaacaccgctcaggtgtgggagcggccttgtgtcgcgaaagggccgcatagcggcccCCGTCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCATAGCAgcccc                                                                                                                                            <  1:279606‑M2/50‑2 (MQ=255)
gacatgagcaactccttttgtcatgtgtcgccgacagcatagcaacaccgctcaggtgtgggagcggccttgtgtcgcgaaagggccgcatagcggcccCCGTCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCATAGCAgcccc                                                                                                                                            <  1:284322‑M2/50‑2 (MQ=255)
   catgagcaactccttttgtcatgtgtcgccgacagcatagcaacaccgctcaggtgtgggagcggccttgtgtcgcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCATAgcagccccagg                                                                                                                                         <  6:22320‑M2/53‑5 (MQ=255)
     tgagcaactccttttgtcatgtgtcgccgacagcatagcaacaccgctcaggtgtgggagcggccttgtgtcgcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTcgc                                                                                                                                       <  2:30669‑M2/55‑1 (MQ=255)
      gagcaactccttttgtcatgtgtcgccgacagcatagcaacaccgctcaggtgtgggagcggccttgtgtagcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTcgcg                                                                                                                                      >  7:339884‑M2/94‑149 (MQ=255)
      gagcaactccttttgtcatgtgtcgccgacagcatagcaacaccgctcaggtgtgggagcggccttgtgtagcgaaagggccgcaaagcggccCCGGGCATTTGCGCGGGTTTGGAAAAAACCAGGGGCGGCGATGCAGCCCGTTcgcg                                                                                                                                      >  3:264120‑M2/94‑149 (MQ=255)
               cTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTTTAGCCAAATGGCCCCATAGCGGCCCCGGGAATTTTCGCGGGTTTTCCACCAACCAGGGGCGGGTATGCCGCCCGTTCGCGGCACACggg                                                                                                                             >  5:202133/1‑148 (MQ=255)
                  ttgtcatgtgtcgccgacagcatagcaacaccgctcaggtgtgggagcggccttgtgtcgcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGccc                                                                                                                           <  7:70982‑M2/67‑1 (MQ=255)
                      cctgtgtcgccgacagcatagcaacaccgctcaggtgtgggagcggccttgtggcgcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACaaa                                                                                                                      <  3:40526‑M2/72‑1 (MQ=255)
                                 gacagcatagcaacaccgctcaggtgtgggagcggccttgtgtcgcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAg                                                                                                           <  5:210647‑M2/83‑1 (MQ=255)
                                  acagcatagcaacaccgctcaggtgtgggagcggccttgtgtcgcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAgg                                                                                                          <  7:241768‑M2/84‑1 (MQ=255)
                                      catagcaacaccgctcaggtgtgggagcggccttgtgtcgcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCtt                                                                                                       >  2:45222‑M2/62‑148 (MQ=255)
                                            aacaccgctcaggtgtgggagcggccttgtgtcgcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGccc                                                                                               <  4:347180‑M2/94‑1 (MQ=255)
                                             acaccgctcaggtgtgggagcggccttgtgtcgcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTc                                                                                                                                         <  5:31993‑M2/53‑1 (MQ=255)
                                             acaccgctcaggtgtgggagcggccttgtgtcgcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTc                                                                                                                                         >  6:31993‑M2/55‑107 (MQ=255)
                                                              gagcggccttgtgtcgcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTgg                                                                             <  3:309942‑M2/112‑1 (MQ=255)
                                                                gcggccttgtgtcgcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGTAGGCTTCAGGCCCGGAACTGGCCCAGGCTGgcg                                                                           <  1:45222‑M2/114‑1 (MQ=255)
                                                                 cggccttgtgccgctaacgggccgcatcgcgcccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGgcgc                                                                          <  5:250424‑M2/115‑1 (MQ=255)
                                                                     cttgtgtcgcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAg                                                                      <  4:264120‑M2/119‑1 (MQ=255)
                                                                             gcgaaagggccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccg                                                              >  6:21730‑M2/23‑149 (MQ=255)
                                                                                      ccgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcat                                                     >  6:254020‑M2/14‑149 (MQ=255)
                                                                                       cgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcatc                                                    <  5:216504‑M2/137‑1 (MQ=255)
                                                                                       cgcatagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcatc                                                    <  7:322852‑M2/137‑1 (MQ=255)
                                                                                          atagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAg                                                 <  6:202133‑M2/140‑1 (MQ=255)
                                                                                           tagcggccCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGc                                                <  8:339884‑M2/141‑1 (MQ=255)
                                                                                                   ccGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCAcc                                             <  4:14886/144‑1 (MQ=255)
                                                                                                   ccGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCAcc                                             >  3:14886/1‑144 (MQ=255)
                                                                                                    cGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGCCGCGCACATGAGCCGCCCGATCATCAAGCACCTTGCt                                        >  3:144803/1‑148 (MQ=255)
                                                                                                     ggCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTAc                                      <  5:21730/149‑1 (MQ=255)
                                                                                                     ggCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTAc                                      >  7:311246/1‑149 (MQ=255)
                                                                                                      gCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg                                                                              <  4:200304/108‑1 (MQ=255)
                                                                                                      gCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg                                                                              >  3:200304/1‑108 (MQ=255)
                                                                                                            tCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCg                                <  7:188837/148‑1 (MQ=255)
                                                                                                             ccGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTgg                                   >  2:273803/1‑144 (MQ=255)
                                                                                                             ccGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTgg                                   <  1:273803/144‑1 (MQ=255)
                                                                                                             ccGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCGTg                              >  5:65685/1‑149 (MQ=255)
                                                                                                                cgTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCGTGGTc                           <  1:173865/149‑1 (MQ=255)
                                                                                                                    tATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCGTGGTCTGc                        <  1:355079/148‑1 (MQ=255)
                                                                                                                      tGCATCAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCTGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCGTGGTCTGCAt                      <  1:103678/148‑1 (MQ=255)
                                                                                                                        caacaaACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCGTGGTCTGCATcac                   <  1:196652/149‑1 (MQ=255)
                                                                                                                         aacaaACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCg                                >  4:277347/1‑135 (MQ=255)
                                                                                                                         aacaaACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCg                                <  3:277347/135‑1 (MQ=255)
                                                                                                                         aacaaACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCGTGGTCTGCATcacc                  >  3:10001/1‑149 (MQ=255)
                                                                                                                               gaaagggCTGCATAGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCGTGGTCTGCATCACCACCt              >  3:149858/5‑147 (MQ=255)
                                                                                                                                cTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCGTGGTCTGCATCACCACCTCAg           >  7:37955/1‑149 (MQ=255)
                                                                                                                                aaagggCTGCATAGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCGAGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCGTGGTCTGCATCACCACCt              >  1:111144/4‑147 (MQ=255)
                                                                                                                                    ggCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCGTGGTCTGCATCACCACCTCAGCCGa       <  5:254020/149‑1 (MQ=255)
                                                                                                                                       tgctATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCGTGGTCTGCATCACCACCTCAGCCGAGCg    >  8:12690/1‑149 (MQ=255)
                                                                                                                                         ctATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTAATGGCCGTGGTCTGCATCACAACCTCAGACGATCGtt  >  6:175009/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
GGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCG‑‑‑CCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCACCTTGCTACTGGCCGTGGTCTGCATCACCACCTCAGCCGAGCGTT  >  NC_002947/5555513‑5555795

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