Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,315,648 C→A 7.8% E171* (GAG→TAG)  PP_2915 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,315,6480CA 7.8% 155.8 / 4.0 51E171* (GAG→TAG) PP_2915membrane protein
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base A (2/2);  total (26/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 = 4.26e-01

ACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGG  >  NC_002947/3315518‑3315794
                                                                                                                                  |                                                                                                                                                  
aCGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGcc                                                                                                                                                 <  1:220897/134‑1 (MQ=255)
aCGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGcc                                                                                                                                                 >  2:220897/1‑134 (MQ=255)
aCGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGccc                                                                                                                                  >  1:77383/1‑149 (MQ=255)
 cGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGcccc                                                                                                                                 <  1:72985/149‑1 (MQ=255)
    aGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAgg                                                                                                                              >  4:124833/1‑149 (MQ=255)
      gCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGcc                                                                                                                            >  2:225667/1‑149 (MQ=255)
                   tGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCACCGCCCCAGGCCGCGGCCGAAAcgc                                                                                                               >  3:176785/1‑149 (MQ=255)
                     ccGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGa                                                                                                             >  3:62361/1‑149 (MQ=255)
                       ggCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGa                                                                                                             >  6:39331/1‑147 (MQ=255)
                       ggCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACa                                                                                                           <  6:111847/149‑1 (MQ=255)
                            gcgcCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGgcg                                                                                                      <  6:83178/149‑1 (MQ=255)
                            gcgcCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGgcg                                                                                                      <  2:77383/149‑1 (MQ=255)
                              gcCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGcc                                                                                                    <  5:39331/149‑1 (MQ=255)
                                cTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCCGGCCGCGGCCGAAACGCGATAGGGCGCCCg                                                                                                  >  6:16917/1‑149 (MQ=255)
                                cTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGc                                                                                                                       <  3:245664/128‑1 (MQ=255)
                                cTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGc                                                                                                                       >  4:245664/1‑128 (MQ=255)
                                   cAGGCCCAGGTGTAGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGt                                                                                                 <  8:276362/147‑1 (MQ=255)
                                          aGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGt                                                                                        <  8:108417/149‑1 (MQ=255)
                                                    gCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAATCAGGCCGATg                                                                              >  1:260745/1‑149 (MQ=255)
                                                       gctgggCGCCGCGCTCGGCGCCTAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTaa                                                                                          <  2:240092/130‑1 (MQ=255)
                                                       gctgggCGCCGCGCTCGGCGCCTAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTaa                                                                                          >  1:240092/5‑134 (MQ=255)
                                                          aGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAgg                                                                        <  7:29588/149‑1 (MQ=255)
                                                              gcCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCAt                                                                    <  3:124833/149‑1 (MQ=255)
                                                               ccGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCcgcg                                                                                                                         <  8:281640/95‑1 (MQ=255)
                                                               ccGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCcgcg                                                                                                                         >  7:281640/1‑95 (MQ=255)
                                                                 gAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGAt                                                                               >  5:141502/1‑135 (MQ=255)
                                                                 gAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGAt                                                                               <  6:141502/135‑1 (MQ=255)
                                                                 gAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCa                                                                 >  3:47206/1‑149 (MQ=255)
                                                                   gcgcGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCAcc                                                               <  3:66190/149‑1 (MQ=255)
                                                                    cgcgGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCt                                                              <  2:260745/149‑1 (MQ=255)
                                                                          gcCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCa                                                        >  6:300342/1‑149 (MQ=255)
                                                                           cccAGCATCATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCAc                                                       <  5:300342/149‑1 (MQ=255)
                                                                            ccAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTAGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAGACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCa                                                        >  6:235878/1‑147 (MQ=255)
                                                                            ccAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTAGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAGACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCa                                                        >  6:231855/1‑147 (MQ=255)
                                                                            ccAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTAGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAGACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCa                                                        <  5:231855/147‑1 (MQ=255)
                                                                            ccAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTAGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAGACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCa                                                        <  5:235878/147‑1 (MQ=255)
                                                                              aGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGAc                                                    >  1:272773/1‑149 (MQ=255)
                                                                                           cccAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATg                                       >  1:186045/1‑149 (MQ=255)
                                                                                             cAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGcc                                     >  2:139790/1‑149 (MQ=255)
                                                                                                 gCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATg                                 >  5:100923/1‑149 (MQ=255)
                                                                                                           cACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCAAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTgg                       >  2:258112/1‑149 (MQ=255)
                                                                                                            aCCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGgc                      >  4:299877/1‑149 (MQ=255)
                                                                                                             cccGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGgcc                     <  1:225667/149‑1 (MQ=255)
                                                                                                              ccGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGgccg                    <  3:108139/149‑1 (MQ=255)
                                                                                                                gCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGgccgcc                  <  4:62361/149‑1 (MQ=255)
                                                                                                                            gccgcTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATg                                       >  4:35179/1‑116 (MQ=255)
                                                                                                                            gccgcTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATg                                       <  3:35179/116‑1 (MQ=255)
                                                                                                                             ccgcTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAggg      >  2:104674/1‑148 (MQ=255)
                                                                                                                               gcTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAg   <  2:272773/149‑1 (MQ=255)
                                                                                                                               gcTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAg   >  4:19461/1‑149 (MQ=255)
                                                                                                                                cTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAgg  <  3:19461/149‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                                  
ACGCAGGCCAGCGCACAGCTGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGG  >  NC_002947/3315518‑3315794

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