Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,507,798 T→C 12.3% V35A (GTG→GCG)  PP_3105 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,507,7980TC12.3% 138.1 / 10.8 57V35A (GTG→GCG) PP_3105membrane protein
Reads supporting (aligned to +/- strand):  ref base T (23/27);  new base C (4/3);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-01

GGGGCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCATCG  >  NC_002947/3507649‑3507928
                                                                                                                                                     |                                                                                                                                  
ggggCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGt                                                                                                                                    >  2:1742051/1‑150 (MQ=255)
ggggCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGt                                                                                                                                    >  2:385241/1‑150 (MQ=255)
    cTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCg                                                                                                                                <  2:1560378/150‑1 (MQ=255)
              cccTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGGGCCGCCATTCCTGGCGGTGCCGGCGATCCTGa                                                                                                                      >  1:889687/1‑150 (MQ=255)
                cTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATc                                                                                                                    <  2:1847082/150‑1 (MQ=255)
                  gggTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCAt                                                                                                                  >  1:39009/1‑150 (MQ=255)
                   ggTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCAtt                                                                                                                 >  1:353506/1‑150 (MQ=255)
                         ttCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGAGGTGCCGGCGATCCTGATCATTATCGtt                                                                                                           <  2:475540/150‑1 (MQ=255)
                          tCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTa                                                                                                          <  1:1012766/150‑1 (MQ=255)
                          tCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTa                                                                                                          <  2:353506/150‑1 (MQ=255)
                           cTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAACTCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCTGGCGATCCTGATCATTATCGTTAt                                                                                                         <  1:1786099/150‑1 (MQ=255)
                           cTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTAt                                                                                                         <  1:15583/150‑1 (MQ=255)
                               aCGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACg                                                                                                     >  1:76123/1‑150 (MQ=255)
                                                gTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTAt                                                                                                               <  1:1599375/123‑1 (MQ=255)
                                                gTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTAt                                                                                                               >  2:1599375/1‑123 (MQ=255)
                                                    aCCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATtgtg                                                                                  >  2:1511964/1‑148 (MQ=255)
                                                    aCCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATtgtg                                                                                  <  1:1511964/148‑1 (MQ=255)
                                                       acaGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCa                                                                                           >  2:1363072/1‑136 (MQ=255)
                                                       acaGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCa                                                                                           <  1:1363072/136‑1 (MQ=255)
                                                         aGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCg                                                                           <  1:1742051/150‑1 (MQ=255)
                                                             cGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTGGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGt                                                                       <  1:1476969/150‑1 (MQ=255)
                                                              gTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAACTGTGCTTACCGCAGTg                                                                      >  1:1931881/1‑150 (MQ=255)
                                                                    ttAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACt                                                                <  2:76123/150‑1 (MQ=255)
                                                                             tCTGGTGCGGGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGACCACTTTGTAAGtt                                                       <  2:889687/150‑1 (MQ=255)
                                                                             tCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGtt                                                       >  2:615280/1‑150 (MQ=255)
                                                                                ggTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAAtt                                                                                     <  2:1047845/117‑1 (MQ=255)
                                                                                ggTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAAtt                                                                                     >  1:1047845/1‑117 (MQ=255)
                                                                                     gtgtGGCTGCTGGTCGTGATATTAGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGaga                                               >  1:1868083/1‑148 (MQ=255)
                                                                                       gtgGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAAttt                                             <  2:1909958/150‑1 (MQ=255)
                                                                                         ggCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGt                                                            >  1:331111/1‑133 (MQ=255)
                                                                                         ggCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGt                                                            <  2:331111/133‑1 (MQ=255)
                                                                                          gctgctGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCACCGCCATGAATGGCGGCGcca                                                                                                                                >  1:421458/1‑63 (MQ=25)
                                                                                          gctgctGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCACCGCCATGAATGGCGGCGcca                                                                                                                                <  2:421458/64‑2 (MQ=21)
                                                                                            tgctgGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCa                                        >  1:1612525/1‑150 (MQ=255)
                                                                                                   cGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGttctt                                 <  1:615280/150‑1 (MQ=255)
                                                                                                       atatTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAat                             <  1:385241/150‑1 (MQ=255)
                                                                                                         atTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAatat                           >  2:986573/1‑150 (MQ=255)
                                                                                                         atTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAACTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAatat                           <  2:1931881/150‑1 (MQ=255)
                                                                                                            ggCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGAGCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTc                        >  2:1061937/1‑150 (MQ=255)
                                                                                                              cccTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTAt                                                                                                         >  1:1044722/1‑67 (MQ=255)
                                                                                                              cccTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTAt                                                                                                         <  2:1044722/67‑1 (MQ=255)
                                                                                                                 tgatgatgaCCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTg                   >  2:752418/1‑150 (MQ=255)
                                                                                                                 tgatgatgaCCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTg                   >  1:1028903/1‑150 (MQ=255)
                                                                                                                  gatgatgaCCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGc                  <  2:39009/150‑1 (MQ=255)
                                                                                                                   atgatgaCCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATc                                                                                                              <  2:322687/57‑1 (MQ=255)
                                                                                                                   atgatgaCCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATc                                                                                                              >  1:322687/1‑57 (MQ=255)
                                                                                                                           cGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGt                         <  1:920664/134‑1 (MQ=255)
                                                                                                                           cGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGt                         >  2:920664/1‑134 (MQ=255)
                                                                                                                               gTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAg                                                                        <  1:727375/83‑1 (MQ=255)
                                                                                                                               gTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAg                                                                        >  2:727375/1‑83 (MQ=255)
                                                                                                                                cggcaccgccATGAATGGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGtt         <  1:636075/140‑1 (MQ=255)
                                                                                                                                cggcaccgccATGAATGGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGtt         >  2:636075/6‑145 (MQ=255)
                                                                                                                                  gcaccgccATGAATGGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTaa                                                          <  2:1435198/91‑1 (MQ=255)
                                                                                                                                  gcaccgccATGAATGGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTaa                                                          >  1:1435198/4‑94 (MQ=255)
                                                                                                                                  gcaccgccATGAATGGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCATCg  >  1:1354713/4‑150 (MQ=255)
                                                                                                                                    gccgccATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGt                                                            >  2:1059142/1‑90 (MQ=255)
                                                                                                                                    gccgccATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGt                                                            <  1:1059142/90‑1 (MQ=255)
                                                                                                                                       gccATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCgg                      <  1:911363/125‑1 (MQ=255)
                                                                                                                                       gccATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCgg                      >  2:911363/1‑125 (MQ=255)
                                                                                                                                                     |                                                                                                                                  
GGGGCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCATCG  >  NC_002947/3507649‑3507928

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: