Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,381,431 G→A 7.4% G26D (GGC→GAC)  PP_2983 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,381,4310GA 7.4% 157.9 / 2.4 54G26D (GGC→GAC) PP_2983hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/25);  new base A (2/2);  total (27/27)
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.68e-01

TCAGGGCACCACTTGCGAACAGCGATGGCGCACCGAGCACATGAACCGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCACCTGCTCACCCATGCGCGTGGTAGTGCGGACATCGATGCG  >  NC_002947/3381284‑3381570
                                                                                                                                                   |                                                                                                                                           
tCAGGGCACCACTTGCGAACAGCGATGGCGCACCGAGCACATGAACCGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTaggc                                                                                                                                            >  8:48759/1‑149 (MQ=255)
tCAGGGCACCACTTGCGAACAGCGATGGCGCACCGAGCACATGAACCGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTaggc                                                                                                                                            >  2:42905/1‑149 (MQ=255)
          acTTGCGAACAGCGATGGCGCACCGAGCACATGAACCGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTg                                                                                                                                  <  3:356700/149‑1 (MQ=255)
            ttGCGAACAGCGATGGCGCACCGAGCACATGAACCGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTgg                                                                                                                                 >  1:191871/1‑148 (MQ=255)
            ttGCGAACAGCGATGGCGCACCGAGCACATGAACCGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTgg                                                                                                                                 >  4:167856/1‑148 (MQ=255)
               cGAACAGCGATGGCGCACCGAGCACATGAACCGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAAt                                                                                                                             <  2:279820/149‑1 (MQ=255)
                    aggtaTGGCGCACCGAGCACATGAACCGGCAAAATTGAGTCCATGTGATGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGTGCGGCGGAACCTGTCCGTTTGACTACACCGCTATGCCAGGCCATTGGTAATCTGtt                                                                                                                        <  8:197317/145‑1 (MQ=255)
                      cGATGGCGCACCGAGCACATGAACCGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCAAGGCCAGGCCATTGGTAATCTGTTCg                                                                                                                      >  1:373468/1‑149 (MQ=255)
                             gcACCGAGCACATGAACCGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTc                                                                                                                       <  6:72257/141‑1 (MQ=255)
                             gcACCGAGCACATGAACCGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTc                                                                                                                       >  5:72257/1‑141 (MQ=255)
                                   aGCACATGAACCGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGcc                                                                                                         <  3:243618/149‑1 (MQ=255)
                                            aCCGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGa                                                                                                >  5:327302/1‑149 (MQ=255)
                                              cGGCACACATGAGTGCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGTGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATg                                                                                              <  2:224322/149‑1 (MQ=255)
                                              cGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATg                                                                                              <  2:191871/149‑1 (MQ=255)
                                              cGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATg                                                                                              <  2:34098/149‑1 (MQ=255)
                                                    acaTGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAAcc                                                                                        <  4:241634/149‑1 (MQ=255)
                                                    acaTGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAAcc                                                                                        >  3:241634/1‑149 (MQ=255)
                                                       tGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTc                                                                                     <  6:327302/149‑1 (MQ=255)
                                                         aGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTc                                                                                     >  2:106068/1‑147 (MQ=255)
                                                            ccATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGgc                                                                                >  3:191923/1‑149 (MQ=255)
                                                             cATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGgct                                                                               >  1:9102/1‑149 (MQ=255)
                                                                    tGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTAGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGgctcgctc                                                                          >  7:260266/1‑147 (MQ=255)
                                                                      aGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAg                                                                       <  4:191923/148‑1 (MQ=255)
                                                                      aGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCAAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAg                                                                       <  2:373468/148‑1 (MQ=255)
                                                                          gAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGgtctt                                                                  >  5:322419/1‑149 (MQ=255)
                                                                          gAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGgtct                                                                   <  2:72338/148‑1 (MQ=255)
                                                                             gTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGgtcttcgt                                                               <  3:2451/149‑1 (MQ=255)
                                                                               gAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGgtcttcgtct                                                             >  7:298796/1‑149 (MQ=255)
                                                                                aaGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGgtcttcgtct                                                             <  8:298796/148‑1 (MQ=255)
                                                                                aaGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTATGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCAGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGgtcttagtct                                                             >  1:89444/1‑148 (MQ=255)
                                                                                    tGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTcggc                                                        <  3:329422/149‑1 (MQ=255)
                                                                                      acacTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGgtcttc                                                                 <  7:54040/138‑1 (MQ=255)
                                                                                      acacTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGgtcttc                                                                 >  8:54040/1‑138 (MQ=255)
                                                                                        acTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGgtc                                                                    >  8:224821/1‑133 (MQ=255)
                                                                                        acTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGgtc                                                                    <  7:224821/133‑1 (MQ=255)
                                                                                          tAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGACCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTgc                                                  >  4:47547/1‑149 (MQ=255)
                                                                                          tAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGACCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTgc                                                  >  4:52418/1‑149 (MQ=255)
                                                                                            aTTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTgc                                                  <  1:106068/147‑1 (MQ=255)
                                                                                                agaCTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTg                                            <  7:33294/149‑1 (MQ=255)
                                                                                                agaCTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGACCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGt                                             <  3:47547/148‑1 (MQ=255)
                                                                                                agaCTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGACCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGt                                             <  3:52418/148‑1 (MQ=255)
                                                                                                  aCTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCa                                          >  6:14580/1‑149 (MQ=255)
                                                                                                  aCTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCa                                          >  1:70159/1‑149 (MQ=255)
                                                                                                   cTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCAc                                         <  3:167856/149‑1 (MQ=255)
                                                                                                    tGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCAcc                                        >  4:260848/1‑149 (MQ=255)
                                                                                                       cAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCACCTGc                                     >  3:14841/1‑149 (MQ=255)
                                                                                                                ggCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTACGGCGGCTGCGCGGTGCACCTGCTCACCCATg                            >  8:213527/1‑149 (MQ=255)
                                                                                                                            tGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCACCTGCTCACCCATGCGCGTGGTAGTg                >  2:336179/1‑149 (MQ=255)
                                                                                                                                 gTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCACCTGCTCACCCATGCGCGTGGTAGTGCGGa            >  7:177214/1‑148 (MQ=255)
                                                                                                                                    tGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCACCTGCTCACCCATGCGCGTGGTAGTGCGGACAt         <  1:336179/148‑1 (MQ=255)
                                                                                                                                     gACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCACCTGCTCACCCATGCGCGTGGTAGTGCGGACATc        >  2:43755/1‑148 (MQ=255)
                                                                                                                                       cTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCACCTGCTCACCCATGCGCGTGGTAGTGCGGACATCGAt     <  2:9102/149‑1 (MQ=255)
                                                                                                                                       cTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCACCTGCTCACCCATGCGCGTGGTAGTGCGGACATCGAt     <  7:48759/149‑1 (MQ=255)
                                                                                                                                        tACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTg                                            <  4:182824/109‑1 (MQ=255)
                                                                                                                                        tACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTg                                            >  3:182824/1‑109 (MQ=255)
                                                                                                                                             cGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCACCTGCTCACCCATGCGCGTGGTAGTGCGGACATCGATGCg  >  3:1537/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
TCAGGGCACCACTTGCGAACAGCGATGGCGCACCGAGCACATGAACCGGCACACATGAGTCCATGTGCTGAGTGGAGGTGAAGATGACACTAATTGAGACTGCCAATGAGGGGGCCGCGGAACCTGTCCGTTTGACTACACCGCTAGGCCAGGCCATTGGTAATCTGTTCGCCCGGGCCTTGCCGCTACTGGATGGCAACCCTCCAGGCTCGCTCAAGGTCTTCGTCTTCGGCGGCTGCGCGGTGCACCTGCTCACCCATGCGCGTGGTAGTGCGGACATCGATGCG  >  NC_002947/3381284‑3381570

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