Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,104,759 T→C 8.1% H567R (CAC→CGC)  PP_5355 ← sodium/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,104,7590TC 8.1% 134.1 / 2.8 49H567R (CAC→CGC) PP_5355sodium/proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (25/20);  new base C (2/2);  total (27/22)
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

AAATCGCCAGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGCGGGA  >  NC_002947/6104612‑6104903
                                                                                                                                                   |                                                                                                                                                
aaaTCGCCAGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCgtgt                                                                                                                                                  <  4:349839/148‑1 (MQ=255)
 aaTCGCCAGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCATCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGc                                                                                                                                                >  7:241071/1‑149 (MQ=255)
         ggCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCTAGATACTTTGCCCGGGGCAGA‑CCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATg                                                                                                                                       <  6:294637/149‑1 (MQ=255)
           cTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGa                                                                                                                                      <  4:339373/149‑1 (MQ=255)
           cTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGa                                                                                                                                      <  2:313462/149‑1 (MQ=255)
                               cGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGc                                                                                                                   >  6:186077/1‑148 (MQ=255)
                                       ccccATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGa                                                                                                                                      >  8:130202/1‑121 (MQ=255)
                                       ccccATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGa                                                                                                                                      <  7:130202/121‑1 (MQ=255)
                                        cccATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCcaccaccac                                                                                                          <  1:287007/148‑1 (MQ=255)
                                         ccATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCaacccccacag                                                                                                        >  4:30667/1‑147 (MQ=255)
                                                ctcGTGCGGGCGAGATACTTTGCCCGGGGCAGACCCGCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTAc                                                                                                 <  4:326246/149‑1 (MQ=255)
                                                   gggcgggcGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTAcc                                                                                                <  3:262825/147‑1 (MQ=255)
                                                        ggcgAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCaa                                                                                         <  5:186077/149‑1 (MQ=255)
                                                                  tttGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACcagcag                                                                               >  1:269201/1‑149 (MQ=255)
                                                                       ccGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAAtt                                                                          >  2:321873/1‑149 (MQ=255)
                                                                       ccGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAAtt                                                                          >  5:242764/1‑149 (MQ=255)
                                                                        cGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTc                                                                         >  1:224569/1‑149 (MQ=255)
                                                                         ggggCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCg                                                                        >  4:150432/1‑149 (MQ=255)
                                                                          gggCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGc                                                                       <  2:224569/149‑1 (MQ=255)
                                                                          gggCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGc                                                                       >  6:180696/1‑149 (MQ=255)
                                                                          gggCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCATCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGc                                                                       >  1:69023/1‑149 (MQ=255)
                                                                                aCCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCtggtgg                                                                 >  1:132721/1‑149 (MQ=255)
                                                                                aCCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCtggtgg                                                                 <  3:150432/149‑1 (MQ=255)
                                                                                        aGTGCGAGTTCTCGCGACTGTTGTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCGACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATg                                                         >  4:240012/1‑149 (MQ=255)
                                                                                             gAGGTCTCGCGACTGTTCTCCAGCGGCTCCAGCAACGCGATCTGCATGCGCGGGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGc                                                    <  8:241071/149‑1 (MQ=255)
                                                                                                 tctcGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGa                                                 >  6:109052/1‑148 (MQ=255)
                                                                                                     gcgACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATg                                            >  1:186050/1‑149 (MQ=255)
                                                                                                                 cAGCGTCTCCAGCAACGCGATCTGCATGCTCGTGTGCAAGCGGATGAGCCAGCGCCAAAACAACGCCACCAACACCGCAGCTACCACGGCAATAACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCAGAcagaag                                  >  7:134108/1‑147 (MQ=255)
                                                                                                                 cAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGa                                        <  1:338782/141‑1 (MQ=255)
                                                                                                                 cAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGa                                        >  2:338782/1‑141 (MQ=255)
                                                                                                                   gCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCcac                                                                                                                           >  8:252535/1‑56 (MQ=255)
                                                                                                                   gCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCcac                                                                                                                           <  7:252535/56‑1 (MQ=255)
                                                                                                                    cGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGTTGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAgg                              >  4:172443/1‑148 (MQ=255)
                                                                                                                      tctcCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGaaa                           <  2:269201/149‑1 (MQ=255)
                                                                                                                                     ggctcATCCGCGTGCACACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACa                 <  2:257270/140‑1 (MQ=255)
                                                                                                                                     ggctcATCCGCGTGCACACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACa                 >  1:257270/5‑144 (MQ=255)
                                                                                                                                        tcATCCGCGTGCACACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGAcagcagc                                 >  2:232969/2‑125 (MQ=255)
                                                                                                                                        tcATCCGCGTGCACACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGAcagcagc                                 <  1:232969/124‑1 (MQ=255)
                                                                                                                                        gCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGacaa          <  3:272848/148‑1 (MQ=255)
                                                                                                                                        gCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGacaa          >  4:272848/1‑148 (MQ=255)
                                                                                                                                         cATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGc                                                    >  5:335487/1‑105 (MQ=255)
                                                                                                                                         cATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGc                                                    <  6:335487/105‑1 (MQ=255)
                                                                                                                                         cATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCg              <  6:242764/143‑1 (MQ=255)
                                                                                                                                         cATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAAGCACCAACACCGacaa          >  3:288021/1‑147 (MQ=255)
                                                                                                                                          aTGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACa                 <  3:117833/139‑1 (MQ=255)
                                                                                                                                          aTGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACa                 >  4:117833/1‑139 (MQ=255)
                                                                                                                                            gcgcGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGCg     >  3:34186/1‑149 (MQ=255)
                                                                                                                                             cgcgTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGCg     <  4:34186/148‑1 (MQ=255)
                                                                                                                                              gcgTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGCgg    >  3:167605/1‑148 (MQ=255)
                                                                                                                                               cgTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGCGGGa  <  1:321873/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AAATCGCCAGGCTTTTTGCGCCTGGCGACATCGCTATAGCCCCATCTCCTCGGGCGGGCGAGATACTTTGCCCGGGGCAGACCCCCTCAGTGCGAGTTCTCGCGACTGTTCTCCAGCGTCTCCAGCAACGCGATCTGCATGCGCGTGTGCACGCGGATGAGCCAGCGCCACAACAACGCCACCACCACCGCAGCTACCACGGCAATCACCAGCAGCAATTCGCTGGTGGGCAGAATGCTGGCGGACAATGCCGACAGCAGCAGGAAAATCACCAACAGCGACAACAGCGGGA  >  NC_002947/6104612‑6104903

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