Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,063,228 C→G 16.6% intergenic (+85/‑122) PP_1838 → / → PP_1839 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,063,2280CG16.6% 123.6 / 16.6 54intergenic (+85/‑122)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/25);  new base G (3/6);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

CAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCATTTT  >  NC_002947/2063081‑2063372
                                                                                                                                                    |                                                                                                                                                 
cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                                                   <  3:174125/148‑1 (MQ=255)
 aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGa                                                                                                                                                 <  5:211340/149‑1 (MQ=18)
  gACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAg                                                                                                                                                <  5:359465/149‑1 (MQ=18)
   aCCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGc                                                                                                                                               <  5:384026/149‑1 (MQ=255)
   aCCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGc                                                                                                                                               >  1:274780/1‑149 (MQ=255)
     cgcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                                                             <  2:386717/149‑1 (MQ=18)
     cgcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                                             <  3:202096/149‑1 (MQ=255)
     cgcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                                             <  2:274780/149‑1 (MQ=255)
     cgcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCGCCCACAGGCATTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                                                             <  8:123269/149‑1 (MQ=25)
       cgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTTCCCGCTGCCACAGGGATTGCACAGGGTTCGGACCATGTTCACGTACTGTTCGAGCagg                                                                                                                                            >  4:116361/1‑148 (MQ=255)
               agCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAACCCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCcgc                                                                                                                                   <  5:71112/149‑1 (MQ=255)
                gCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCcgc                                                                                                                                   >  6:127688/1‑148 (MQ=255)
                 ccTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCGTGTGCAGGTCCTGTGCGCGCGGGCCAGCCcgcg                                                                                                                                  >  5:6204/1‑148 (MQ=255)
                    gctgGATACCACCGCCGGCGTGGGGCACGGCACACCGAATTACACCGGCGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGCCAGTGCACAGGGTTCGGCCCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATg                                                                                                                              <  1:30098/149‑1 (MQ=37)
                    gctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTGCCACAGGCATTTCACAGGGTTCGGACCTTGTGCAGGTCCTTTGCGAGCGGGCAAGCCCGCGAATg                                                                                                                              >  8:182584/1‑149 (MQ=255)
                    gctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGCAAACCGGGGCCGCGTGGCGGCCCATTCGCGGGCTTGCCCGCTCCCACCGCCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATg                                                                                                                              <  3:18942/149‑1 (MQ=255)
                     ctgGATAACCCCGCCGGCGTGGTGCACGGAACACCGAAGTAACCCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACCGGGTTCGGACCTTGTGCAGGTCCTGTGCTAGCGGGCAAGCCCGCGAATg                                                                                                                              <  6:6204/148‑1 (MQ=38)
                     ctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCAGGCATTGGACAGGGGTGGGCCCTTGTGCAGGGCCGGTGCGAGCGGGCTAGCCCGCGAATg                                                                                                                              >  2:30098/1‑148 (MQ=38)
                     ctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATg                                                                                                                              <  5:127688/148‑1 (MQ=255)
                        gATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATGCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATgg                                                                                                                             <  1:251588/147‑1 (MQ=255)
                             cacCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGACCGCTCCCACAGGCATTGCACAGGGTTAGGACCTTGTGGAGGTCCTGTGCGAGCCGGCAAGCCCGCCAAAGGGCCCCGcc                                                                                                                     >  6:380169/1‑148 (MQ=255)
                              acCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTTCCCGCTCCCCCAGGGCTTGGACAGGGTTCGGGCCTTGTGGAGGGGCGGTGGGCGCGGGCCAGCCCGCGGATGGGGCCCGcc                                                                                                                     >  1:239375/1‑147 (MQ=255)
                                     gCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTc                                                                                                             <  4:85425/149‑1 (MQ=255)
                                            gCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCAAAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACt                                                                                                      >  2:107453/1‑149 (MQ=255)
                                             cACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTCCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTg                                                                                                     <  5:380169/149‑1 (MQ=255)
                                             cACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCACCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACt                                                                                                      >  6:311707/1‑148 (MQ=255)
                                                     acCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACACGCAGTGCACAGGGGTCGGCCCTTGTGCCGGGCCTGGGCGAGCGGGCACGCCCGCGAATGGGGCCCGCCCAAGTTTCCCGGCCGGACCAGTAt                                                                                            >  7:126387/1‑149 (MQ=25)
                                                        gAAGTAACCGGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCGCCCACAGCCATTGCACCGGGTGCGCCCCTTGTGCAGGTCCTGGGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATccc                                                                                         <  8:213814/149‑1 (MQ=37)
                                                                   gggCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGt                                                                              <  7:182584/149‑1 (MQ=255)
                                                                   gggCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGt                                                                              <  5:311707/149‑1 (MQ=255)
                                                                     cccGCGTTGCGGCCCATTCGCGGGCGTGCCCGCTCCCACAGGCATGCCACAGGGTGCGCCCCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATTCCCGCCCCAAGGTCa                                                                            <  6:301386/148‑1 (MQ=38)
                                                                     cccGCGTGGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGCACCTTGTCCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGa                                                                                                        <  7:120849/121‑1 (MQ=255)
                                                                                     ttCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGCCTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATcac                                                            >  5:140953/1‑149 (MQ=14)
                                                                                        gcgGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATcacaca                                                         >  1:354211/1‑149 (MQ=255)
                                                                                              ttGCCCGCTCCCACAGGCTTTGCACTGGGTGCGGACGTTGTGCAGGTCGTGTGCGAGCGGGAAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgt                                                   <  5:206649/149‑1 (MQ=255)
                                                                                                gCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgtgt                                                 >  8:206854/1‑149 (MQ=255)
                                                                                                  ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                                                    <  1:153087/128‑1 (MQ=14)
                                                                                                  ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                                                    >  2:153087/1‑128 (MQ=14)
                                                                                                    gCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt                                             <  2:239375/149‑1 (MQ=255)
                                                                                                     cTCCCATAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGCAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt                                             <  1:107453/149‑1 (MQ=255)
                                                                                                         cacaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                                        >  8:107256/1‑149 (MQ=255)
                                                                                                           caGGCCCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                                      <  6:140953/149‑1 (MQ=21)
                                                                                                            aGGCATTGCACCGGGTTCGGACCTTGTGCTGGTCCTGTGCGAGCGGGAAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgag                                     <  2:354211/149‑1 (MQ=255)
                                                                                                            aGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgag                                     >  5:43139/1‑149 (MQ=255)
                                                                                                                  tGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTc                               >  6:68140/1‑149 (MQ=255)
                                                                                                                    cacaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCCGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGTTCCg                             >  6:215879/1‑149 (MQ=255)
                                                                                                                     acaGGGTTCGGACCTTCTGCAGGTCCTGACCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAAATTTCCCGGACTGAACATTAACCCCGCCCCAAGGTCAGTTTTTCAGGCATCACACACCAGGTGTACGTTGACAGAGAGCGTCCGt                            >  2:390978/1‑149 (MQ=255)
                                                                                                                      ctgggTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATTGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt                           <  3:116361/147‑1 (MQ=255)
                                                                                                                             cGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGa                     <  7:107256/148‑1 (MQ=255)
                                                                                                                              ggcccTTGTGCAGGGCCGGGGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGaaa                   <  5:6482/146‑1 (MQ=255)
                                                                                                                                aCCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAaga                 <  7:206854/149‑1 (MQ=255)
                                                                                                                                  cTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCAGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGt               >  3:285633/1‑149 (MQ=255)
                                                                                                                                       gCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTggg            >  7:23241/1‑147 (MQ=255)
                                                                                                                                         cAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGaa          <  3:342235/148‑1 (MQ=255)
                                                                                                                                                tgtgCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCAtttt  >  1:336083/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCATTTT  >  NC_002947/2063081‑2063372

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