Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,107,357 Δ1 bp 8.5% coding (40/429 nt) PP_3614 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,107,3570T. 8.5% 138.3 / 8.3 47coding (40/429 nt)PP_3614membrane protein
Reads supporting (aligned to +/- strand):  ref base T (20/22);  new base . (2/2);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CCCACAGAAGACCAGGGCCCCGGCGATGAGCCAGGCGGCGAAGTTGCTCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGCCCCTTGACGCCCCGATAGT  >  NC_002947/4107211‑4107497
                                                                                                                                                   |                                                                                                                                            
cccACAGAAGACCAGGGCCCCGGCGATGAGCCAGGCGGCGAAGTTGCTCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTgg                                                                                                                                            <  3:236539‑M2/149‑3 (MQ=255)
  cacaGAAGACCAGGGCCCCGGCGATGGGCCCGGCGGCGAAGTTGCGCAACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCCGCAGCAGGGCATGCAGCGGGCGgggc                                                                                                                                          <  1:10415‑M2/149‑5 (MQ=255)
   acaGAAGACCAGGGCCCCGGCGATGAGCCAGGCGGCGAAGTTGCTCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCg                                                                                                                                         >  8:219821/1‑149 (MQ=255)
            cAGGGCCCCGGCGATGAGCCAGGCGGCGAAGTTGCTCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCgg                                                                                                                                 >  5:123026/1‑148 (MQ=255)
                  cccGGCGATGAGCCAGGCGGCGAAGTTGCTCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGATGGGCGTACAGCGGTAGAggg                                                                                                                          >  7:224292/1‑149 (MQ=255)
                   ccGGCGATGAGCCAGGCGGCGAAGTTGCTCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGt                                                                                                                         >  2:265715/1‑149 (MQ=255)
                    cGGCGATGAGCCAGGCGGCGAAGTTGCTCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGtt                                                                                                                        <  4:236539/149‑1 (MQ=255)
                       cGATGAGCCAGGCGGCGAAGTTGCTCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGa                                                                                                                     <  1:265715/149‑1 (MQ=255)
                              ccAGGCGGCGAAGTTGCTCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTg                                                                                                              <  2:10415/149‑1 (MQ=255)
                                     gcgAAGTTGCTCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCAc                                                                                                       <  3:135463/149‑1 (MQ=255)
                                     gcgAAGTTGCTCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCAc                                                                                                       <  2:49537/149‑1 (MQ=255)
                                       gAAGTTGCTCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACgg                                                                                                     >  8:284186/1‑149 (MQ=255)
                                               tCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCCAGCCCGCTGCATGCCCTGCTGCTGGCggcagtgtgccgttgttcct                                                                                                                           <  2:339660‑M2/121‑21 (MQ=255)
                                               tCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCCAGCCCGCTGCATGCCCTGCTGCTGGCggcagtgtgccgttgttcct                                                                                                                           >  1:339660‑M2/1‑101 (MQ=255)
                                                 cACTGAATCTGATAGCTGTTGTAGTAGACGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGcctcct                                                                                             <  5:257810/147‑1 (MQ=255)
                                                      aaTCTGATAGCTGTTGTAGTAGGCGATATAGCTCAGCAGCGAGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGTTCACGGTGCCTCCTTttcatt                                                                                      <  6:361473/149‑1 (MQ=255)
                                                           gATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCa                                                                                                                  >  3:196266/1‑116 (MQ=255)
                                                           gATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCa                                                                                                                  <  4:196266/116‑1 (MQ=255)
                                                            aTAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTttcattca                                                                                    <  6:123026/145‑1 (MQ=255)
                                                                   ttgtAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAggg                                                                         <  7:219821/149‑1 (MQ=255)
                                                                   ttgtAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAgg                                                                          <  7:280621/148‑1 (MQ=255)
                                                                         tagGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGt                                                                                                    >  1:62836/1‑116 (MQ=255)
                                                                         tagGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGt                                                                                                    <  2:62836/116‑1 (MQ=255)
                                                                                    gCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTgcgc                                                         >  1:43041/1‑148 (MQ=255)
                                                                                         gcagcGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAg                                                   <  7:284186/149‑1 (MQ=255)
                                                                                           agcGCGCCAAGGAACATCGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCCGGTGCGGTCACGGTGCCTGCTTTTCATTCAGTGACGCAGTGGCTGATGGTCTTGCGCCATCAGGt                                                 <  5:358355/149‑1 (MQ=255)
                                                                                                  cAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAgcg                                            >  6:82415/1‑147 (MQ=255)
                                                                                                   aaGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATg                                                                  <  5:202611/124‑1 (MQ=255)
                                                                                                   aaGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATg                                                                  >  6:202611/1‑124 (MQ=255)
                                                                                                   aaGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAg                                         >  3:174213/1‑149 (MQ=255)
                                                                                                         caacGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTgcg                                   >  2:228654/1‑149 (MQ=255)
                                                                                                             ggCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTGCGCTTa                               <  2:43041/149‑1 (MQ=255)
                                                                                                               cacaCTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTgggcagtgtgccgttgttccttggcgcgctgctgagcgatatcgcctactacaacagctatcagattcagtggagcaacttcgccgcctggctcatcgccggggccctggtct                             <  4:105239‑M2/149‑114 (MQ=255)
                                                                                                                  acTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTGCGCTTAGCa                            >  5:330195/1‑147 (MQ=255)
                                                                                                                      ccc‑GCCAGCAGCAGGGCATGCAGCGGGCgtgtgccgttgttccttggcgcgctgctgagcgatatcgcctactacaacagctatcagattcagtggagcaacttcgc                                                                <  5:288230‑M2/108‑80 (MQ=255)
                                                                                                                      ccc‑GCCAGCAGCAGGGCATGCAGCGGGCgtgtgccgttgttccttggcgcgctgctgagcgatatcgcctactacaacagctatcagattcagtggagcaacttcgc                                                                >  6:288230‑M2/1‑29 (MQ=255)
                                                                                                                       cc‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGc                     >  8:185410/1‑149 (MQ=255)
                                                                                                                       cc‑GCCAGCAGCAGGGCATGCAGCGGCCTGGGCGTACACCGGTAGAGGGATTGATCAGCTGCGGTCACCGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGc                     >  1:278332/1‑149 (MQ=255)
                                                                                                                          gCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGccc                   <  5:82415/149‑1 (MQ=255)
                                                                                                                             agcagcagGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCACCCTGCGCTTac                              <  2:306761/135‑2 (MQ=255)
                                                                                                                             agcagcagGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCACCCTGCGCTTac                              >  1:306761/1‑134 (MQ=255)
                                                                                                                                    gggCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGCCCCTTGACGccc         >  4:279619/1‑149 (MQ=255)
                                                                                                                                        aTGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGCCCCTTGACGCCCCGa      >  8:144167/1‑148 (MQ=255)
                                                                                                                                           cAGCGGGATGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCATCCTTTTCATTTAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGCCCCTTGACGCACCGATAGt  >  1:12041/1‑149 (MQ=255)
                                                                                                                                            aGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTGCGCTTAg                              <  2:296149/120‑1 (MQ=255)
                                                                                                                                                  cTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTGCGCtt                                >  6:243129/1‑112 (MQ=255)
                                                                                                                                                  cTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTGCGCtt                                <  5:243129/112‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
CCCACAGAAGACCAGGGCCCCGGCGATGAGCCAGGCGGCGAAGTTGCTCCACTGAATCTGATAGCTGTTGTAGTAGGCGATATCGCTCAGCAGCGCGCCAAGGAACAACGGCACACTGCCC‑GCCAGCAGCAGGGCATGCAGCGGGCTGGGCGTACAGCGGTAGAGGGTTTGATCAGGTGCGGTCACGGTGCCTCCTTTTCATTCAGTGACCCAGGGGCTGATGGTCTTGCGCCATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGCCCCTTGACGCCCCGATAGT  >  NC_002947/4107211‑4107497

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