Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,786,403 G→A 13.8% S89L (TCG→TTG)  alx ← alkaline‑induced transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,786,4030GA13.8% 67.6 / 3.1 29S89L (TCG→TTG) alxalkaline‑induced transporter
Reads supporting (aligned to +/- strand):  ref base G (11/14);  new base A (2/2);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

AGCGCCGCGCCCACGCCGATCATGATGGCCCGCAGGAACACCACGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAACTGGTACCAGACCCCCGCGCCGAACAGTACACCGACACTGAAATAGCCTGAATACAGCAGCAGGCTCTCGCGCA  >  NC_002947/2786259‑2786529
                                                                                                                                                |                                                                                                                              
aGCGCCGCGCCCACGCCGATCATGATGGCCCGCAGGAACACCACGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTg                                                                                                                            >  7:190965/1‑149 (MQ=255)
 gcgccgcgccCACGCCGATCATGATGGCCCGCAGGAACACCACGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGCAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGc                                                                                                                           <  5:95027/149‑1 (MQ=255)
  cgccgcgccCACGCCGATCATGATGGCCCGCAGGAACACCACGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCt                                                                                                                          >  4:168471/1‑149 (MQ=255)
   gccgcgccGACGCCGATCATGATGGCCCGCAGGAACACCACGGCCACGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTc                                                                                                                         <  8:253320/149‑1 (MQ=255)
   gccgcgccCACGCCGATCATGATGGCCCGCAGGAACACCACGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTc                                                                                                                         <  5:294845/149‑1 (MQ=255)
                cGATCATGATGGCCCGCAGGAACACCACGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCg                                                                                                            <  3:168471/149‑1 (MQ=255)
                        atgGCCCGCAGGAACACCACGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGaa                                                                                                    <  8:190965/149‑1 (MQ=255)
                          ggCCCGCAGGAACACCACGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACt                                                                                                  >  1:163655/1‑149 (MQ=255)
                              cGCAGGAACACCACGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACt                                                                                                  <  3:239432/145‑1 (MQ=255)
                                cAGGAACACCACGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTg                                                                                            >  6:345377/1‑149 (MQ=255)
                                        ccacGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTgc                                                                                    <  2:183368/149‑1 (MQ=255)
                                         cacGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCAATTGCTCGACCAAAAAACCGgtgt                                                                                                        >  3:27111/1‑128 (MQ=255)
                                         cacGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCAATTGCTCGACCAAAAAACCGgtgt                                                                                                        <  4:27111/128‑1 (MQ=255)
                                             ccaaGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATTGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACAGCGAGcc                                                                               >  1:115093/4‑149 (MQ=255)
                                                 ggATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAAc                                                                           <  2:163655/149‑1 (MQ=255)
                                                 ggATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAAc                                                                           <  5:345377/149‑1 (MQ=255)
                                                      ccccAGAACAGAAGGCGGTGCTGGTTGGGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACTGCGATTGCTCGACGAAAAAACCTGTGTATAACTCCAGTGCCGACTGCGCGCCCAACTGGTa                                                                      <  2:115093/149‑1 (MQ=255)
                                                      ccccAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAACTGGTa                                                                      >  1:294716/1‑149 (MQ=255)
                                                           gaacagaacGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTgcgcgc                                                                                <  3:175232/134‑1 (MQ=255)
                                                                 aacGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAACTgg                                                                        <  6:34142/136‑1 (MQ=255)
                                                                 aacGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAACTgg                                                                        >  5:34142/1‑136 (MQ=255)
                                                                  acGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAACTGGTACCAGACCCCcgc                                                          <  2:294716/149‑1 (MQ=255)
                                                                           tGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAACTGGTACCAGACCCCCGCGCCGAACAg                                                 >  5:212348/1‑149 (MQ=255)
                                                                                       gggATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAACTGGTACCAGACCCCCGCGCCGAACAGTACACCGACACt                                     >  1:16519/1‑149 (MQ=255)
                                                                                                     gCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAACTGGTACCAGACCCCCGCGCCGAACAGTACACCGACACTGAAATAGCCTGaa                        <  6:212348/148‑1 (MQ=255)
                                                                                                       tGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCAATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAACTGGTACCAGACCCCCGCGCCGAACAGTACACCGACACTGa                                   >  1:61692/1‑135 (MQ=255)
                                                                                                       tGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCAATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAACTGGTACCAGACCCCCGCGCCGAACAGTACACCGACACTGa                                   <  2:61692/135‑1 (MQ=255)
                                                                                                         aaGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAACTGGTACCAGACCCCCGCGCCGAACAGTACACCGACACTGAAATAGCCTGAATAcag                   >  1:317864/1‑149 (MQ=255)
                                                                                                                          gAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAACTGGTACCAGACCCCCGCGCCGAACAGTACACCGACACTGAAATAGCCTGAATACAGCAGCAGGCTCTCGCGCa  >  3:325943/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                              
AGCGCCGCGCCCACGCCGATCATGATGGCCCGCAGGAACACCACGCCCAGGATGCCCCAGAACAGAACGCGGTGCTGGTAGCGGCGGGGGATGGCGAAGTAGCTGAAGATCATCGCCATGACGAACACGTTGTCCATGGACAGCGATTGCTCGACCAAAAAACCGGTGTAGAACTCCAGTGCCGACTGCGCGCCCAACTGGTACCAGACCCCCGCGCCGAACAGTACACCGACACTGAAATAGCCTGAATACAGCAGCAGGCTCTCGCGCA  >  NC_002947/2786259‑2786529

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