Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,777,233 Δ1 bp 13.3% coding (315/585 nt) PP_2429 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,777,2330G.13.3% 83.1 / 9.1 30coding (315/585 nt)PP_2429membrane protein
Reads supporting (aligned to +/- strand):  ref base G (11/15);  new base . (2/2);  total (13/17)
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.91e-01

GCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGCTCATCAGCTGTTGCGCCAGGGGGTAGCGCAGCAGCAGTTCGCCTAGCCCCAGCCCGA  >  NC_002947/2777092‑2777381
                                                                                                                                             |                                                                                                                                                    
gCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTtg                                                                                                                                              >  2:337201/1‑150 (MQ=255)
    ggCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGcc                                                                                                                                          >  2:268396/1‑150 (MQ=255)
       agagCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTg                                                                                                                                       <  2:777742/150‑1 (MQ=255)
                    aGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCg                                                                                                                          <  2:557996/150‑1 (MQ=255)
                          ccAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGggcgg                                                                                                                    <  1:667002/150‑1 (MQ=255)
                             gcgcCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGggcggcgg                                                                                                                 <  1:268396/150‑1 (MQ=255)
                                ccAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCAc                                                                                                              >  1:386668/1‑150 (MQ=255)
                                       ccACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCgcag                                                                                                       >  2:750746/1‑150 (MQ=255)
                                              ggCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTc                                                                                                <  2:238714/150‑1 (MQ=255)
                                                     acgGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCa                                                                                         <  2:775654/150‑1 (MQ=255)
                                                          gcgGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGcc                                                                                                                                                 <  1:49019/89‑1 (MQ=255)
                                                          gcgGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGcc                                                                                                                                                 >  2:49019/1‑89 (MQ=255)
                                                                 aTACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGcc                                                                             <  1:436730/150‑1 (MQ=255)
                                                                    ccccGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAg                                                                         >  1:841149/1‑150 (MQ=255)
                                                                    ccccGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAg                                                                         >  1:841158/1‑150 (MQ=255)
                                                                            accGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCg                                                                  <  1:491404/150‑1 (MQ=255)
                                                                                     gCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGCtca                                                        <  2:841149/150‑1 (MQ=255)
                                                                                     gCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGCtca                                                        <  2:841158/150‑1 (MQ=255)
                                                                                        aTCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGGAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGCtcatc                                                      >  1:296498/1‑150 (MQ=255)
                                                                                             ccAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGCTCATCAGCTg                                                 >  1:632022/1‑150 (MQ=255)
                                                                                               aGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGCTCATCAGCTGtt                                               <  1:477128/150‑1 (MQ=255)
                                                                                                gcccTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCCTCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGCTCATCAGCTGTTg                                              <  1:62031/148‑1 (MQ=255)
                                                                                                               caccTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGTTGATCCACAGCACCCCGGCCCAGCTCATCAGCTGTTGCGCCAGGGGGTAgcg                               >  1:692553/1‑150 (MQ=255)
                                                                                                                accTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGCTCATCAGCTGTTGCGCCAGGGGGTAgcgc                              <  1:750746/150‑1 (MQ=255)
                                                                                                                   tGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGCTCATCAGCTGTTGCGCCAGGGGGTAGCgcagc                           >  1:318170/1‑150 (MQ=255)
                                                                                                                              cGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGCTCATCAGCTGTTGCGCCAGGGGGTAGCGCAGCAGCAGTTCGcc                <  2:318170/150‑1 (MQ=255)
                                                                                                                                 aCTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGc                                                           >  1:715245/1‑104 (MQ=255)
                                                                                                                                 aCTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGc                                                           <  2:715245/104‑1 (MQ=255)
                                                                                                                                  ttGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGCTCATCAGCTGTTGCGCCAGGGGGTAGCGCAGCAGCAGTTCGCCTAGc            <  2:442026/149‑1 (MQ=255)
                                                                                                                                            tGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGCTCATCAGCTGTTGCGCCAGGGGGTAGCGCAGCAGCAGTTCGCCTAGCCCCAGCCcga  >  2:383788/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTTCACCACCTGCAACGAGGCCGCACTGAGCGCACTGAAGCCTTGTGCCGTGGCCGCATGGAGCGGGGCGGCGGCACTGCGCAGCATCTTCCACGCCAGGCTGCTCAGCCACAGCACCCCGGCCCAGCTCATCAGCTGTTGCGCCAGGGGGTAGCGCAGCAGCAGTTCGCCTAGCCCCAGCCCGA  >  NC_002947/2777092‑2777381

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