Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,063,221 T→G 23.1% intergenic (+78/‑129) 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,2210TG23.1% 60.9 / 10.7 38intergenic (+78/‑129)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/14);  new base G (4/5);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCAT  >  NC_002947/2063079‑2063369
                                                                                                                                              |                                                                                                                                                    
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:112005/1‑149 (MQ=255)
                                                                                                    ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                            <  3:82370/149‑1 (MQ=255)
                                                                                                       cTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                                         >  1:212892/1‑149 (MQ=18)
                                                                                                         cccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACaa                                                                                                              <  1:266511/78‑1 (MQ=255)
                                                                                                         cccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACaa                                                                                                              >  2:266511/1‑78 (MQ=255)
                                                                                                          ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                            <  2:295883/143‑1 (MQ=255)
                                                                                                          ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                            >  1:295883/1‑143 (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                                                                                                       >  1:157525/1‑65 (MQ=255)
                                                                                                                             ttCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTccc                                                                                                       <  2:157525/65‑1 (MQ=255)
                                                                                                                                    cTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGt            >  8:92848/1‑149 (MQ=255)
                                                                                                                                       gtgCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTg           <  4:159909/147‑1 (MQ=255)
                                                                                                                                              tCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCAt  >  3:284813/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGGGAAGTCAT  >  NC_002947/2063079‑2063369

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: