Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,315,627 A→C 21.7% L192R (CTC→CGC)  PP_2915 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,315,6270AC21.7% 60.7 / 17.6 37L192R (CTC→CGC) PP_2915membrane protein
Reads supporting (aligned to +/- strand):  ref base A (17/12);  new base C (2/6);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-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: