Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,271,710 A→G 8.1% intergenic (+51/‑49) PP_t32 → / → PP_t33 tRNA‑Asp/tRNA‑Asp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,271,7100AG 8.1% 142.3 / 3.0 49intergenic (+51/‑49)PP_t32/PP_t33tRNA‑Asp/tRNA‑Asp
Reads supporting (aligned to +/- strand):  ref base A (21/24);  new base G (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

AGCCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCTAG  >  NC_002947/2271562‑2271856
                                                                                                                                                    |                                                                                                                                                  
aGCCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACAt                                                                                                                                                   <  2:1085365/150‑1 (MQ=255)
       gagGTACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTgg                                                                                                                                            <  1:810083/150‑1 (MQ=255)
       gagGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTgg                                                                                                                                            <  2:1414600/150‑1 (MQ=255)
             cGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGAt                                                                                                                                      <  2:1391317/150‑1 (MQ=255)
                        gCGGTAGTTCAGTCGGTTACAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCNCCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCaa                                                                                                                           <  1:1102365/150‑1 (MQ=255)
                                        ttAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACAt                                                                                                                                                   <  2:1044505/110‑1 (MQ=255)
                                        ttAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACAt                                                                                                                                                   >  1:1044505/1‑110 (MQ=255)
                                         tAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                          <  1:1255725/150‑1 (MQ=255)
                                          cgacccccGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTATCAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAAc                                                                                                         <  1:1361467/144‑1 (MQ=255)
                                                 cGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGc                                                                                                  >  2:231785/1‑150 (MQ=255)
                                                  ggCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATCGGAGCCTGGCACCGGTCTGAAGCCACAAAGCACATCCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATAAGCa                                                                                                 >  2:1554341/1‑150 (MQ=38)
                                                   gCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGc                                                                                                                                          <  1:1560434/108‑1 (MQ=255)
                                                   gCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGc                                                                                                                                          >  2:1560434/1‑108 (MQ=255)
                                                   gCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAg                                                                                                <  2:1358432/150‑1 (MQ=255)
                                                        tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCTCCGGTCAGAAACAACAAGGCACATGCCTTTGGCTGATCTCAGGTTCAATCGGTCACACGCCTGCACGATACGCAGCGGTa                                                                                           >  2:1618103/1‑150 (MQ=37)
                                                        tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGCTCTGAAGCC‑CAAGGCATGTGCCTTgtggc                                                                                                                                         <  1:374773/103‑5 (MQ=25)
                                                        tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGCTCTGAAGCC‑CAAGGCATGTGCCTTgtggc                                                                                                                                         >  2:374773/1‑99 (MQ=25)
                                                           cGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGtt                                                                                        <  2:1338387/150‑1 (MQ=255)
                                                                ggggTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAggc                                                                                   <  1:231785/150‑3 (MQ=255)
                                                                ggggTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTGCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTc                                                                                   <  1:827673/150‑1 (MQ=255)
                                                                              gAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGcaca                                                                                                                                                    >  1:209082/1‑71 (MQ=255)
                                                                              gAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGcaca                                                                                                                                                    <  2:209082/71‑1 (MQ=255)
                                                                                   ccGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTgg                                                                                                                                            >  2:1061545/1‑74 (MQ=255)
                                                                                   ccGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTgg                                                                                                                                            <  1:1061545/74‑1 (MQ=255)
                                                                                   ccGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTg                                                                                                                                             <  2:482930/73‑1 (MQ=255)
                                                                                   ccGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTg                                                                                                                                             >  1:482930/1‑73 (MQ=255)
                                                                                    cGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTc                                                               >  1:1560534/1‑150 (MQ=255)
                                                                                      tCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGgc                                                                                                                     >  2:399991/1‑94 (MQ=255)
                                                                                      tCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGgc                                                                                                                     <  1:399991/94‑1 (MQ=255)
                                                                                      tCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCa                                                              >  2:482414/1‑149 (MQ=255)
                                                                                      tCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCa                                                              <  1:482414/149‑1 (MQ=255)
                                                                                        cGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGc                                                           <  2:1560534/150‑1 (MQ=255)
                                                                                             cgcCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCATGTGCCTTgtg                                                                                                                                           >  1:297810/1‑63 (MQ=255)
                                                                                             cgcCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCATGTGCCTTgtg                                                                                                                                           <  2:297810/65‑3 (MQ=255)
                                                                                                         gTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCg                                          >  1:624506/1‑150 (MQ=255)
                                                                                                         gTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCg                                          >  2:1387346/1‑150 (MQ=255)
                                                                                                          tAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGa                                         <  2:624506/150‑1 (MQ=255)
                                                                                                              aGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCACAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTcc                                     >  2:703011/1‑150 (MQ=255)
                                                                                                                    gTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCg                               >  2:479759/1‑150 (MQ=255)
                                                                                                                           gCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGcag                                                                                             <  2:917762/81‑3 (MQ=255)
                                                                                                                           gCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGcag                                                                                             >  1:917762/1‑79 (MQ=255)
                                                                                                                           gCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGcc                        >  2:574013/1‑150 (MQ=255)
                                                                                                                               cGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACAGTACGCAGCCGCGGCGCACTCGGATAGAATACCAGCCTGTCACTCCGGGCCGCGCGGCCCCGTGTCCCGTCCCCAACGCCATCt                    >  2:124223/1‑150 (MQ=16)
                                                                                                                                ggTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCtt                   <  1:1387346/150‑1 (MQ=255)
                                                                                                                                    tGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATa                                                                                                     <  1:912933/64‑1 (MQ=255)
                                                                                                                                    tGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATa                                                                                                     >  2:912933/1‑64 (MQ=255)
                                                                                                                                      aaGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATAc                                                                        <  2:1339084/91‑1 (MQ=255)
                                                                                                                                      aaGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATAc                                                                        >  1:1339084/1‑91 (MQ=255)
                                                                                                                                           acaAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGAt                                                                                                      >  2:570749/1‑56 (MQ=255)
                                                                                                                                           acaAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGAt                                                                                                      <  1:570749/56‑1 (MQ=255)
                                                                                                                                           acaAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGt        >  2:877498/1‑150 (MQ=255)
                                                                                                                                                 cacaTGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCTAg  <  1:877498/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGCCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCCTTGGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCCTAG  >  NC_002947/2271562‑2271856

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑