Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,050,473 T→C 8.1% P196P (CCA→CCG PP_5687 ← MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,050,4730TC 8.1% 134.3 / 2.0 49P196P (CCA→CCGPP_5687MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (22/23);  new base C (2/2);  total (24/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 = 9.44e-01

GCAGCCTCAAGCGTTTTTCTGCGATCCCCTTCCCATCAAGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGACATAATCGCCTTTGGCTGGGGTCAGTAAGTCTATGATGAGCCCGGCA  >  NC_002947/5050327‑5050617
                                                                                                                                                  |                                                                                                                                                
gcagcCTCAAGCGTTTTTCTGCGATCCCCTTCCCATCAAGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACCgg                                                                                                                                                >  3:201266/1‑149 (MQ=255)
    ccTCAAGCGTTTTTCTGCGATCCCCTTCCCATCAAGCCGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCCACACCTGAAATTTGCAGATTTCAAAAGATCCCAAGCACGCGCCTTAATCGCAACTGGGacc                                                                                                                                            >  4:264080/1‑148 (MQ=255)
      tCAAGCGTTTTTCTGCGATCCCCTTCCCATCAAGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACt                                                                                                                                          >  2:58904/1‑149 (MQ=255)
        aaGCGTTTTTCTGCGATCCCCTTCCCATCAAGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTg                                                                                                                                        <  6:120065/149‑1 (MQ=255)
          gCGTTTTTCTGCGATCCCCTTCCCATCAAGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGc                                                                                                                                       >  8:45585/1‑148 (MQ=255)
               tttCTGCGATCCCCTTCCCATCAAGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTgc                                                                                                                                  <  4:241323/148‑1 (MQ=255)
                   tGCGATCCCCTTCCCATCAAGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATc                                                                                                                             >  7:276563/1‑149 (MQ=255)
                          cccTTCCCATCAAGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGAtt                                                                                                                       <  1:58904/148‑1 (MQ=255)
                            cTTCCCATCAAGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTTAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATc                                                                                                                    <  2:254727/149‑1 (MQ=255)
                            cTTCCCATCAAGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACCGGGACACTTGCACTGCGCATCGCGATTATc                                                                                                                    >  6:286330/1‑149 (MQ=255)
                                  aTCAAGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTg                                                                                                                                        <  4:197991/123‑1 (MQ=255)
                                  aTCAAGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACCGGGACACTTGCACTGCGCATCGCGATTATCTTACgg                                                                                                              <  4:201266/149‑1 (MQ=255)
                                     aaGACGATTGCAAAATTCCTATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGTAATTTGCAGATTTCAATAGATCCCAAGCCGGCTCCTTAATCGCATCTGGGACACTTGCACTGCGCATCGCGATTATCTTAc                                                                                                                <  1:239913/144‑1 (MQ=255)
                                     aaGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGCACACTTGCACTGCGAATCGCCATTATCTAAc                                                                                                                >  2:239913/1‑144 (MQ=255)
                                       gACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGaa                                                                                                         >  6:37850/1‑149 (MQ=255)
                                                 aaaGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGAc                                                                                               >  6:5858/1‑149 (MQ=255)
                                                       cGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAAc                                                                                           <  5:5858/147‑1 (MQ=255)
                                                              gAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATaaaa                                                                                  <  6:280355/149‑1 (MQ=255)
                                                               aaaTATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTa                                                                                                                      >  3:12560/1‑112 (MQ=255)
                                                               aaaTATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTa                                                                                                                      <  4:12560/112‑1 (MQ=255)
                                                                    tCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCa                                                                            <  5:37850/149‑1 (MQ=255)
                                                                       ttGAAGACCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGa                                                                          >  7:30090/1‑148 (MQ=255)
                                                                          aaGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGTACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAAt                                                                       <  2:308230/148‑1 (MQ=255)
                                                                          aaGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAAt                                                                       <  2:199268/148‑1 (MQ=255)
                                                                          aaGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATg                                                                      <  8:276563/149‑1 (MQ=255)
                                                                             gCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACCGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCaa                                                                   <  5:286330/149‑1 (MQ=255)
                                                                               ccTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGaaaa                                                                  >  7:261535/1‑148 (MQ=255)
                                                                                 tCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAAc                                                                >  1:93528/1‑148 (MQ=255)
                                                                                       ccTCAAGCCCTGAAATTTGCAGATTTCAAAAGATCCCAATCAGGCGCCTTAATCGCATCTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAAGCAAGACGATCAGGTTATAAGTGCCAGAAATGCAAAACAATACt                                                          <  3:264080/148‑1 (MQ=255)
                                                                                       ccTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTg                                                                                                           <  4:160855/99‑1 (MQ=255)
                                                                                       ccTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTg                                                                                                           >  3:160855/1‑99 (MQ=255)
                                                                                       ccTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATaaaa                                                                                  <  4:30619/124‑1 (MQ=255)
                                                                                       ccTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATaaaa                                                                                  >  3:30619/1‑124 (MQ=255)
                                                                                         tCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATa                                                            <  2:283336/144‑1 (MQ=255)
                                                                                         tCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATAcaga                                                        >  1:283336/1‑145 (MQ=255)
                                                                                           aaCCCCTGAAATTTGCTGATTTCATAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGaaa                                                     <  1:237335/149‑1 (MQ=255)
                                                                                               ccTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGa                                                 >  6:119431/1‑149 (MQ=255)
                                                                                                             aTTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGACATAATCGCCTTTg                                   >  7:198558/1‑149 (MQ=255)
                                                                                                              tttCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGACATAATCGCCTTTgg                                  >  7:42182/1‑149 (MQ=255)
                                                                                                                tCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATTCAAAACAATACTGGAAAGCGACATAATCGCCTTTGGc                                 <  8:30090/148‑1 (MQ=255)
                                                                                                                  aaaaGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGACATAATCGCCTTTGGCTg                               <  2:60030/148‑1 (MQ=255)
                                                                                                                   aaaGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGACATAATCGCCTTTGGCTggg                             <  7:177584/149‑1 (MQ=255)
                                                                                                                         cccAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGACATAATCGCCTTTGGCTGGGGTCAg                        >  8:203194/1‑148 (MQ=255)
                                                                                                                          ccAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGACATAATCGCCTTTGGCTgg                              <  7:235590/141‑1 (MQ=255)
                                                                                                                          ccAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGACATAATCGCCTTTGGCTgg                              >  8:235590/1‑141 (MQ=255)
                                                                                                                           cAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGACATAATCGCCTTTGGCTGGGGTCAGt                       >  1:24892/1‑147 (MQ=255)
                                                                                                                                       ttAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGACATAATCGCCTTTGGCTGGGGTCAGTAAGTCTATGATGAg         >  7:41506/1‑149 (MQ=255)
                                                                                                                                              cAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGACATAATCGCCTTTGGCTGGGGTCAGTAAGTCTATGATGAGCCCgg    >  1:261483/1‑147 (MQ=255)
                                                                                                                                              cAACTCGTACACTTGAACTGCGCATCGCGCTTATCTTACTGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGACATAATCGCCTTTGGCTGGGGTCAGTAAGTCTATGATGAGCCCGGCa  <  8:261535/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCAGCCTCAAGCGTTTTTCTGCGATCCCCTTCCCATCAAGACGATTGCAAAAGTCCGATTCGGAAATATCGTTGAAGGCCCTCATTGCCTCAACCCCTGAAATTTGCAGATTTCAAAAGATCCCAAGCAGGCGCCTTAATCGCAACTGGGACACTTGCACTGCGCATCGCGATTATCTTACGGGTGAATAACCAAGACGAACAGGATAAAAGTGCCAGAAATGCAAAACAATACTGGAAAGCGACATAATCGCCTTTGGCTGGGGTCAGTAAGTCTATGATGAGCCCGGCA  >  NC_002947/5050327‑5050617

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: