Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,063,219 G→A 9.1% intergenic (+76/‑131) 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,2190GA 9.1% 154.3 / 4.6 55intergenic (+76/‑131)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base A (3/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

ATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAG  >  NC_002947/2063078‑2063365
                                                                                                                                             |                                                                                                                                                   
atCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                                                              >  1:257664/1‑149 (MQ=255)
atCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                                                              >  1:269086/1‑149 (MQ=255)
   cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                                            >  7:375721/1‑148 (MQ=255)
   cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                                            <  3:317187/148‑1 (MQ=255)
   cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCATGCCCCCGACCCCAGGCATTGCACAGGGTTCGGAAGATGTGGAGGTCCtgtggg                                                                                                                                           >  4:32047/1‑147 (MQ=255)
        cgcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                                     <  8:341630/149‑1 (MQ=255)
         gcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCggg                                                                                                                                    <  5:127239/149‑1 (MQ=255)
               aaGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGcc                                                                                                                              <  3:32047/149‑1 (MQ=255)
                  agCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                                     <  1:65120/139‑1 (MQ=255)
                  agCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                                     >  2:65120/1‑139 (MQ=255)
                        ctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAAt                                                                                                                       <  8:375721/147‑1 (MQ=255)
                                   cgccgGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCGCCCACAGGCATTGCCCAGGGGTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCCAGCCCGCGCATGGGCCCCGCCCaa                                                                                                          >  3:78841/1‑149 (MQ=38)
                                           tggtgCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCAGGCATTGCACAGGGGTCGGACCTTGTGCAGGTCGTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCg                                                                                                  <  8:220912/149‑1 (MQ=255)
                                                cACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCCTTGCCCAGGGTTAGGACCTTTTGCAGGTACTGTGCGAGCGGGCAAGCCCGCGAAAGGGCCCCGCACAATTTTGCCGGGCt                                                                                              >  7:226444/1‑148 (MQ=37)
                                                  cGGCACACCGAAGTAAACCGGGGCCGGGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCAGGCATTGCACAGGGTTGGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGaa                                                                                           <  2:137058/149‑1 (MQ=255)
                                                  cGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGcc                                                                                                                  <  6:83040/126‑1 (MQ=255)
                                                  cGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGcc                                                                                                                  >  5:83040/1‑126 (MQ=255)
                                                    gCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCCGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTGCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACa                                                                                         <  1:122867/149‑1 (MQ=255)
                                                        acCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGa                                                                                                                                          >  7:364721/1‑97 (MQ=255)
                                                        acCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGa                                                                                                                                          <  8:364721/97‑1 (MQ=255)
                                                          cGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCt                                                                                                                                                <  6:130064/89‑1 (MQ=255)
                                                          cGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCt                                                                                                                                                >  5:130064/1‑89 (MQ=255)
                                                            aaGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTTCCCGCTCCCACAGGCATTGCACAGGGTTAGGACCATTTGGAGGTGCGGTGCGAGCCGGCAAGCCCGCGAATGGGCCCCGCACAATTTTACCGGACTGAAAATTATgccc                                                                                 >  4:361909/1‑149 (MQ=37)
                                                                     ggggCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTGCTGTGCGAGCGGGCAAGCCCGCGAATGGGGCCCGCGCCATTTTGCCGGCCTGCACATTATGCCCGCCCCAAgg                                                                        >  1:254057/1‑149 (MQ=37)
                                                                          cgcgTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGt                                                                   <  2:82698/149‑1 (MQ=255)
                                                                                 cGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                                             <  2:257664/148‑1 (MQ=255)
                                                                                  ggCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGCGCGGGCAAGCCCGCGAATGGGTCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAgg                                                           <  2:269086/149‑1 (MQ=255)
                                                                                     acATTCGCGGGCTTGCCCGCTCCCACAGGCTTTGCACAGGGTTCGGACCTTGTGCTGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCAt                                                        <  4:78841/148‑1 (MQ=255)
                                                                                            cgGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCacacac                                                 <  3:361909/149‑1 (MQ=255)
                                                                                                   gCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACa                                                                                         >  8:290349/1‑102 (MQ=255)
                                                                                                   gCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACa                                                                                         <  7:290349/102‑1 (MQ=255)
                                                                                                     ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCGCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                        >  6:138023/1‑149 (MQ=17)
                                                                                                     ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGTGCGGGAAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                        <  8:226444/149‑1 (MQ=255)
                                                                                                     ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATc                                                       <  1:56285/134‑1 (MQ=255)
                                                                                                     ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATc                                                       >  2:56285/1‑134 (MQ=255)
                                                                                                     ccGCTCCCACAGGCATTGCACAGGGTTAGGACCTTGTGCAGGACCTGTGCAAGCGGGCAAACCCGCGAAAGGGCCCCGCACAATTTTACCGGACTGAACATTATCCCCGACCCAAGGTCAGTTATTCAGGCATAACAAAACTGGTGttg                                        >  7:61095/1‑148 (MQ=255)
                                                                                                       gCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                        >  4:172627/1‑147 (MQ=18)
                                                                                                       gCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                        <  3:172627/147‑1 (MQ=18)
                                                                                                            cacaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCTCTAATGGGCCCCGTACAATTTTCCCGTACTGAACATTATCCCCGTTTCAAGGTTAGTTATTCATGCATCATACACTTGGTGTTCGTTGACt                                 <  8:61095/149‑1 (MQ=255)
                                                                                                            cacaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                                 >  6:424110/1‑149 (MQ=255)
                                                                                                             acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTc                                                              >  7:162137/1‑119 (MQ=255)
                                                                                                             acaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTc                                                              <  8:162137/119‑1 (MQ=255)
                                                                                                                ggCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                               >  4:215799/1‑147 (MQ=255)
                                                                                                                  cATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg                           >  4:218255/1‑149 (MQ=255)
                                                                                                                    ttGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGTGCGGGAAAGCCCGCGAATGGGCCCCGCTCAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGt                         <  2:254057/149‑1 (MQ=255)
                                                                                                                        acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGt                     >  6:71852/1‑149 (MQ=255)
                                                                                                                        acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGt                     >  8:260801/1‑149 (MQ=255)
                                                                                                                        acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAACCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTAACCCCGCCCCAAGGTAAGTTATTCAGGCATCAAACACCTGGTGTTCGATGACTTAGAGGGTCCGt                     >  8:253288/1‑149 (MQ=255)
                                                                                                                          aaggTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt                    >  1:274085/3‑149 (MQ=37)
                                                                                                                          aGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCGTCACACACCTGGTGTTc                                        <  6:337490/128‑1 (MQ=255)
                                                                                                                          aGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCGTCACACACCTGGTGTTc                                        >  5:337490/1‑128 (MQ=255)
                                                                                                                             gTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt                <  5:424110/149‑1 (MQ=255)
                                                                                                                                 gaaCCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGaaa            <  2:274085/147‑1 (MQ=38)
                                                                                                                                        gtgCAGGTCCTGTGCGAGCGGGCAAGCGCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTg       <  7:253288/147‑1 (MQ=255)
                                                                                                                                           cAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATGCAGACATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAg  >  7:84665/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
ATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAG  >  NC_002947/2063078‑2063365

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