Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,097,126 G→A 8.7% A44T (GCC→ACC)  PP_2705 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,097,1260GA 8.7% 157.8 / 4.5 57A44T (GCC→ACC) PP_2705hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/24);  new base A (3/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.95e-01

GAAGCGCGTCGTGACAGGCTTCATCCCGTTATCGTCGCTCTGGCGATTCACACCGACGGACCAGGCTCGGTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGCCAGCTGTTTTATGCCGAAGTGCGCAAGAGGG  >  NC_002947/3096987‑3097273
                                                                                                                                             |                                                                                                                                                   
gAAGCGCGTCGTGACAGGCTTCATCCCGTTATCGTCGCTCTGGCGATTCACACCGACGGACCAGGCTCGGTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGCTTCAACctgctgca                                                                                                                                            <  8:137893/149‑3 (MQ=255)
        tcgtGACAGGCGTCATCCCGGTTTCGTCGCTCGTGCGATTCAAACCGACGGACCAGGCTCGGGCCCAACTCGGTGCTGCAATGGCGTCTATCGAGTACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGta                                                                                                                                     <  1:251659/149‑1 (MQ=255)
               aGGCTTCATCCCGTTATCGTCGCTCTGGCGATTCACACCGACGGACCAGGCTCGGTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGt                                                                                                                             >  1:245014/1‑149 (MQ=255)
                 gCTTCATCCCGTTATCGTCGCTCTGGCGATTCACACCGACGGACCAGGCTCGGTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAGGGGGCTGCCAGCAGGCTTCAACctgctgcagc                                                                                                                                          <  1:250601/135‑5 (MQ=255)
                 gCTTCATCCCGTTATCGTCGCTCTGGCGATTCACACCGACGGACCAGGCTCGGTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAGGGGGCTGCCAGCAGGCTTCAACctgctgcagc                                                                                                                                          >  2:250601/1‑131 (MQ=255)
                       tCCCGTTATCGTCGCTCTGGCGATTCACACCGACGGACCAGGCTCGGTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCgg                                                                                                                     >  6:108602/1‑149 (MQ=255)
                        cccGTTATCGTCGCTCTGGCGATTCACACCGACGGACCAGGCTCGGTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGCCGTATACCTCGTCACAGCGGc                                                                                                                    >  8:209392/1‑149 (MQ=255)
                                 gtcgCTCTGGCGATTCACACCGACGGACCAGGCTCGGTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGaa                                                                                                            <  5:108602/148‑1 (MQ=255)
                                    gCTCTGGCGATTCACACCGACGGACCAGGCTCGGTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGcc                                                                                                        >  8:250807/1‑149 (MQ=255)
                                          gCGATTCACACCGACGGACCAGGCTCGGTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTcgcg                                                                                                   <  7:250807/148‑1 (MQ=255)
                                                   acCGACGGACCAGGCTCGGTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAg                                                                                                           >  4:143811/1‑131 (MQ=255)
                                                   acCGACGGACCAGGCTCGGTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAg                                                                                                           <  3:143811/131‑1 (MQ=255)
                                                         ggACCAGGCTCGGTCCAAAGATGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGaa                                                                                                            <  6:85693/124‑1 (MQ=255)
                                                                    ggTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGa                                                                         >  6:142548/1‑148 (MQ=255)
                                                                     gTCCAAAGAGGGTGCTGCAATGGCGACAATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGAGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGaaa                                                                       >  1:183456/1‑149 (MQ=255)
                                                                         aaaGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGAc                                                                                          >  6:153105/1‑126 (MQ=255)
                                                                         aaaGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGAc                                                                                          <  5:153105/126‑1 (MQ=255)
                                                                              gggTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCAtt                                                                           <  6:220875/136‑1 (MQ=255)
                                                                              gggTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCAtt                                                                           >  5:220875/1‑136 (MQ=255)
                                                                               ggTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAACCCGCCAACcg                                                             >  8:306520/1‑149 (MQ=255)
                                                                               ggTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGACACCcg                                                             >  8:20519/1‑149 (MQ=255)
                                                                                 tgctgcAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACcgc                                                            >  7:308557/1‑148 (MQ=255)
                                                                                    tgcAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGa                                                          >  3:229376/1‑147 (MQ=255)
                                                                                    tgcAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGa                                                          <  4:229376/147‑1 (MQ=255)
                                                                                     gcAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATAcc                                                                                                                                 <  7:44105/75‑1 (MQ=255)
                                                                                     gcAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATAcc                                                                                                                                 >  8:44105/1‑75 (MQ=255)
                                                                                     gcAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGAttt                                                       >  7:98738/1‑149 (MQ=255)
                                                                                        aTGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCAc                                                    <  5:142548/149‑1 (MQ=255)
                                                                                             gACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGt                                               >  2:33455/1‑149 (MQ=255)
                                                                                                  tCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCg                                                                                                                                       >  2:286407/1‑56 (MQ=255)
                                                                                                  tCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCg                                                                                                                                       <  1:286407/56‑1 (MQ=255)
                                                                                                  tCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGttt                                             >  3:110897/1‑146 (MQ=255)
                                                                                                     acacCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCa                                           >  4:183717/1‑145 (MQ=255)
                                                                                                     acacCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCa                                           <  3:183717/145‑1 (MQ=255)
                                                                                                     acacCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAAc                                       <  7:209392/149‑1 (MQ=255)
                                                                                                          cGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCACCCGCGAt                                                         >  8:254181/1‑127 (MQ=255)
                                                                                                            agaACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCg                                                                                        <  1:257192/94‑1 (MQ=255)
                                                                                                            agaACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCg                                                                                        >  2:257192/1‑94 (MQ=255)
                                                                                                                  ttGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGCCAGCTg                           >  2:128349/1‑149 (MQ=255)
                                                                                                                      aaaGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGa                                                          <  5:306038/114‑1 (MQ=255)
                                                                                                                      aaaGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGa                                                          >  6:306038/1‑114 (MQ=255)
                                                                                                                      aaaGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCa                                           >  3:111790/1‑129 (MQ=255)
                                                                                                                      aaaGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCa                                           <  4:111790/129‑1 (MQ=255)
                                                                                                                      aaaGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGt                                   <  5:86536/137‑1 (MQ=255)
                                                                                                                      aaaGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGt                                   >  6:86536/1‑137 (MQ=255)
                                                                                                                        aGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACGGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGCCAGCTGTTTTAt                     <  5:143385/149‑1 (MQ=255)
                                                                                                                        aGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGCCAGCTGTTTTAt                     <  8:308557/149‑1 (MQ=255)
                                                                                                                             ctgccagcagGCTTCAACCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGCCCGCTGTTTTATGcc                  >  4:10774/5‑148 (MQ=255)
                                                                                                                                ccagcagGCTTCAACCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTTCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGCCAGCTGTTTTATGCCGAAg              >  3:223001/2‑149 (MQ=255)
                                                                                                                                 cagcagGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGCCAGCTGTTTTATGCCGAAGt             <  8:98738/149‑1 (MQ=255)
                                                                                                                                   gcagGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGTCAGCTGttt                        >  6:274344/1‑136 (MQ=255)
                                                                                                                                   gcagGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGTCAGCTGttt                        <  5:274344/136‑1 (MQ=255)
                                                                                                                                      ggTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAACCCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGCCAGCTGTTTTa                      >  6:331121/1‑135 (MQ=255)
                                                                                                                                      ggTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGCCAGCTGTTTTa                      <  5:331121/135‑1 (MQ=255)
                                                                                                                                           aaGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGCCAGCTGTTTTATGCCGAAGTGCGCAagag    <  1:33455/148‑1 (MQ=255)
                                                                                                                                           aaGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGCCAGCTGTTTTATGCCGAAGTGCGCAAGAgg   >  5:59074/1‑149 (MQ=255)
                                                                                                                                            aGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGCCAGCTGTTTTATGCCGAAGTGCGCAAGAggg  <  1:128349/149‑1 (MQ=255)
                                                                                                                                             gCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCGCCAGGTTTCAGAACCAGTGCCAGCTGTTTTATGCCGAAGTGCGCAAGAggg  >  3:227376/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GAAGCGCGTCGTGACAGGCTTCATCCCGTTATCGTCGCTCTGGCGATTCACACCGACGGACCAGGCTCGGTCCAAAGAGGGTGCTGCAATGGCGACTATCGACACCGAGAACCTTGCAAAGGGGCTGCAGCAGGTTGAAGCCTGCTGGCAGGCGTATACCTCGTCACAGCGGCTTGTCGAAGGCCTCGCGGGCAGCGACCGGAGCAGTGCCATTGAAACCGCCAACCGCGATTTCACCAGGTTTCAGAACCAGTGCCAGCTGTTTTATGCCGAAGTGCGCAAGAGGG  >  NC_002947/3096987‑3097273

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: