Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,315,630 G→T 20.6% A191E (GCG→GAG)  PP_2915 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,315,6300GT20.6% 75.1 / 11.1 34A191E (GCG→GAG) PP_2915membrane protein
Reads supporting (aligned to +/- strand):  ref base G (16/11);  new base T (2/5);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

ACCCGCGCATCGAATGATCAGCGCGGCATGTAGCCACGCTGCCAGGACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCAC  >  NC_002947/3315483‑3315775
                                                                                                                                                   |                                                                                                                                                 
aCCCGCGCATCGAATGATCAGCGCGGCATGTAGCCACGCTGCCAGGACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGcgcg                                                                                                                                                 >  2:763367/1‑150 (MQ=255)
 cccGCGCATCGAATGATCAGCGCGGCATGTAGCCACGCTGCCAGGACCGCGGACTGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCgcgg                                                                                                                                                <  1:763367/150‑1 (MQ=255)
  ccGCGCATCGAATGATCAGCGCGGCATGTAGCCACGCTGCCAGGACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGCGCtcggc                                                                                                                                               <  1:200419/150‑1 (MQ=255)
       cATCGAATGATCAGCGCGGCATGTAGCCACGCTGCCAGGACCGCGGAATGAACCGCGAGATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGCGCTCGGCGcctg                                                                                                                                          <  1:350646/150‑3 (MQ=255)
        aTCGAATGATCCGCGCGGCATGTAGCCACGCTGCCAGGACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGCGCTCGGCGcctgg                                                                                                                                         <  1:745333/150‑4 (MQ=255)
        aTCGAATGATCAGCGCGGCATGTAGCCACGCTGCCAGGACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGTGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAg                                                                                                                                         <  2:592393/150‑1 (MQ=255)
        aTCGAATGATCAGCGCGGCATGTAGCCACGCTGCCAGGACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAg                                                                                                                                         >  2:756632/1‑150 (MQ=255)
          cGAATGATCAGCGCGGCATGTAGCCACGCTGCCAGGACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCGGGCGCCGCGCTCGGCGcctggtt                                                                                                                                       <  1:374991/150‑6 (MQ=255)
                        ggCATGTAGCCACGCTGCCAGGACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCaa                                                                                                                         <  1:756632/150‑1 (MQ=255)
                                gCCACGCTGCCAGGACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGCGCTCGGCGcctggtt                                                                                                                                       >  2:759205/1‑123 (MQ=255)
                                gCCACGCTGCCAGGACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGCGCTCGGCGcctggtt                                                                                                                                       <  1:759205/128‑6 (MQ=255)
                                gCCACGCTGCCAGGACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCaa                                                                                                                 >  2:254691/1‑150 (MQ=255)
                                           aGGACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGcc                                                                                                      >  2:439099/1‑150 (MQ=255)
                                            ggACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCa                                                                                                     >  2:770572/1‑150 (MQ=255)
                                               ccGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGcc                                                                                                  >  1:470678/1‑150 (MQ=255)
                                                cgcgGCAGGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGCCCCAAGACACCAACTGCACCCGCCAGCCg                                                                                                 >  2:211616/1‑150 (MQ=255)
                                                         aaCCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGcgccgc                                                                                        <  1:527624/150‑1 (MQ=255)
                                                            cgcgATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGCGCTCGGCGcctgg                                                                                                                                         >  2:657520/1‑95 (MQ=255)
                                                            cgcgATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGCGCTCGGCGcctgg                                                                                                                                         <  1:657520/98‑4 (MQ=255)
                                                                   tCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCcagca                                                                              >  1:260767/1‑150 (MQ=255)
                                                                     aGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTCGCGCCGCTCGCCCagcagc                                                                            <  1:254691/150‑1 (MQ=255)
                                                                         ccAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAAc                                                                        >  1:603043/1‑150 (MQ=255)
                                                                               cGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCa                                                                  <  2:473454/150‑1 (MQ=255)
                                                                                gCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAg                                                                 >  1:326664/1‑150 (MQ=255)
                                                                                     ccAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCcgc                                                            <  2:470678/150‑1 (MQ=255)
                                                                                               gCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAAcg                                                  <  2:16197/150‑1 (MQ=255)
                                                                                                          gcgcCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGgcgc                                       >  2:737636/1‑150 (MQ=255)
                                                                                                            gcCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGccc                                     <  2:326664/150‑1 (MQ=255)
                                                                                                                      ccAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTaa                            >  2:287061/1‑149 (MQ=255)
                                                                                                                      ccAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAg                           >  1:637981/1‑150 (MQ=255)
                                                                                                                      ccAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCACAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTaa                            <  1:287061/149‑1 (MQ=255)
                                                                                                                          gtgtTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAg                       >  2:196307/1‑150 (MQ=255)
                                                                                                                               gttgCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGgc                                         >  2:63864/1‑127 (MQ=255)
                                                                                                                               gttgCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGgc                                         <  1:63864/127‑1 (MQ=255)
                                                                                                                               gttgCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGa                  >  1:173232/1‑150 (MQ=255)
                                                                                                                                              cGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCa   <  1:737636/150‑1 (MQ=255)
                                                                                                                                               gAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCAc  >  2:728218/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ACCCGCGCATCGAATGATCAGCGCGGCATGTAGCCACGCTGCCAGGACCGCGGAATGAACCGCGATATCAGCGCCAGCACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCAC  >  NC_002947/3315483‑3315775

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