Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,007,080 C→G 10.0% intergenic (+40/+38) PP_2627 → / ← PP_2628 membrane protein/ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,007,0800CG10.0% 148.5 / 6.6 50intergenic (+40/+38)PP_2627/PP_2628membrane protein/ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base G (3/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

AATCCCATTACCCAGCAACCTTTCAAAGGCTTCCGCTTGCCCCGAGGCTTGTTGCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCTTTTTCGCCAACAGTTCGGCATGGGT  >  NC_002947/3006933‑3007225
                                                                                                                                                    |                                                                                                                                                  
aaTCCCATTACCCAGCAACCTTTCAAAGGCTTCCGCTTGCCCCGAGGCTTGTTTCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTc                                                                                                                                                    <  6:141463/148‑1 (MQ=255)
aaTCCCATTACCCAGCAACCTTTCAAAGGCTTCCGCTTGCCCCGAGGCTTGTTGCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATCGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTc                                                                                                                                                    <  1:94773/149‑1 (MQ=255)
       ttACCCAGCAACCTTTCAAAGGCTTCCGCTTGCCCCGAGGCTTGTTGCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTg                                                                                                                                            <  8:121009/149‑1 (MQ=255)
          cccAGCAACCTTTCAAAGGCTTCCGCTTGCCCCGAGGCTTGTTGCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGt                                                                                                                                         <  2:305921/149‑1 (MQ=255)
                aaCCTTTCAAAGGCTTCCGCTTGCCCCGAGGCTTGTTGCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATgcgc                                                                                                                                   >  5:56134/1‑149 (MQ=255)
                   cTTTCAAAGGCTTCCGCTTGCCCCGAGGCTTGTTGCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCgctgg                                                                                                                                >  2:142412/1‑149 (MQ=255)
                    tttCAAAGGCTTCCGCTTGCCCCGAGGCTTGTTGCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCt                                                                                                                                             <  1:97314/135‑1 (MQ=255)
                    tttCAAAGGCTTCCGCTTGCCCCGAGGCTTGTTGCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCt                                                                                                                                             >  2:97314/1‑135 (MQ=255)
                              ttCCGCTTGCCCCGAGGCTTGTTGCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGc                                                                                                                     >  1:181749/1‑149 (MQ=255)
                                      gCCCCGAGGCTTGTTGCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGc                                                                                                               <  1:142412/147‑1 (MQ=255)
                                                ttgttgCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTAc                                                                                                          >  4:6379/1‑141 (MQ=255)
                                                ttgttgCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTAc                                                                                                          <  3:6379/141‑1 (MQ=255)
                                                ttgttgCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGc                                                                                                  >  2:110167/1‑149 (MQ=255)
                                                   ttgCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGc                                                                                                                                               >  8:314051/1‑102 (MQ=255)
                                                   ttgCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGc                                                                                                                                               <  7:314051/102‑1 (MQ=255)
                                                   ttgCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCcac                                                                                                     <  7:15588/143‑1 (MQ=255)
                                                   ttgCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCcac                                                                                                     >  8:15588/1‑143 (MQ=255)
                                                     gCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCg                                                                                              <  6:56134/148‑1 (MQ=255)
                                                          agccagccAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCAc                                                                                         >  8:277479/1‑148 (MQ=255)
                                                          agccagccAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCAc                                                                                         <  7:277479/148‑1 (MQ=255)
                                                                      tggtgCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCcaca                                                                                                    <  7:25149/125‑1 (MQ=255)
                                                                      tggtgCAGGCGGTGCTACCGGATGGCCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCcaca                                                                                                    >  8:25149/1‑125 (MQ=255)
                                                                         tgCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACCTCCACAGCCTCGGCCACGCCCTGATCCTCAcc                                                                          <  1:40941/149‑1 (MQ=255)
                                                                          gCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGTCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCATCAcc                                                                          >  8:71079/1‑147 (MQ=255)
                                                                          gCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCg                                                                                              <  3:229208/127‑1 (MQ=255)
                                                                          gCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCg                                                                                              >  4:229208/1‑127 (MQ=255)
                                                                               cGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAAc                                                                   <  2:181749/149‑1 (MQ=255)
                                                                                          gATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCActct                                                         >  3:312455/1‑148 (MQ=255)
                                                                                              gggCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGAtg                                                     >  3:86606/1‑149 (MQ=255)
                                                                                                gCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGAtgct                                                   >  3:273498/1‑149 (MQ=255)
                                                                                                gCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGAtgct                                                   <  5:63006/149‑1 (MQ=255)
                                                                                                       cTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCAccc                                                                         >  2:244434/1‑120 (MQ=255)
                                                                                                       cTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCAccc                                                                         <  1:244434/120‑1 (MQ=255)
                                                                                                        tATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAAc                                           >  1:130256/1‑149 (MQ=255)
                                                                                                                         gcgcAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCtt                          >  5:277094/1‑149 (MQ=255)
                                                                                                                            cAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACTGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCttttt                       <  1:110167/149‑1 (MQ=255)
                                                                                                                                 cATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCTTTTTCGcc                   <  6:277094/148‑1 (MQ=255)
                                                                                                                                   ttGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCCCCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCttttt                       <  2:35120/142‑1 (MQ=255)
                                                                                                                                   ttGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCTCAACCGTGCATAGGTACCCTTTTTCGCCaa                 >  5:61393/1‑148 (MQ=255)
                                                                                                                                   ttGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCttttt                       >  1:35120/1‑142 (MQ=255)
                                                                                                                                         acccAGGCTTCGAAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCTTTTTCGCCAACa               >  4:13459/4‑144 (MQ=255)
                                                                                                                                         acccAGGCTTCGAAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCTTTTTCGCCAACa               <  3:13459/141‑1 (MQ=255)
                                                                                                                                           ccAGGCTTCGAAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCCCTCTGATGCTGCCACAACCGTGCATAGGTACCCTTTTTCGCCAACAGTTCgg         >  2:166690/2‑148 (MQ=255)
                                                                                                                                           ccAGGCTTCGAAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGAtgctgc                                                 <  4:71236/107‑1 (MQ=255)
                                                                                                                                           ccAGGCTTCGAAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGAtgctgc                                                 >  3:71236/2‑108 (MQ=255)
                                                                                                                                             aGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCTCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCTTTTTCGCCAACAGTTCGGCAt      <  3:335792/149‑1 (MQ=255)
                                                                                                                                             aGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACGCAGTATCCCGCCAAGCTGATGCTGCCACAACCGTGCATAGGTACCCTTTTTCGCCAACAGTTCGGCAt      <  5:239656/149‑1 (MQ=255)
                                                                                                                                               gCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCTTTTTCGCCAACAGTTCGGCATgg    <  2:130256/149‑1 (MQ=255)
                                                                                                                                                 tttCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCTTTTTCGCCAACAGTTCGGCATgg    >  3:62317/1‑147 (MQ=255)
                                                                                                                                                 tttCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCTTTTTCGCCAACAGTTCGGCATGGGt  >  8:175721/1‑149 (MQ=255)
                                                                                                                                                 tttCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCTTTTTCGCCAACAGTTCGGCATGGGt  >  1:205592/1‑149 (MQ=255)
                                                                                                                                                   tCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCTTTTTCGCCAACAGTTCGGCATgg    <  4:62317/145‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AATCCCATTACCCAGCAACCTTTCAAAGGCTTCCGCTTGCCCCGAGGCTTGTTGCAACAGCCAGCCAGTGTGGTGCAGGCGGTGCTACCGGATGGGCAGCGACTATGAACCTCGTCAGGGGCGCAGTGCATTGACTGTGCAGGCTTTCGAAGCCTGGGTATGCGCTGGCTGCCCCCGGCCTCTGCTTACTCCACAGCCTCGGCCACGCCCTGATCCTCACCCAGGAACCCGCCACTCTGATGCTGCCACAACCGTGCATAGGTACCCTTTTTCGCCAACAGTTCGGCATGGGT  >  NC_002947/3006933‑3007225

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: