Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,543,049 G→T 8.8% G263W (GGG→TGG)  PP_1353 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,543,0490GT 8.8% 133.8 / 2.5 46G263W (GGG→TGG) PP_1353membrane protein
Reads supporting (aligned to +/- strand):  ref base G (19/22);  new base T (2/2);  total (21/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 = 8.55e-01

AAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCATTTCTGG  >  NC_002947/1542906‑1543198
                                                                                                                                                |                                                                                                                                                     
aaGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCg                                                                                                                                                 >  2:477481/1‑150 (MQ=255)
   aCCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACa                                                                                                                                              <  2:1081701/150‑1 (MQ=255)
      cGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTc                                                                                                                                           <  1:611180/150‑1 (MQ=255)
               aGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGc                                                                                                                                  >  2:216313/1‑150 (MQ=255)
                    gTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCgg                                                                                                                             <  1:216313/150‑1 (MQ=255)
                            gACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGCTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTcg                                                                                                                     >  1:1254776/1‑150 (MQ=255)
                             aCGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTcgc                                                                                                                    <  1:1076943/150‑1 (MQ=255)
                                 tGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGCTTCACACTTGCGCCAGCGCAAATGGAACAAGTATGCGCGTGACCAATTACAGGCAGCACAAATCTACATGCCTTTTCCGCAGCGGGTTATTCGCgtgg                                                                                                                <  1:1391493/150‑1 (MQ=255)
                                           gtCGATCACTGTAGCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTGAACGAGTATGTGCGTGACCGGTTGAAGGGCGAAGGGATCGACATGCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGAc                                                                                                      <  2:328009/149‑1 (MQ=255)
                                           cTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAATGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGAc                                                                                                      >  2:527841/1‑150 (MQ=255)
                                                tCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCaa                                                                                                 <  1:477481/150‑1 (MQ=255)
                                                    cGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCa                                                                                             <  2:1236773/150‑1 (MQ=255)
                                                               gcTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCcgcg                                                                                  >  1:272753/1‑150 (MQ=255)
                                                                     ggACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATg                                                                                       >  2:433335/1‑139 (MQ=255)
                                                                     ggACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATg                                                                                       <  1:433335/139‑1 (MQ=255)
                                                                      gACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAATGGa                                                                                                                                                    <  1:508484/77‑1 (MQ=255)
                                                                      gACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAATGGa                                                                                                                                                    >  2:508484/1‑77 (MQ=255)
                                                                         cAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAATGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCgg                                                                        <  1:527841/150‑1 (MQ=255)
                                                                         aaaGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCgg                                                                        >  1:276225/2‑150 (MQ=255)
                                                                            gACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGgc                                                                     >  2:7230/1‑150 (MQ=255)
                                                                               cTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCa                                                                  >  1:424489/1‑150 (MQ=255)
                                                                               cTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCa                                                                  >  1:483547/1‑150 (MQ=255)
                                                                                          tGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACg                                                                                               <  2:1380354/111‑1 (MQ=255)
                                                                                          tGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGCTATTCGCGTGGTGCAGGAGACTGCAACg                                                                                               >  1:1380354/1‑111 (MQ=255)
                                                                                                 tGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAACGCGGCCCCAGCtttt                                                 >  1:23374/1‑150 (MQ=255)
                                                                                                  gATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGtttttt                                                >  1:405040/1‑150 (MQ=255)
                                                                                                         aTGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTc                                               <  1:1419392/144‑1 (MQ=255)
                                                                                                         aTGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTc                                               >  2:1419392/1‑144 (MQ=255)
                                                                                                           gCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCg                                                             <  1:1356121/128‑1 (MQ=255)
                                                                                                           gCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCg                                                             >  2:1356121/1‑128 (MQ=255)
                                                                                                            cTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGa                                      <  2:1138179/150‑1 (MQ=255)
                                                                                                               aaCGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTa                                   >  2:1069163/1‑150 (MQ=255)
                                                                                                                aCGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTac                                  <  2:483547/150‑1 (MQ=255)
                                                                                                                 cGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTaca                                 <  1:1374706/150‑1 (MQ=255)
                                                                                                                 cGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTaca                                 <  2:276225/150‑1 (MQ=255)
                                                                                                                    gTATGCGCGTAACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACtt                              <  2:975226/150‑1 (MQ=255)
                                                                                                                     tATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTg                             <  2:23374/150‑1 (MQ=255)
                                                                                                                      aTGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGc                            <  2:996668/150‑1 (MQ=255)
                                                                                                                         cgcgTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTca                         <  1:1069163/150‑1 (MQ=255)
                                                                                                                              gACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAg                                                                                            >  2:1182082/1‑78 (MQ=255)
                                                                                                                              gACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAg                                                                                            <  1:1182082/78‑1 (MQ=255)
                                                                                                                                ccGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAAt                  >  1:584800/1‑150 (MQ=255)
                                                                                                                                      aaaGGCCGAAGGGATCGACATTCCGTTTCTGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACATTTGCGCACAGCAATTTCCAg            <  2:1254776/149‑1 (MQ=255)
                                                                                                                                       aaGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGCCCCTGATTaca                                 <  2:707703/128‑1 (MQ=255)
                                                                                                                                       aaGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTaca                                 >  1:707703/1‑128 (MQ=255)
                                                                                                                                             gAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCATTTc     >  1:258124/1‑150 (MQ=255)
                                                                                                                                                gggATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCATTTCTgg  <  1:7230/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
AAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCATTTCTGG  >  NC_002947/1542906‑1543198

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