Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,008,318 C→T 7.7% S259S (TCG→TCA)  PP_2628 ← ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,008,3180CT 7.7% 144.2 / 2.2 52S259S (TCG→TCA) PP_2628ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (22/26);  new base T (2/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

TCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATGAA  >  NC_002947/3008176‑3008464
                                                                                                                                               |                                                                                                                                                   
tCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCgg                                                                                                                                               >  1:402738/1‑149 (MQ=255)
  aCCAGGCGCATCTGGCGAAATCCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCggg                                                                                                                                              >  5:170041/1‑148 (MQ=255)
     aGGCGCCTCTGGCGAAAACCGGTCTGCTTGAAATCATCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCAtt                                                                                                                                          <  6:316350/149‑1 (MQ=255)
       gcgcATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCAGTgg                                                                                                                                        >  4:96713/1‑149 (MQ=255)
              tGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GATc                                                                                                                                >  5:427915/1‑149 (MQ=255)
                gCGACAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTc                                                                                                                                <  2:412049/147‑1 (MQ=255)
                gCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCt                                                                                                                              >  3:136697/1‑149 (MQ=255)
                   aaaaCCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCggg                                                                                                                                              >  1:346562/1‑131 (MQ=255)
                   aaaaCCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCggg                                                                                                                                              <  2:346562/131‑1 (MQ=255)
                       ccGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGc                                                                                                                            <  8:8246/144‑1 (MQ=255)
                       ccGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGc                                                                                                                            >  7:8246/1‑144 (MQ=255)
                        cGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGAcc                                                                                                                      <  2:233032/149‑1 (MQ=255)
                        cGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGAc                                                                                                                       <  6:170041/148‑1 (MQ=255)
                               cTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCAGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCa                                                                                                               >  4:158982/1‑149 (MQ=255)
                                ttGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAg                                                                                                              >  6:18257/1‑149 (MQ=255)
                                  gAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGc                                                                                                             <  8:184951/148‑1 (MQ=255)
                                       cTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGCGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCAGCCCGACCTTGCCCAGCCgggg                                                                                                        >  5:232700/1‑148 (MQ=255)
                                        ttCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGa                                                                                                      >  5:228192/1‑149 (MQ=255)
                                        ttCCATGGCTGCGCGCGCAAACTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGGGGGTATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCTAGCTCCTGCCCGACCTTGCCCAGCCggggg                                                                                                       <  1:254327/149‑1 (MQ=255)
                                          ccATGGCTGCGTGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAACACCCGCCCGGTCATCGACGACCGGGCATTgg                                                                                                                                        <  1:45969/115‑1 (MQ=255)
                                          ccATGGCTGCGTGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTgg                                                                                                                                        >  2:45969/1‑114 (MQ=255)
                                                cTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATAc                                                                                              <  4:136697/149‑1 (MQ=255)
                                                 tGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATa                                                                                               >  3:42236/1‑147 (MQ=255)
                                                 tGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATa                                                                                               >  7:169228/1‑147 (MQ=255)
                                                    gcgcgcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCa                                                                                          <  2:402738/149‑1 (MQ=255)
                                                     cgcgcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCAt                                                                                         <  8:169228/149‑1 (MQ=255)
                                                      gcgcAAAGTGCGCTTCGCGGGTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATg                                                                                        <  2:334616/149‑1 (MQ=255)
                                                                cgcTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCg                                                                               >  7:320147/1‑148 (MQ=255)
                                                                 gcTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGaa                                                                             <  6:120618/149‑1 (MQ=255)
                                                                 gcTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGaa                                                                             <  5:18257/149‑1 (MQ=255)
                                                                   ttCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACa                                                                           >  8:225529/1‑149 (MQ=255)
                                                                       cgGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCAc                                                                       >  1:192112/1‑149 (MQ=255)
                                                                                           ttttACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAggg                                                     >  2:259436/1‑147 (MQ=255)
                                                                                                ccGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATc                                              <  4:42236/149‑1 (MQ=255)
                                                                                                        aTGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCa                                       <  3:158982/148‑1 (MQ=255)
                                                                                                        aTGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCa                                       <  7:225529/148‑1 (MQ=255)
                                                                                                         tGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTGATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAgg                                     <  6:228192/149‑1 (MQ=255)
                                                                                                         tGTTGGGGTTAGCTTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GATCCTGCCCGACCTTGCCCAGCCGTGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAgg                                     <  6:232700/149‑1 (MQ=255)
                                                                                                            tGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGc                                   <  2:192112/148‑1 (MQ=255)
                                                                                                              gtgtAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCaaa                                <  1:198195/149‑1 (MQ=255)
                                                                                                               tgtAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAg                               <  6:417604/149‑1 (MQ=255)
                                                                                                               tgtAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAg                               <  8:320147/149‑1 (MQ=255)
                                                                                                                 tAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAgcc                             >  8:10045/1‑149 (MQ=255)
                                                                                                                   aGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGTACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAgccgc                           <  1:328307/149‑1 (MQ=255)
                                                                                                                   aGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAAAAAAAGATGGCCAACTGGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAgccgc                           >  1:360018/1‑149 (MQ=255)
                                                                                                                           ggccaatgcccgGTCGTCGATGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAgccgcc                          >  8:20965/8‑143 (MQ=255)
                                                                                                                           ggccaatgcccgGTCGTCGATGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAgccgcc                          <  7:20965/136‑1 (MQ=255)
                                                                                                                                  gcccgGTCGTCGATGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAgg                                     >  5:420831/1‑125 (MQ=255)
                                                                                                                                  gcccgGTCGTCGATGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAgg                                     <  6:420831/125‑1 (MQ=255)
                                                                                                                                        tCATCGACGACCGGCCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGACCAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTACAGGGCAAAGCCGCCGGGCAGGGCCATGAgcgc        >  8:307623/1‑148 (MQ=255)
                                                                                                                                          aTCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGa      <  7:10045/148‑1 (MQ=255)
                                                                                                                                           tCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGAt     >  5:111105/1‑148 (MQ=255)
                                                                                                                                             gacgacCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATGaa  <  7:237106/149‑1 (MQ=255)
                                                                                                                                             gacgacCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATGaa  <  3:96713/149‑1 (MQ=255)
                                                                                                                                             gacgacCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATGaa  >  1:156108/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAA‑ACCCGCCCGGTCATCGACGACCGGGCATTGGCCT‑GCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATGAA  >  NC_002947/3008176‑3008464

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