Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,895,712 T→A 10.4% L56Q (CTG→CAG)  PP_4304 → VIC family cation transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,895,7120TA10.4% 131.2 / 5.5 48L56Q (CTG→CAG) PP_4304VIC family cation transporter
Reads supporting (aligned to +/- strand):  ref base T (20/23);  new base A (2/3);  total (22/26)
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.25e-01

CGCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCG  >  NC_002947/4895567‑4895855
                                                                                                                                                 |                                                                                                                                               
cgcgAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCctgct                                                                                                                                              <  6:177020/149‑1 (MQ=255)
  cgAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTCATCCACGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTgg                                                                                                                                            >  2:58936/1‑149 (MQ=255)
   gAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGc                                                                                                                                           >  5:204776/1‑149 (MQ=255)
   gAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGc                                                                                                                                           >  5:216137/1‑149 (MQ=255)
       ggCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGt                                                                                                                                       <  8:4489/149‑1 (MQ=255)
        gcgctTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCAGCAGGCCggcg                                                                                                                                      <  1:128302/146‑3 (MQ=255)
        gCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCg                                                                                                                                         >  1:45191/1‑146 (MQ=255)
        gCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCg                                                                                                                                         <  2:45191/146‑1 (MQ=255)
          tctTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGAATCGATGAGGTGCACCAGGGCTACGCCGGCATTCTTGCCGTTATc                                                                                                                                    >  4:166342/1‑149 (MQ=255)
          tctTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCATCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATc                                                                                                                                    >  1:96198/1‑149 (MQ=255)
              cGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGa                                                                                                                                  <  3:147663/147‑1 (MQ=255)
              cGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCAGCAGGCCggcgtagc                                                                                                                                  <  7:180210/147‑7 (MQ=255)
               gTCATCGTCTTCCAGACCTACACCGTAGCCGGTCGCCGCTTCGACAAGATCATGTTGCTGATCATCCTGGCCAGCCTGGTAACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGt                                                                                                                                <  5:252279/148‑1 (MQ=255)
                  aTCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTgggg                                                                                                                            <  2:154060/149‑1 (MQ=255)
                  aTCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTggg                                                                                                                             <  2:47528/148‑1 (MQ=255)
                           cAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACcgc                                                                                                                   >  2:11649/1‑149 (MQ=255)
                           cAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACcgc                                                                                                                   <  6:119567/149‑1 (MQ=255)
                            aGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCctg                                                                                                                                                <  7:52396/119‑1 (MQ=255)
                            aGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCctg                                                                                                                                                >  8:52396/1‑119 (MQ=255)
                            aGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCAGCAGGCCggc                                                                                                                                       <  7:97177/128‑2 (MQ=255)
                            aGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCAGCAGGCCggc                                                                                                                                       >  8:97177/1‑127 (MQ=255)
                                  acacCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTc                                                                                                            >  2:32005/1‑149 (MQ=255)
                                            cGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAAt                                                                                                  <  3:166342/149‑1 (MQ=255)
                                                          aCAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAAt                                                                                                  >  4:102700/1‑135 (MQ=255)
                                                          aCAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAAt                                                                                                  <  3:102700/135‑1 (MQ=255)
                                                                 ccTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTg                                                                               <  1:188539/147‑1 (MQ=255)
                                                                 ccTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTg                                                                               >  2:188539/1‑147 (MQ=255)
                                                                     ttGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTcc                                                                                                           >  8:105666/1‑115 (MQ=255)
                                                                     ttGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTcc                                                                                                           <  7:105666/115‑1 (MQ=255)
                                                                     ttGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCg                                                                           <  3:72978/147‑1 (MQ=255)
                                                                     ttGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGc                                                                          <  6:204776/148‑1 (MQ=255)
                                                                       gCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGt                                                                                                                                >  8:119986/1‑92 (MQ=255)
                                                                       gCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGt                                                                                                                                <  7:119986/92‑1 (MQ=255)
                                                                        cTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGccc                                                                        >  8:32834/1‑147 (MQ=255)
                                                                           atcatcCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACAt                                                                                              <  7:102685/122‑1 (MQ=255)
                                                                           atcatcCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACAt                                                                                              >  8:102685/1‑122 (MQ=255)
                                                                                    gCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCt                                                                <  4:205655/143‑1 (MQ=255)
                                                                                    gCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCt                                                                >  3:205655/1‑143 (MQ=255)
                                                                                                   gTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTa                                           <  6:216137/149‑1 (MQ=255)
                                                                                                     gATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACg                                         >  8:211250/1‑149 (MQ=255)
                                                                                                          tCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCAc                                                                                           >  2:158474/1‑94 (MQ=255)
                                                                                                          tCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCAc                                                                                           <  1:158474/94‑1 (MQ=255)
                                                                                                                gCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGAt                              >  4:240694/1‑149 (MQ=255)
                                                                                                                  aTCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGAtt                             >  4:182673/1‑148 (MQ=255)
                                                                                                                     gATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTAc                          <  7:211250/148‑1 (MQ=255)
                                                                                                                     gATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTAc                          <  1:58936/148‑1 (MQ=255)
                                                                                                                       tGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTgg                       <  6:45123/149‑1 (MQ=255)
                                                                                                                       tGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCAGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGACTGGTCGCTTTACTgg                       >  6:271093/1‑149 (MQ=255)
                                                                                                                                   gggCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCgg           <  3:182673/149‑1 (MQ=255)
                                                                                                                                         cGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGTCTGGTCGATTTACTGGCGATCGTGCCGGGGATc      <  1:32005/148‑1 (MQ=255)
                                                                                                                                            cGGCCAGCAGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTcg  >  1:115050/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CGCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCG  >  NC_002947/4895567‑4895855

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