Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,063,220 G→T 13.2% intergenic (+77/‑130) 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,2200GT13.2% 98.0 / 4.4 36intergenic (+77/‑130)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/18);  new base T (2/3);  total (15/21)
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

GTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGG  >  NC_002947/2063074‑2063361
                                                                                                                                                   |                                                                                                                                             
gTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTc                                                                                                                                             <  3:31944/149‑1 (MQ=255)
     tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgtg                                                                                                                                        <  7:91199/149‑1 (MQ=255)
     tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgtg                                                                                                                                        <  5:112049/149‑1 (MQ=255)
     tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                                                         <  5:92568/148‑1 (MQ=255)
          aCCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGc                                                                                                                                   <  8:55391/149‑1 (MQ=255)
             gcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCggg                                                                                                                                <  2:248532/149‑1 (MQ=255)
                aCCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGCCATTTCACAGGGTTCGGACCTTGTGCAGTTCCTGTGCGAGCGGGCaa                                                                                                                             <  5:130367/149‑1 (MQ=255)
                      agCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCcgc                                                                                                                       <  3:311305/149‑1 (MQ=255)
                             tgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATggg                                                                                                                >  5:329609/1‑149 (MQ=255)
                               gATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCAATCGCGGGCTTGCCCGCTCCCAGAGGCATTGCACAGGGTTCGGTCCTTGTGCAGGTGCTGTGCGAGCGGGCAAGCCCGCGAATggggc                                                                                                              >  1:283802/1‑147 (MQ=255)
                                aTAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTTCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCCCGAATCGGccc                                                                                                             >  6:113363/1‑149 (MQ=255)
                                  aaCACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCg                                                                                                           >  1:82222/1‑149 (MQ=255)
                                      ccgccgGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCCCAGGGGTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCCAGCCCGCGCATGGGCCCCGcac                                                                                                        >  7:327321/1‑148 (MQ=255)
                                                     ccGGAACACCTAAGTAACCGGGGGCCGCTTGGGGGCCCTTGCGGGGGCTTGCCCGCTCCCACTGGCTTTCCACTGGGTTCGCACCTTGTGCTGTCCCTGTGCTCGCGGGAATGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGa                                                                                        <  4:133039/148‑1 (MQ=255)
                                                             ccGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGGGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATc                                                                                <  1:50686/149‑1 (MQ=255)
                                                                 aGTACACCGGGGCCGCGTTGCGGCCCATGCGCGGGCGTGCCCGCGCCCCCAGGCATTGCACCGGGTTCGGCCCTTGTGCAGGGCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCg                                                                            <  8:327321/149‑1 (MQ=38)
                                                                  gTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGCGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGc                                                                           >  5:330297/1‑149 (MQ=255)
                                                                     aaCCGGGGCCGCGTTGCGGCCCATTCGCGGGCATTCCCCCGCCCCCAGGCCATGGACAGGGGTAGGCCCATTTTGAGGTACAGTTCGCGCCGGCCAGCCCGCCCATTGGCCCCCCACAAATTTAACCGACGGAAAATTTTAACCGccc                                                                         >  4:101156/1‑148 (MQ=255)
                                                                             ccGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTTTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAg                                                                <  5:113363/149‑1 (MQ=255)
                                                                               gcgTTGCGGCCCATTCGCGGGCTTGCCCGCGCCCACAGGCATTGCACAGGGTTCGGGGCCTTGTGCAGGGCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTGCCGGACTGAACATTATGCCCGCCCCAAgggccgg                                                               >  1:98742/1‑144 (MQ=37)
                                                                                  ttGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTAtt                                                           <  8:72986/149‑1 (MQ=255)
                                                                                  ttGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTAtt                                                           <  7:274036/149‑1 (MQ=255)
                                                                                  ttGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTAtt                                                           <  2:82222/149‑1 (MQ=255)
                                                                                         ccATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCGCCACAATTTTACCGGACTGAACATTATAACCCCCCCAAGGTCAGTTATTCAGGCa                                                     >  6:292081/1‑148 (MQ=37)
                                                                                           aTTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATca                                                  >  7:19313/1‑149 (MQ=255)
                                                                                                   gCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTg                                          <  8:19313/149‑1 (MQ=255)
                                                                                                         ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                    >  2:338255/1‑149 (MQ=18)
                                                                                                         ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                                    >  5:160940/1‑149 (MQ=255)
                                                                                                            cccccACCGGCATTGCACCGGTCTCGCCCTTCGTGCCGGTCCTGGGCGCGCGGGCATGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt                                 <  2:283802/147‑1 (MQ=255)
                                                                                                               ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGAc                              >  3:334420/1‑149 (MQ=255)
                                                                                                                cacaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCACCGCACAATTTTCACGGAATGAACATTATCCCCGACCCAAGGTCAGTTATTAAGGCATCACACACCTGGTGTTAGTTGaat                             >  5:613/1‑147 (MQ=255)
                                                                                                                     gacctGCACAAGGTCCGAACCCTGTGCAATGCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCAGGACTGAACATTATCCCCGCCCCAAGTTCAGTTATTCAGGCATCACACACCTGGTTTTCGTTTACAgatcg                        >  4:10275/5‑146 (MQ=21)
                                                                                                                     gaactGCACAAGGTCCGAACCCTGTGCAATGCCTGTGCGAGCTGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGa                               <  3:10275/138‑1 (MQ=25)
                                                                                                                                 gTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt            <  6:329609/149‑1 (MQ=255)
                                                                                                                                             gtgCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTgg  <  6:330297/148‑1 (MQ=255)
                                                                                                                                             gtgCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTg   <  3:101156/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGTGG  >  NC_002947/2063074‑2063361

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: