Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,763,092 A→T 14.0% K112M (AAG→ATG)  PP_3326 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,763,0920AT14.0% 126.3 / 12.3 50K112M (AAG→ATG) PP_3326hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base T (3/4);  total (25/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 = 9.43e-01

ATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAGGCAAGACGCTG  >  NC_002947/3762944‑3763234
                                                                                                                                                    |                                                                                                                                              
aTGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTTGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCaa                                                                                                                                                <  7:265830/149‑1 (MQ=255)
 tGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCATTCAGGATCGCACCCATGTCAAg                                                                                                                                               <  2:112324/149‑1 (MQ=255)
   aTCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGa                                                                                                                                             <  7:54043/149‑1 (MQ=255)
    tCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGAt                                                                                                                                            >  5:84649/1‑149 (MQ=255)
        gggcACGGTCGAGCGCCTGGCGGAGCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGGCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCATCCTTGACATGGGTGCGa                                                                                                                                        <  5:109939/146‑1 (MQ=255)
          tCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCATCCTTGACATGGGTGCGATc                                                                                                                                      <  4:250323/149‑1 (MQ=255)
                         tGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGc                                                                                                                       >  3:22237/1‑149 (MQ=255)
                            cGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTg                                                                                                                    <  5:274425/149‑1 (MQ=255)
                               aTCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGaaa                                                                                                                 <  4:24260/149‑1 (MQ=255)
                                    gCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAg                                                                                                            >  1:26892/1‑149 (MQ=255)
                                      cGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCATCCTTGACATGGGTGCGa                                                                                                                                        >  4:1547/1‑119 (MQ=255)
                                      cGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCATCCTTGACATGGGTGCGa                                                                                                                                        <  3:1547/119‑1 (MQ=255)
                                           tgcAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGc                                                                                                      >  5:250708/1‑148 (MQ=255)
                                           tgcAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGc                                                                                                      >  5:273834/1‑148 (MQ=255)
                                                  aCTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTAcc                                                                                              <  4:22237/149‑1 (MQ=255)
                                                         ccGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCgg                                                                                       >  2:94008/1‑149 (MQ=255)
                                                                  gCACTTGCCGTGCTGCTGGACAATGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCt                                                                              >  2:181784/1‑149 (MQ=255)
                                                                    aCTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGcc                                                                                                          <  2:211737/119‑1 (MQ=255)
                                                                    aCTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGcc                                                                                                          >  1:211737/1‑119 (MQ=255)
                                                                        gCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGa                                                                         >  1:198982/1‑148 (MQ=255)
                                                                        gCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGa                                                                         >  3:98166/1‑148 (MQ=255)
                                                                          cGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAAc                                                                      <  8:84220/149‑1 (MQ=255)
                                                                          cGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAAc                                                                      <  2:26892/149‑1 (MQ=255)
                                                                            tgctgcTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTg                                                                    >  2:177977/1‑149 (MQ=255)
                                                                                gctgGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAAc                                                                      <  8:212991/143‑1 (MQ=255)
                                                                                gctgGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAAc                                                                      >  7:212991/1‑143 (MQ=255)
                                                                                 ctgGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAAAGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGAc                                                                 <  5:272017/147‑1 (MQ=255)
                                                                                    gACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACAtgct                                                            >  3:13880/1‑149 (MQ=255)
                                                                                       aaCGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCg                                                         <  1:94008/149‑1 (MQ=255)
                                                                                          gggCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGttt                                                      >  3:88907/1‑149 (MQ=255)
                                                                                              cTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGc                                                                                                           >  4:1003/1‑92 (MQ=255)
                                                                                              cTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGc                                                                                                           <  3:1003/92‑1 (MQ=255)
                                                                                                 ggTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGAt                                                <  1:177977/148‑1 (MQ=255)
                                                                                                   tCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACGTGCTGCGTTTGGCGa                                                 >  5:293004/1‑145 (MQ=255)
                                                                                                   tCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACGTGCTGCGTTTGGCGa                                                 <  6:293004/145‑1 (MQ=255)
                                                                                                       aGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACAtgctgc                                                          >  8:297999/1‑132 (MQ=255)
                                                                                                       aGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACAtgctgc                                                          <  7:297999/132‑1 (MQ=255)
                                                                                                       aGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTccc                                         <  6:250708/149‑1 (MQ=255)
                                                                                                                aCAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGAc                                <  6:84649/149‑1 (MQ=255)
                                                                                                                  aaGCAATGGGTCAATCTGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGAccc                              >  6:158405/1‑149 (MQ=255)
                                                                                                                              aaTCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCgg                  >  6:156995/1‑149 (MQ=255)
                                                                                                                                tCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACAtgct                                                            >  6:7463/1‑105 (MQ=255)
                                                                                                                                tCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCg                                                         >  5:149259/1‑108 (MQ=255)
                                                                                                                                tCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCg                                                         <  6:149259/108‑1 (MQ=255)
                                                                                                                                    gATCGCATCCTTGACATGGGTGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGa                                                 >  7:230006/1‑112 (MQ=255)
                                                                                                                                    gATCGCATCCTTGACATGGGTGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGa                                                 <  8:230006/112‑1 (MQ=255)
                                                                                                                                    gATCGCATCCTTGACATGGGTGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAgg            >  2:97240/1‑149 (MQ=255)
                                                                                                                                    gATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAgg            <  4:13880/149‑1 (MQ=255)
                                                                                                                                      tCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAGGc           >  2:13304/1‑148 (MQ=255)
                                                                                                                                          aCCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAGGCAAGAc      <  5:158405/149‑1 (MQ=255)
                                                                                                                                              aTGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAGGCAAGACGCTg  <  5:74067/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAGGCAAGACGCTG  >  NC_002947/3762944‑3763234

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