Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,056,321 C→T 7.4% G14D (GGC→GAC)  PP_2668 ← ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,056,3210CT 7.4% 157.4 / 2.3 54G14D (GGC→GAC) PP_2668ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (2/2);  total (28/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 = 8.99e-01

GCACACGCGGATGTCGCCACGCTGCAGGTCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTCGCATCGATGACCGAAACGTCACCGCTGACACCGTTGGTGG  >  NC_002947/3056173‑3056467
                                                                                                                                                    |                                                                                                                                                  
gCACACGCGGATGTCGCCACGCTGCAGGTCGTACAGGCTCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGTGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGc                                                                                                                                                    <  3:42189/149‑1 (MQ=255)
gCACACGCGGATGTCGCCACGCTGCAGGTCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGc                                                                                                                                                    >  7:273302/1‑149 (MQ=255)
       cgGATGTCGCCACGCTGCAGGTCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGc                                                                                                                                             <  8:273302/149‑1 (MQ=255)
        ggATGTCGCCACGCTGCAGGTCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCTACTTTCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCa                                                                                                                                            <  5:101760/149‑1 (MQ=255)
                ccACGCTGCAGGTCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCAc                                                                                                                                    >  1:192772/1‑149 (MQ=255)
                ccACGCTGCAGGTCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCAc                                                                                                                                    >  5:92637/1‑149 (MQ=255)
                     cTGCAGGTCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCAc                                                                                                                               <  8:98456/149‑1 (MQ=255)
                      tGCAGGTCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACt                                                                                                                              >  3:107531/1‑149 (MQ=255)
                      tGCAGGTCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTACCGGAGTGCCTCTCGCTGGCCATAGGCAAAGCTCACCTCACt                                                                                                                              >  4:13183/1‑149 (MQ=255)
                         aGGTCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGAc                                                                                                                           >  1:128874/1‑149 (MQ=255)
                           gTCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGTCATAGGCAAAGCTCACGTCACTGACGt                                                                                                                         >  8:189828/1‑149 (MQ=255)
                           gTCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGTCATAGGCAAAGCTCACGTCACTGACGt                                                                                                                         >  8:194692/1‑149 (MQ=255)
                            tCGTACAGGCGCGTGAGCAGGGTGAACAAGGTCGAATTGCCGGCACAATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGATGAAACCCACCTGCCGGAGTGCCACGCGCTGGCCATAGGCAAAGCTCACGTCACTGagaac                                                                                                                        >  7:47111/1‑145 (MQ=255)
                            tCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACt                                                                                                                              <  8:47111/143‑1 (MQ=255)
                                  acgcgcGAGAGCAGGGCGATCAAGCTCGACTTGCCGGCACCATTGGGCCCCAGCAAAGCGGTGAAGCGCCCGGCTCCCAGGCTGAAACCGACCTGCCGCAGTGCCTCGCGCTGCCCATACGCAACGCTCa                                                                                                                                     >  4:119971/3‑130 (MQ=255)
                                  aGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGTGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCa                                                                                                                                     <  3:119971/130‑1 (MQ=255)
                                   ggCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCg                                                                                                                 <  7:152646/149‑1 (MQ=255)
                                    gcgcGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCa                                                                                                                                            <  6:239364/121‑1 (MQ=255)
                                    gcgcGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCa                                                                                                                                            >  5:239364/1‑121 (MQ=255)
                                      gcgTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGtt                                                                                                               >  4:156106/1‑148 (MQ=255)
                                           aGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGgt                                                                                                         <  2:128874/149‑1 (MQ=255)
                                                         ggTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGcc                                                                                           >  4:130316/1‑149 (MQ=255)
                                                           tCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAgg                                                                                                                    >  1:312010/1‑122 (MQ=255)
                                                           tCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAgg                                                                                                                    <  2:312010/122‑1 (MQ=255)
                                                           tCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGcccc                                                                                         <  2:192772/149‑1 (MQ=255)
                                                                     ggCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAgcg                                                                               >  6:323734/1‑149 (MQ=255)
                                                                          cATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGAc                                                                                                                           <  3:5012/100‑1 (MQ=255)
                                                                          cATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGAc                                                                                                                           >  4:5012/1‑100 (MQ=255)
                                                                          cATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGa                                                                          >  1:131051/1‑149 (MQ=255)
                                                                           aTTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGAc                                                                         >  8:117002/1‑149 (MQ=255)
                                                                              ggggCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAg                                                                                       <  3:187314/132‑1 (MQ=255)
                                                                              ggggCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAg                                                                                       >  4:187314/1‑132 (MQ=255)
                                                                                 gCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGa                                                                                                                      >  1:63290/1‑98 (MQ=255)
                                                                                 gCCCACCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGa                                                                                                                      <  2:63290/98‑1 (MQ=255)
                                                                                     aGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGAc                                                               >  3:50646/1‑149 (MQ=255)
                                                                                     aGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGAc                                                               >  1:78172/1‑149 (MQ=255)
                                                                                       cAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACtt                                                             <  2:78172/149‑1 (MQ=255)
                                                                                         aTGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGa                                                           <  2:131051/149‑1 (MQ=255)
                                                                                             ggcgAAGCGCCCGGGTGCCAGGCTGAAACGGACCTACCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGGCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCAc                                                       <  3:13183/149‑1 (MQ=255)
                                                                                               cgAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACtt                                                     >  6:210624/1‑149 (MQ=255)
                                                                                                gAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACttt                                                    <  4:107531/149‑1 (MQ=255)
                                                                                                    ctcccGTGTGCCAGGCTGAATCCGACCTGCCGGTGTGCCTCGCGCTGTCCATAGGCAAAGGTAACGTCACTGACGTCGAGGGCGTTCTTGGTTTCACCACCACGCCCCAGGTGTAGCGCCCGAGCTTCTCCGACTTGAGGTCTTTGat                                                 <  4:50646/146‑2 (MQ=255)
                                                                                                          ggTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTc                                          <  5:323734/149‑1 (MQ=255)
                                                                                                               cAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTCGCATc                                     <  7:117002/149‑1 (MQ=255)
                                                                                                                aGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTCGCATc                                     >  2:245640/1‑148 (MQ=255)
                                                                                                                 ggCTGAAACCGACTTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTCGCATCGa                                   >  4:128306/1‑149 (MQ=255)
                                                                                                                 ggCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTCGCATCGa                                   <  5:210624/149‑1 (MQ=255)
                                                                                                                 ggCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTCGCATCGa                                   <  1:245640/149‑1 (MQ=255)
                                                                                                                   cTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGTCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTCGCATCGATg                                 <  7:194692/149‑1 (MQ=255)
                                                                                                                   cTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGTCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTCGCATCGATg                                 <  7:189828/149‑1 (MQ=255)
                                                                                                                    tGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAg                                                                                       >  2:289871/1‑94 (MQ=255)
                                                                                                                                 ccGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCCCGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCACAGGGGTAGCGCCCGACCTTCACCGCCTTGAGCACTTTGAGCTGCTTCGCATCGATGCCCGAAGCGTCAcc                   >  8:164518/1‑149 (MQ=255)
                                                                                                                                    gAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTCGCATCGATGACCGAAACGTCACCGCt                <  3:130316/149‑1 (MQ=255)
                                                                                                                                         ccTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTCGCATCGATGACCGAAACGTCACCGCTGacac           >  2:238570/1‑149 (MQ=255)
                                                                                                                                         ccTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTCGCATCGATGACCGAAACGTCACCGCTGacac           >  5:290785/1‑149 (MQ=255)
                                                                                                                                                   gCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTCGCATCGATGACCGAAACGTCACCGCTGACACCGTtggtgg  <  1:238570/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCACACGCGGATGTCGCCACGCTGCAGGTCGTACAGGCGCGTGAGCAGGGCGATCAAGGTCGACTTGCCGGCACCATTGGGGCCCAGCAATGCGGCGAAGCGCCCGGGTGCCAGGCTGAAACCGACCTGCCGGAGTGCCTCGCGCTGGCCATAGGCAAAGCTCACGTCACTGACGTCGAGGGCGTTCATGGTGTCACCACCACGCCCCAGGGGTAGCGCCCGACCTTCACCGACTTGAGCACTTTGAGCTGCTTCGCATCGATGACCGAAACGTCACCGCTGACACCGTTGGTGG  >  NC_002947/3056173‑3056467

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