Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 813,451 G→A 9.3% A312V (GCC→GTC)  PP_0700 ← transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947813,4510GA 9.3% 116.3 / 3.2 43A312V (GCC→GTC) PP_0700transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (19/20);  new base A (2/2);  total (21/22)
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.27e-01

GATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGGCGCTCCAGCACCTGCCTCAATGGCAT  >  NC_002947/813304‑813596
                                                                                                                                                   |                                                                                                                                                 
gATGGCGCGCATTATTCTCTTCAACCCCCCCATTGACCGTTTTCGTTTCAAAAAAAATTTGCACGCTTTCAACAGCTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCt                                                                                                                                                 <  1:137978/150‑1 (MQ=255)
      gcgcATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGcc                                                                                                                                           >  2:99642/1‑150 (MQ=255)
        gcATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCCGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAg                                                                                                                                         <  2:328417/150‑1 (MQ=255)
        gcATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAg                                                                                                                                         >  1:541768/1‑150 (MQ=255)
           ttattCTCTTGAACGCCCCCTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGa                                                                                                                                      <  1:77055/150‑1 (MQ=255)
              ttCTCTTGAACGCCCCCTCTGTCAACACTCATTTCAAGAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGACTGTAGCCAGGGATTg                                                                                                                                   >  1:390569/1‑150 (MQ=255)
              ttCTCTTGAACGCCCCCTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTg                                                                                                                                   >  2:37857/1‑150 (MQ=255)
              ttCTCTTGAACGCCCCCTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGACTGCAGCCAGGGATTg                                                                                                                                   >  1:19056/1‑150 (MQ=255)
                ctctTGAACGCCCCCTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCa                                                                                                                                 <  1:99642/150‑1 (MQ=255)
                   ttGAACGCCCCCTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGc                                                                                                                              >  2:352194/1‑150 (MQ=255)
                       aCGCCCCCTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTg                                                                                                                          <  2:541768/150‑1 (MQ=255)
                        cGCCCCCTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTgg                                                                                                                         >  1:208522/1‑150 (MQ=255)
                         gCCCCCTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGCGTTGCAGGCGTTGGc                                                                                                                        >  2:171836/1‑150 (MQ=255)
                            cccTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAg                                                                                                                     <  2:208522/150‑1 (MQ=255)
                             ccTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAgg                                                                                                                                        >  2:489299/1‑130 (MQ=255)
                             ccTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAgg                                                                                                                                        <  1:489299/130‑1 (MQ=255)
                                        acTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGc                                                                                                                   <  2:501117/140‑1 (MQ=255)
                                        acTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGc                                                                                                                   >  1:501117/1‑140 (MQ=255)
                                            aTTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGaa                                                                                                     >  1:185522/1‑150 (MQ=255)
                                             tttCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCa                                                                                                                                              >  2:130100/1‑108 (MQ=255)
                                             tttCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCa                                                                                                                                              <  1:130100/108‑1 (MQ=255)
                                             tttCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCGAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAg                                                                                                    >  1:32834/1‑150 (MQ=255)
                                               tCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTc                                                                                                  <  1:463466/150‑1 (MQ=255)
                                                caagaaaaaTTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGACTGTAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTcg                                                                                                 <  2:390569/146‑1 (MQ=255)
                                                cAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGACTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTcg                                                                                                 <  2:19056/150‑1 (MQ=255)
                                                    aaaaaTTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTcgccgc                                                                                             >  1:554888/1‑150 (MQ=255)
                                                            gCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACgg                                                                                     <  1:37857/150‑1 (MQ=255)
                                                                  tttCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGa                                                                               <  2:27293/150‑1 (MQ=255)
                                                                                cGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTtcatc                                                                 <  2:554888/150‑1 (MQ=255)
                                                                                    cAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTtcatccatc                                                             >  2:511216/1‑150 (MQ=255)
                                                                                        cccGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGa                                                                                  <  2:249848/125‑1 (MQ=255)
                                                                                        cccGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGa                                                                                  >  1:249848/1‑125 (MQ=255)
                                                                                                 aGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAg                                                                                                                                         <  2:255143/61‑1 (MQ=255)
                                                                                                 aGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAg                                                                                                                                         >  1:255143/1‑61 (MQ=255)
                                                                                                   cgcACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGc                                              <  1:295216/150‑1 (MQ=255)
                                                                                                     caccacAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGccc                                            <  1:492954/150‑1 (MQ=255)
                                                                                                        cacaAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTg                                         >  2:17197/1‑150 (MQ=255)
                                                                                                        cacaAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTg                                         >  2:423092/1‑150 (MQ=255)
                                                                                                         acaAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTgg                                        <  1:511216/150‑1 (MQ=255)
                                                                                                                       gAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGgcg                          >  2:157066/1‑150 (MQ=255)
                                                                                                                                          aactgcTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGGCGCTCCAGCACCTGCCTCAAt       <  1:17197/150‑1 (MQ=255)
                                                                                                                                             tgctgGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGGCGCTCCAGCACCTGCCTCAATGGc    >  2:223873/1‑150 (MQ=255)
                                                                                                                                               ctgGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCTCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGGCGCTCCAGCACCTGCCTCAATGGCAt  <  2:75836/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTCAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGGCGCTCCAGCACCTGCCTCAATGGCAT  >  NC_002947/813304‑813596

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