Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,063,212 T→C 15.4% intergenic (+69/‑138) 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,2120TC15.4% 91.2 / 9.2 40intergenic (+69/‑138)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (16/17);  new base C (2/4);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

GAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGT  >  NC_002947/2063069‑2063359
                                                                                                                                               |                                                                                                                                                   
gAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgtg                                                                                                                                                 <  6:41560/148‑1 (MQ=255)
gAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGc                                                                                                                                                <  1:302606/149‑1 (MQ=255)
   gCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAAGGTCCGAACCCTGTGCaat                                                                                                                                             <  8:248427/149‑3 (MQ=255)
          tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCtgtg                                                                                                                                      >  1:80931/1‑149 (MQ=21)
                 cgcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                               <  7:111340/149‑1 (MQ=255)
                  gcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCggg                                                                                                                              <  8:152305/149‑1 (MQ=255)
                   cgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGc                                                                                                                             <  3:142415/149‑1 (MQ=18)
                          gagCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCcg                                                                                                                      <  2:80931/149‑1 (MQ=17)
                              ctgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGaa                                                                                                                  <  1:119621/149‑1 (MQ=255)
                                    gATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTGGGCCCTTGTGCAGGTCCTGTGCGAGCGGGCCAGCCCGGGAATGGGcc                                                                                                            >  7:6445/1‑149 (MQ=255)
                                        acacCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCATTTCACAGGGTTCGGACCTTGTGGAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCg                                                                                                         >  3:79656/1‑148 (MQ=255)
                                         cacCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTAGCGGGCTTGCCCCCGCCCCCCGGCAATTAACACGGTTGGGACCTTGTGGAGGTTCAGTGGGAGCGGGCAAGCCCGCGCAAGGGGCCCGcc                                                                                                       >  4:264903/1‑148 (MQ=255)
                                          acCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCCCAGGGTTCGGACCTTGTGGAGGTCCTGTGCGAGCGGGCAAGCCCGCGGATGGGCCCCGcac                                                                                                      >  4:27481/1‑149 (MQ=255)
                                                     ggtgCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCGCCCACAGGCATTTCACACGGTTCCGACCATGTTGAGGGACTGTGCGCGCGGGCAAGCCCCCGAATGGGCCCCGCACAATTTTacccg                                                                                           >  4:45485/1‑147 (MQ=37)
                                                                       gTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGc                                                                         <  3:45485/149‑1 (MQ=255)
                                                                          aaCCGGGGCCGCGTTGCGGCCCATTAGCGGGCGTTCCCCCTCCCACAGGCAATGCACAGGGTTACGACCTTGTGCAGGGGCAGTTCGAGGCGGCAAGCCCGCGAAAGGGCCCCGCACAAATTTACCGGAAAGAACATTATACCCGccc                                                                       >  4:82370/1‑148 (MQ=255)
                                                                                    gggTTGGGGCCCATTCCGGGGGTTCCCCGCTCCCACCGGAATTGCACAGGGTTCGCCCTTGGTGCTGGGCGGGGGCGCGGGGGAAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGtt                                                            <  4:159477/147‑1 (MQ=255)
                                                                                          cGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                                       <  6:290050/148‑1 (MQ=255)
                                                                                               cATTCGCGGGCTTGCCCGCTCCCCCAGGCATTGCACAGGGTTCGGACCTTGTGGAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATGCCCGCCCCAAGGTCAGTTATTCAGGCATc                                                 >  7:132423/1‑149 (MQ=255)
                                                                                                   cgcgGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAgg                                                                  >  5:236599/1‑128 (MQ=255)
                                                                                                   cgcgGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGGCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAgg                                                                  <  6:236599/128‑1 (MQ=255)
                                                                                                             cccGCTCCCACAGGCCTTGCACAGGGCTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGtt                                   <  8:6445/149‑1 (MQ=255)
                                                                                                              ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                  <  3:82370/149‑1 (MQ=255)
                                                                                                              ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                  >  3:112005/1‑149 (MQ=255)
                                                                                                                 cTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                               >  1:212892/1‑149 (MQ=18)
                                                                                                                   cccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACaa                                                                                                    <  1:266511/78‑1 (MQ=255)
                                                                                                                   cccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACaa                                                                                                    >  2:266511/1‑78 (MQ=255)
                                                                                                                    ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                  >  1:295883/1‑143 (MQ=255)
                                                                                                                    ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                  <  2:295883/143‑1 (MQ=255)
                                                                                                                        aGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGTGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAg                                                      <  2:100668/119‑1 (MQ=255)
                                                                                                                        aGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTGCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAg                                                      >  1:100668/1‑119 (MQ=255)
                                                                                                                           cATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg                     <  4:112005/149‑1 (MQ=255)
                                                                                                                             ttGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGt                   >  7:294249/1‑149 (MQ=255)
                                                                                                                               gCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTcc                 >  2:181580/1‑149 (MQ=255)
                                                                                                                                   aGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt              >  3:159909/1‑148 (MQ=255)
                                                                                                                                     ggTGCGGACCTTGTGCAGGTCCTGTGGGAGCGGTCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCaa           <  3:201451/149‑1 (MQ=255)
                                                                                                                                      gtccgaaCCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt          <  3:34494/143‑1 (MQ=37)
                                                                                                                                       ttCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTccc                                                                                             <  2:157525/65‑1 (MQ=255)
                                                                                                                                       ttCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTccc                                                                                             >  1:157525/1‑65 (MQ=255)
                                                                                                                                              cTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGt  >  8:92848/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGT  >  NC_002947/2063069‑2063359

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