Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,237,665 G→A 14.9% A148A (GCC→GCT)  PP_1078 ← ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,237,6650GA14.9% 77.4 / 5.0 34A148A (GCC→GCT) PP_1078ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (16/12);  new base A (3/2);  total (20/14)
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

CTGGTCATGGGTGACCAGAATGGCGCTGGTACCGCGGCTTTTGAGAATGTCGCGAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCTTCGATGACGGCCGCCTGGCGTGGGTGTTTGGCGATACCAAAGGCAATG  >  NC_002947/1237518‑1237803
                                                                                                                                                   |                                                                                                                                          
cTGGTCATGGGTGACCAGAATGGCGCTGGTACCGCGGCTTTTGAGAATGTCGCGAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCgg                                                                                                                                            <  1:17580/148‑1 (MQ=255)
 tGGTCATGGGTGACCAGAATGGCGCTGGTACCGCGGCTTTTGAGAATGTCGCGAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCggg                                                                                                                                           >  5:174694/1‑148 (MQ=255)
   gTCATGGGTGACCAGAATGGCGCTGGTACCGCGGCTTTTGAGAATGTCGCGAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCggg                                                                                                                                           <  6:174694/146‑1 (MQ=255)
           tGACCAGAATGGCGCTGGTACCGCGGCTTTTGAGAATGTCGCGAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGaca                                                                                                                                <  4:150377/149‑1 (MQ=255)
                   aTGGCGCTGGTACCGCGGCTTTTGAGAATGGCGCGAACCTCATGGCGCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTtgctg                                                                                                                        <  5:238751/149‑1 (MQ=255)
                    tGGCGCTGGTACCGCGGCTTTTGAGAATGTCGCGAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTtgctgc                                                                                                                       <  4:231027/149‑1 (MQ=255)
                                  cgGCTTTTGAGAATGTCGCGAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAg                                                                                                         >  7:4836/1‑149 (MQ=255)
                                   ggCTTTTGAGAATGTCGCGAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTCGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGc                                                                                                        <  5:31688/149‑1 (MQ=255)
                                     cTTTTGAGAATGTCGCGAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTc                                                                                                      <  5:98794/149‑1 (MQ=255)
                                            gaATGTCGCGAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCggg                                                                                               >  3:196854/1‑149 (MQ=255)
                                               tGTCGCGAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTaa                                                                                            >  1:271431/1‑149 (MQ=255)
                                                   gcgAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCCCGGGCGCTGGCGCCCGAGccgcaa                                                                                                                                     >  1:227065/1‑100 (MQ=38)
                                                   gcgAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCCCGGGCGCTGGCGCCCGAGccgcaa                                                                                                                                     <  2:227065/104‑5 (MQ=38)
                                                         cTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGc                                                                                                  >  6:307311/1‑133 (MQ=255)
                                                         cTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAAcgcccgcc                                                                                    >  3:132366/1‑147 (MQ=255)
                                                         cTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATCGCTCGTCAAGCAGCAGCAGTTTCGGCACCGGCGCCAGCGCCCGCGCCAGGGCGCCACGTTGCTGCTGCACGCCGGACAGCTCATGCTGGTAACGCCAGCCCa                                                                                  >  2:246115/1‑149 (MQ=255)
                                                                cTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCccgcc                                                                           >  5:72829/1‑149 (MQ=255)
                                                                     gcgcCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGtt                                                                      >  3:50922/1‑149 (MQ=255)
                                                                           gcgcAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCa                                                                <  4:50922/149‑1 (MQ=255)
                                                                             gcAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGc                                                              >  1:45671/1‑149 (MQ=255)
                                                                                         cGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATctct                                                  >  1:7530/1‑149 (MQ=255)
                                                                                                ggAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCTTCGATGa                                           >  4:51019/1‑149 (MQ=255)
                                                                                                  agaATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCATCGATGACg                                         >  6:126153/1‑149 (MQ=255)
                                                                                                   gaATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGAGCAGCTCCAGCATCTCTTCGATGACgg                                        >  6:322273/1‑149 (MQ=255)
                                                                                                            gTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCTTCGATGACGGCCGCCTGGc                               <  3:51019/149‑1 (MQ=255)
                                                                                                                agcagcagcAGTTGCGGCTCGGGCGCCAGCGACCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCTTCGATGACGGCCGCCTGGCGTgg                           <  4:196854/149‑1 (MQ=255)
                                                                                                                    gcagcagTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCtt                                                 <  2:90162/123‑1 (MQ=255)
                                                                                                                     cagcagTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCTTCGATGACGGCCGCCTGGCGTGGGTGtt                      >  8:113912/1‑149 (MQ=255)
                                                                                                                      agcagTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCTTCGATGACGGCCGCCTGGCGTGGGTGttt                     <  6:72829/149‑1 (MQ=255)
                                                                                                                        cagTTGCTGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCGAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCTTCGATGACGGCCGCCTGGCGTTGGTGTTTgg                   <  2:271431/149‑1 (MQ=255)
                                                                                                                             gCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGa                                                                         >  4:110579/1‑90 (MQ=255)
                                                                                                                               ggCCCGGGCGCTGGCGCCCGAGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCTTCGATGACGGCCGCCTGGCGTGGGTGTTTGGCGATAc             >  8:234507/1‑148 (MQ=255)
                                                                                                                                gcccGGGCGCTGGCGCCCGAGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCTTCGATGACGGCCGCCTGGCGTGGGTGTTTGGCGATAcc            >  4:296811/4‑148 (MQ=255)
                                                                                                                                   cGGGCGCTGGCGCCCGAGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCTTCGATGACGGCCGCCTGGCGTGGGTGTTTGGCGATACCAAAg        <  3:296811/149‑1 (MQ=255)
                                                                                                                                   cGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCTTCGATGACGGCCGCCTGGCGTGGGTGTTTGGCGATACCAAAg        >  3:48802/1‑149 (MQ=255)
                                                                                                                                      gcgcCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCTTCGATGa                                           >  4:179812/1‑111 (MQ=255)
                                                                                                                                         caagggcCCGGGCCAGGGCGACAGGTGGCTGCTGGCGGCCGGACAGCTAATGCGGGTAACGCCGGCCCAGCCGGCGGAGTTTTACCAGCTCCAGTATCTCTTCGATGACGGCCGCCTGGCGTGGGTGTTTGGCGATACCAAAGGCAATg  <  1:246115/144‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CTGGTCATGGGTGACCAGAATGGCGCTGGTACCGCGGCTTTTGAGAATGTCGCGAACCTCATGGCTCAGGCGCCGGCGCAGCTCCACGTCGAGGTTGGAGAATGGCTCGTCAAGCAGCAGCAGTTGCGGCTCGGGCGCCAGCGCCCGGGCCAGGGCGACACGTTGCTGCTGGCCGCCGGACAGCTCATGCGGGTAACGCCCGCCCAGCCCGCCGAGTTTGACCAGCTCCAGCATCTCTTCGATGACGGCCGCCTGGCGTGGGTGTTTGGCGATACCAAAGGCAATG  >  NC_002947/1237518‑1237803

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑