Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,104,816 G→A 6.9% H567Y (CAC→TAC)  PP_5355 ← sodium/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,104,8160GA 6.9% 166.7 / 2.0 58H567Y (CAC→TAC) PP_5355sodium/proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (31/23);  new base A (2/2);  total (33/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.60e-01

CAGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGCGG  >  NC_002947/6104675‑6104957
                                                                                                                                             |                                                                                                                                             
cAGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCAcgc                                                                                                                                          <  2:272769/147‑1 (MQ=255)
cAGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCAcgc                                                                                                                                          >  1:272769/1‑147 (MQ=255)
 aGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGCACACGCGCAt                                                                                                                                      <  2:226703/150‑1 (MQ=255)
    cTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCgg                                                                                                                                        >  2:1419974/1‑145 (MQ=255)
    cTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCgg                                                                                                                                        <  1:1419974/145‑1 (MQ=255)
         tGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAgc                                                                                                                              >  1:1014443/1‑150 (MQ=255)
                   gACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCAcaacaa                                                                                                                    >  2:1099949/1‑150 (MQ=255)
                      aTCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATACGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGc                                                                                                                 >  1:785449/1‑150 (MQ=255)
                       tCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGcc                                                                                                                >  1:412661/1‑150 (MQ=255)
                                ccccATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCaccaccacc                                                                                                       >  1:199636/1‑150 (MQ=255)
                                  ccATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGc                                                                                                     >  1:294619/1‑150 (MQ=255)
                                   cATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCa                                                                                                    <  1:1099949/150‑1 (MQ=255)
                                         ctcGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCatcaccacc                                                                                                       >  2:1235266/1‑141 (MQ=255)
                                         ctcGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCatcaccacc                                                                                                       <  1:1235266/141‑1 (MQ=255)
                                         ctcGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTAcc                                                                                              >  1:1293605/1‑150 (MQ=255)
                                              gcgggcgAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCa                                                                                                                                >  1:1159137/1‑111 (MQ=255)
                                              gcgggcgAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCa                                                                                                                                <  2:1159137/111‑1 (MQ=255)
                                                  gcgAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATc                                                                                     >  1:483086/1‑150 (MQ=255)
                                                     agaTACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCAcc                                                                                  >  1:1305613/1‑150 (MQ=255)
                                                       aTACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACcag                                                                                <  2:785449/150‑1 (MQ=255)
                                                         aCTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACcagca                                                                              <  2:696681/150‑1 (MQ=255)
                                                           tttGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCagcagc                                                                            <  2:1207128/150‑1 (MQ=255)
                                                                ccGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTc                                                                       >  1:802998/1‑150 (MQ=255)
                                                                   gggCAGGCCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGa                                                                                                                                    <  1:310662/86‑1 (MQ=255)
                                                                   gggCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGa                                                                                                                                    >  2:310662/1‑86 (MQ=255)
                                                                    ggCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCtg                                                                   >  1:494284/1‑150 (MQ=255)
                                                                      cAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCtggt                                                                 <  2:1293605/150‑1 (MQ=255)
                                                                         aCCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTggg                                                              <  2:802998/150‑1 (MQ=255)
                                                                           ccccTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCa                                                            <  2:199636/150‑1 (MQ=255)
                                                                              cTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATCCGCGTGCACACGCGCATgca                                                                                                                                   >  2:1364717/1‑74 (MQ=37)
                                                                              cTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATCCGCGTGCACACGCGCATgca                                                                                                                                   <  1:1364717/76‑3 (MQ=37)
                                                                               tCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCAGCAGCAGCAATTCGCGGGTTGGCAGAAt                                                        >  1:314770/1‑150 (MQ=255)
                                                                                  gTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCCCCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCt                                                     <  1:746676/150‑1 (MQ=255)
                                                                                        gTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGa                                               <  2:412661/150‑1 (MQ=255)
                                                                                         ttCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGa                                                                                                                                    >  2:332747/1‑64 (MQ=255)
                                                                                         ttCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGa                                                                                                                                    <  1:332747/64‑1 (MQ=255)
                                                                                         ttCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCtggtg                                                                <  2:207653/132‑1 (MQ=255)
                                                                                         ttCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCtggtg                                                                >  1:207653/1‑132 (MQ=255)
                                                                                         ttCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCCGGCgttc                                              <  2:494284/150‑4 (MQ=255)
                                                                                              gcgACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGc                                         >  1:839517/1‑150 (MQ=255)
                                                                                                 ccTGTTCTCCAGCCGCTCCACCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGa                                      <  2:1014443/149‑1 (MQ=255)
                                                                                                    gTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGAcag                                   <  2:839517/150‑1 (MQ=255)
                                                                                                          cAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGTATGAGCCAGCGCCACAACAACGCCACCACCACCGCCACCACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGAcagcagcag                             >  1:1095482/1‑150 (MQ=255)
                                                                                                             cGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCGATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGaa                          >  1:1309620/1‑150 (MQ=255)
                                                                                                                 tcCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATc                      >  1:1185695/1‑150 (MQ=255)
                                                                                                                    aGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCGATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCAcc                   <  2:1309620/150‑1 (MQ=255)
                                                                                                                    aGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGAcagcagcag                             <  2:33783/140‑1 (MQ=255)
                                                                                                                    aGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGAcagcagcag                             >  1:33783/1‑140 (MQ=255)
                                                                                                                    aGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCAcc                   >  2:855987/1‑150 (MQ=255)
                                                                                                                    aGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGG‑AAATCACCa                  >  1:682396/1‑150 (MQ=255)
                                                                                                                       aaCGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCt                                                                    >  1:731633/1‑98 (MQ=255)
                                                                                                                       aaCGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCt                                                                    <  2:731633/98‑1 (MQ=255)
                                                                                                                       aaCGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAAc                >  2:1306424/1‑150 (MQ=255)
                                                                                                                            gATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGa           >  2:50057/1‑150 (MQ=255)
                                                                                                                                 tcatccgcgTGCACACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGacaaca      >  1:426324/6‑150 (MQ=255)
                                                                                                                                    tGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCAc                                                                                            <  2:1276192/61‑1 (MQ=255)
                                                                                                                                    tGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCAc                                                                                            >  1:1276192/1‑61 (MQ=255)
                                                                                                                                     gcgcGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGCgg  >  2:896640/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                             
CAGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGCGG  >  NC_002947/6104675‑6104957

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