Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,724,160 C→T 9.5% L144L (CTG→CTA)  azlC ← branched‑chain amino acid transport protein AzlC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,724,1600CT 9.5% 111.9 / 2.7 42L144L (CTG→CTA) azlCbranched‑chain amino acid transport protein AzlC
Reads supporting (aligned to +/- strand):  ref base C (16/22);  new base T (2/2);  total (18/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 = 6.96e-01

GGCTACGCACACCACGGTGGCCAGGTTGCGCACCAGGGGCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAGCCCAGGCCCAGGCGCCAGCGCAGCGGCTGGCTGGA  >  NC_002947/2724022‑2724306
                                                                                                                                          |                                                                                                                                                  
ggCTACGCACACCACGGTGGCCAGGTTGCGCACCAGGGGCGCGATCAGGGCCACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAgg                                                                                                                                          >  3:261343/1‑149 (MQ=255)
   tACGCACACCACGGTGGCCAGGTTGCGCACCAGGGGCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGa                                                                                                                                       >  2:295965/1‑149 (MQ=255)
     cGCACACCACGGTGGCCAGGTTGCGCACCAGGGGCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGGCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATg                                                                                                                                     <  5:92152/149‑1 (MQ=255)
      gCACACCACGGTGGCCAGGTTGCGCACCAGGGGCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATg                                                                                                                                     >  2:107739/1‑148 (MQ=255)
       cacaCCACGGTGGCCAGGTTGCGCACCAGGGGCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTa                                                                                                                                   <  4:19421/149‑1 (MQ=255)
          accacGGGGGCCAGGTTGCGCACCAGGGGCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGaa                                                                                                                                <  2:145051/149‑1 (MQ=255)
           ccacGGTGGCCAGGTTGCGCACCAGGGGCGCGAGCAGGGCGACAACAGTGGCGACAATGGAGAAATCCAGACCGAACTGGGCCAGGGGGGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGaat                                                                                                                               <  2:176215/149‑2 (MQ=255)
            cacGGTGGCCAGGTTGCGCACCAGGGGCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAgg                                                                                                                              >  2:59975/1‑149 (MQ=255)
                      aGGTTGCGCACCAGGGGCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCAc                                                                                                                    <  1:11599/149‑1 (MQ=255)
                              cACCAGGGGCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAgg                                                                                                            <  4:224455/149‑1 (MQ=255)
                              cACCAGGGGCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAgg                                                                                                            <  1:168074/149‑1 (MQ=255)
                               aCCAGGGGCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAgg                                                                                                            >  5:202997/1‑148 (MQ=255)
                                ccAGGGGCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCa                                                                                                          >  8:26899/1‑149 (MQ=255)
                                   ggggCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTGCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATAc                                                                                                       >  7:241447/1‑149 (MQ=255)
                                              aGGGCGACAAAAGTGGCGACAATCGAGAAATCCAGACCGACCTGGTGCAGGTGTGGAATGGTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCTATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGaa                                                                                            <  8:241447/149‑1 (MQ=255)
                                                         aGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGt                                                                                 <  2:264785/149‑1 (MQ=255)
                                                               gACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTAct                                                                           >  2:137970/1‑149 (MQ=255)
                                                                   aTGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAATAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCa                                                                                                          >  7:66631/1‑114 (MQ=255)
                                                                   aTGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAATAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCa                                                                                                          <  4:307858/114‑1 (MQ=255)
                                                                   aTGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAATAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCa                                                                                                          <  8:66631/114‑1 (MQ=255)
                                                                   aTGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAATAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCa                                                                                                          >  3:307858/1‑114 (MQ=255)
                                                                        gaAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGt                                                                    >  4:107124/1‑147 (MQ=255)
                                                                                gACCGAACTGGTCCAGGTGTGGAATGTTCGTGCCCAGCAGTATGCCGGCGAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCgaag                                                          <  4:212215/149‑1 (MQ=255)
                                                                                          gTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGGCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATtcgtc                                                  <  2:108924/147‑1 (MQ=255)
                                                                                            ccAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGt                                               >  7:11553/1‑148 (MQ=255)
                                                                                             cAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTca                                             <  1:295965/149‑1 (MQ=255)
                                                                                               ggTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTcagc                                           <  6:130432/149‑1 (MQ=255)
                                                                                                      aaTGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAgc                                    >  1:152518/1‑149 (MQ=255)
                                                                                                                cccAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAgcccaggcccag                          <  6:202997/149‑1 (MQ=255)
                                                                                                                  cagcagTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAGCccaggcccaggc                        <  4:261343/149‑1 (MQ=255)
                                                                                                                           gCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAGCCCAGGCCCAGgcg                       >  6:104019/1‑141 (MQ=255)
                                                                                                                           gCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAGCCCAGGCCCAGgcg                       <  5:104019/141‑1 (MQ=255)
                                                                                                                                    ggTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAgcccagg                               >  2:139367/1‑124 (MQ=255)
                                                                                                                                    ggTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAgcccagg                               <  1:139367/124‑1 (MQ=255)
                                                                                                                                    ggTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAGCCCAGGCCCAGGCGCCAGCGCAGCggctggc      <  1:137970/149‑1 (MQ=255)
                                                                                                                                     gTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAGCCCAGGCCCAGGCGCCAGCGCAGCggctggct     <  7:26899/149‑1 (MQ=255)
                                                                                                                                     gTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAGCCCAGGCCCAGGCGCCAGCGCAGCggctggct     <  8:11553/149‑1 (MQ=255)
                                                                                                                                      tGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAGCCCAGGCCCAGGCGCCAGCGCAGCGgctggctg    >  8:227119/1‑149 (MQ=255)
                                                                                                                                        aaCAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATtc                                                     <  7:172360/98‑1 (MQ=255)
                                                                                                                                        aaCAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATtc                                                     >  8:172360/1‑98 (MQ=255)
                                                                                                                                        aaCAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAGCCCAGGCCCAGGCGCCAGCGCAGCGGCTGGCTGGa  <  7:227119/149‑1 (MQ=255)
                                                                                                                                        aaCAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAGCCCAGGCCCAGGCGCCAGCGCAGCGGCTGGCTGGa  >  8:249546/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
GGCTACGCACACCACGGTGGCCAGGTTGCGCACCAGGGGCGCGATCAGGGCGACAAAAGTGGCGACAATGGAGAAATCCAGACCGAACTGGTCCAGGTGTGGAATGTTCTGGCCCAGCAGTATGCCGGCCAGGGTGAACAGGTTCCAGGCGATGTAGAAGGTCAGGCCCACCCCCAAGGCATACCAGCGGTTGAATTGCTGCCGGTCGTACTGGCTGGTCAGGGCGAAGAATTCGTCGGTCAGCAGAAAGCCCAGGCCCAGGCGCCAGCGCAGCGGCTGGCTGGA  >  NC_002947/2724022‑2724306

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: ‑