Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,063,201 G→A 12.2% intergenic (+58/‑149) 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,2010GA12.2% 123.4 / 6.8 49intergenic (+58/‑149)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base A (3/3);  total (25/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 = 7.16e-01

TGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTT  >  NC_002947/2063055‑2063347
                                                                                                                                                    |                                                                                                                                                  
tGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAggg                                                                                                                                                  <  4:244599/149‑1 (MQ=255)
  cTCAAGAAGCTGGAAGAGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAgtgt                                                                                                                                                 <  4:35633/148‑1 (MQ=255)
   tCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAAGGtcc                                                                                                                                               >  2:76163/1‑147 (MQ=255)
          gCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                                        <  2:161994/149‑1 (MQ=255)
          gCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGACCTGCACAAGGTCCGAAcc                                                                                                                                          >  3:259464/1‑147 (MQ=38)
            tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgt                                                                                                                                      >  7:190096/1‑149 (MQ=255)
             ggAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgtg                                                                                                                                     >  2:93259/1‑149 (MQ=255)
                     atatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGt                                                                                                                                                 <  4:60462/129‑1 (MQ=255)
                     atatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGt                                                                                                                                                 >  3:60462/1‑129 (MQ=255)
                      tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCaag                                                                                                                                 >  3:93579/1‑142 (MQ=255)
                      tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGGGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCaag                                                                                                                                 <  4:93579/144‑3 (MQ=255)
                       atCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                                           <  8:36224/149‑1 (MQ=255)
                        tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgtg                                                                                                                          <  8:190096/149‑1 (MQ=255)
                        tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgtg                                                                                                                          <  4:242233/149‑1 (MQ=255)
                        tCCAGACCGCGACCAAGAGCCTGCTCGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTATGACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCCCGCCCACAGGCATTGCCCAGGGTTAGGACCTTGTGCAGGTCCagt                                                                                                                           >  1:298058/1‑149 (MQ=255)
                         ccAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGGCCTGTGc                                                                                                                         >  5:127111/1‑149 (MQ=255)
                          cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCtgtgg                                                                                                                         <  4:259464/148‑2 (MQ=21)
                             aCCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCCCGGGCTTTCCCGCTCCCACAGGCATTGCACCGGGTTCGGACCATGTGCAGGTCCTGTgcgcgc                                                                                                                     >  4:236611/1‑149 (MQ=255)
                             aCCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCAT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGc                                                                                                                     <  1:69321/149‑1 (MQ=18)
                              ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCg                                                                                                                    >  3:45540/1‑149 (MQ=255)
                               cgcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                                   <  6:156739/149‑1 (MQ=18)
                                     cAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGc                                                                                                             >  3:248965/1‑149 (MQ=255)
                                       agagCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCGCCCCCAGGCCTTGCACAGGGTTCGGGCCTTGTGCAGGTGCTGTGCGCGCGGGCCAGccc                                                                                                           >  3:146470/1‑149 (MQ=255)
                                             tgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTTCCCGCGCCCCCAGGCATTGCCCAGGGGTAGGACCATGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAAt                                                                                                     >  3:67382/1‑149 (MQ=255)
                                               ctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACTGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAAt                                                                                                     <  5:88268/147‑1 (MQ=255)
                                                     aaCACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCGCCCACAGGCATTGCACAGGGTTAGGACCTTGTGCAGGTACTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCg                                                                                             >  8:176725/1‑149 (MQ=255)
                                                                cGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTcc                                                                                  <  1:93259/149‑1 (MQ=255)
                                                                 gtggtgCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCAGGCATTGCACAGGGTTAGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTccc                                                                                 >  8:100691/1‑149 (MQ=255)
                                                                 gtggtgCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTccc                                                                                 <  4:45540/149‑1 (MQ=255)
                                                                    gtgCACGCCACACCGAAGTA‑AACCGGGGCCGCGT‑GGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGCACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGa                                                                              <  5:303888/149‑1 (MQ=255)
                                                                       cACGGCACACCGAAGTA‑AACCGGGGCCGCGTATGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACATTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACt                                                                            >  1:33803/1‑149 (MQ=255)
                                                                               acCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACa                                                                       <  6:81218/145‑1 (MQ=255)
                                                                               acCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACa                                                                       >  5:81218/1‑145 (MQ=255)
                                                                                 cGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATcc                                                                 >  3:269080/1‑149 (MQ=255)
                                                                                  gAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATccc                                                                <  4:248965/149‑1 (MQ=255)
                                                                                              gggCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGTGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGt                                                     <  5:140122/149‑1 (MQ=255)
                                                                                                gCCGCGT‑TGCGGCCCATGCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTGCGGACCTGGTCCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCa                                                   <  1:43949/149‑1 (MQ=255)
                                                                                                          cGGCCCATTCGCGGGCTTGCCGGCTCCCACCGGCATGGCACAGGGTTCGGACCTTTTGCTGGGCCTGTGCTCGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCTAGGTCAGTTATTCa                                           <  5:220218/148‑1 (MQ=38)
                                                                                                          cGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGCGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                           <  6:127111/148‑1 (MQ=255)
                                                                                                              ccATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCTGGGCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCAt                                      <  7:135292/149‑1 (MQ=255)
                                                                                                                         cTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTgg                           <  2:33803/149‑1 (MQ=255)
                                                                                                                            gCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgtgt                        >  3:74106/1‑149 (MQ=255)
                                                                                                                              ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                      >  7:157520/1‑149 (MQ=255)
                                                                                                                                gCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt                    >  1:328441/1‑149 (MQ=18)
                                                                                                                                     cacaGGCATTGCACAGGGATCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAACCCGGCGAATGGGCCCAGCACAATATTCCCCGACTCAACATTATCCCCGCCCCAACGTCAGTTATACAGGCATCACACACCTGGTGTTCGTAGACt               >  4:103200/1‑149 (MQ=255)
                                                                                                                                       caGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt               >  7:326843/1‑147 (MQ=255)
                                                                                                                                       caGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga             <  1:139141/149‑1 (MQ=255)
                                                                                                                                        aGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCGGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgag            <  4:269080/149‑1 (MQ=255)
                                                                                                                                          gctttGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga           <  3:236611/145‑1 (MQ=255)
                                                                                                                                                  caGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACTCCTGGTGTTCGTTGTCTGAGAGGGTCCGt   >  7:334959/1‑148 (MQ=255)
                                                                                                                                                   aGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt  >  2:25067/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTT  >  NC_002947/2063055‑2063347

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