Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,113,722 A→C 8.0% Y203D (TAT→GAT)  PP_1882 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,113,7220AC 8.0% 129.2 / 4.8 50Y203D (TAT→GAT) PP_1882membrane protein
Reads supporting (aligned to +/- strand):  ref base A (24/22);  new base C (2/2);  total (26/24)
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.10e-01

GCTCACGCACAACAACCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGAGGTCGATCCTGTGCCCCGTGGCCGCGTAGA  >  NC_002947/2113576‑2113871
                                                                                                                                                  |                                                                                                                                                     
gCTCACGCACAACAACCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAgg                                                                                                                                                    <  1:723795/150‑1 (MQ=255)
  tCACGCACAACAACCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGgc                                                                                                                                                  >  1:154769/1‑150 (MQ=255)
        acaacaACCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTcgc                                                                                                                                            <  1:196716/150‑1 (MQ=255)
            caacCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCaa                                                                                                                                        <  2:674285/150‑1 (MQ=255)
             aacCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAAt                                                                                                                                       <  2:790264/150‑1 (MQ=255)
             aacCCGTACACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAAt                                                                                                                                       <  2:154769/150‑1 (MQ=255)
                      aCCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCgg                                                                                                                              >  2:546095/1‑150 (MQ=255)
                      aCCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCgg                                                                                                                              >  2:386018/1‑150 (MQ=255)
                       ccAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTTATCAGGgcgacaacg                                                                                                                                           >  2:419635/1‑130 (MQ=255)
                       ccAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTTATCAGGgcgacaacg                                                                                                                                           <  1:419635/136‑7 (MQ=255)
                            gCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTAtcg                                                                                                                        >  2:563058/1‑150 (MQ=255)
                            gCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTAtcg                                                                                                                        >  2:656155/1‑150 (MQ=255)
                                    cAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTc                                                                                                                >  2:710376/1‑150 (MQ=255)
                                     aGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAAt                                                                                                                                       >  2:838026/1‑126 (MQ=255)
                                     aGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAAt                                                                                                                                       <  1:838026/126‑1 (MQ=255)
                                        cgcCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGAtt                                                                                                            <  1:623788/150‑1 (MQ=255)
                                         gcCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGAttt                                                                                                           >  1:373935/1‑150 (MQ=255)
                                                   cAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGAt                                                                                                             >  1:143190/1‑138 (MQ=255)
                                                   cAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGAt                                                                                                             <  2:143190/138‑1 (MQ=255)
                                                      cTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGAtctc                                                                                              <  1:563058/150‑1 (MQ=255)
                                                       tCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGTTATCGTCGCGTTCGATTTCCAGCAGGATctcc                                                                                             >  1:15128/1‑150 (MQ=255)
                                                             gTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGc                                                                                       <  2:373935/150‑1 (MQ=255)
                                                                         atcgatGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGa                                                                                                                                   <  2:448662/94‑1 (MQ=255)
                                                                         atcgatGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGa                                                                                                                                   >  1:448662/1‑94 (MQ=255)
                                                                              tgtgACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCt                                                                      <  1:386018/150‑1 (MQ=255)
                                                                                  aCCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGc                                                                  >  1:234469/1‑150 (MQ=255)
                                                                                    ccAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGcat                                                                >  1:335480/1‑150 (MQ=255)
                                                                                           gtgtTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGGAGGATCTTCTCACGCTCGTCGCTGATCGAGCTGCGCATCATGTaa                                                         <  1:412501/150‑1 (MQ=255)
                                                                                                 aTCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAg                                                   >  2:9242/1‑150 (MQ=255)
                                                                                                  tCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCTCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATg                                                            >  1:619940/1‑140 (MQ=255)
                                                                                                  tCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATg                                                            <  2:619940/140‑1 (MQ=255)
                                                                                                           tCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCTCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGaa                                              <  2:174089/145‑1 (MQ=255)
                                                                                                           tCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGaa                                              >  1:174089/1‑145 (MQ=255)
                                                                                                           tCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGCGGCTGCTCACGCTCGTCGATGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGtt                                         >  1:438393/1‑150 (MQ=255)
                                                                                                               aTAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGAtctc                                                                                              <  2:702693/93‑1 (MQ=255)
                                                                                                               aTAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGAtctc                                                                                              >  1:702693/1‑93 (MQ=255)
                                                                                                                tAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATg                                                            >  1:475465/1‑126 (MQ=255)
                                                                                                                tAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATg                                                            <  2:475465/126‑1 (MQ=255)
                                                                                                                    gggCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGa                                >  2:799123/1‑150 (MQ=255)
                                                                                                                    gggCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCCATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGa                                <  1:9242/150‑1 (MQ=255)
                                                                                                                        aaCGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCTTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGa                                                                    >  1:125401/1‑110 (MQ=255)
                                                                                                                        aaCGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGa                                                                    <  2:125401/110‑1 (MQ=255)
                                                                                                                         aCGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCAGCATGTAAGACGAGTTGAACGGTTGGTACTGGAGGTCg                           >  1:830588/1‑150 (MQ=255)
                                                                                                                                    gTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGa                                                                             >  2:290185/1‑89 (MQ=255)
                                                                                                                                    gTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGa                                                                             <  1:290185/89‑1 (MQ=255)
                                                                                                                                       agcGCCCTTATCAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGa                                                    >  1:95207/3‑111 (MQ=255)
                                                                                                                                       agcGCCCTTATCAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGa                                                    <  2:95207/109‑1 (MQ=255)
                                                                                                                                         cGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGAGGTCGATCCTGTGCCCCGTgg           <  1:546095/150‑1 (MQ=255)
                                                                                                                                           cccTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGAGGTCGATCCTGTGCCCCGTGGcc         <  2:335480/150‑1 (MQ=255)
                                                                                                                                           cccTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGAGGTCGATCCTGTGCCCCGTGGcc         >  1:306211/1‑150 (MQ=255)
                                                                                                                                                  aaGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGAGGTCGATCCTGTGCCCCGTGGCCGCGTAGa  >  1:95806/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCTCACGCACAACAACCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGAGGTCGATCCTGTGCCCCGTGGCCGCGTAGA  >  NC_002947/2113576‑2113871

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