Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,550,444 A→C 7.2% N98T (AAC→ACC)  PP_3134 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,550,4440AC 7.2% 159.1 / 4.3 56N98T (AAC→ACC) PP_3134hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/29);  new base C (2/2);  total (25/31)
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.17e-01

TCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCTTCATCGGCCATGGTGCGATCGTGATGCC  >  NC_002947/3550303‑3550591
                                                                                                                                             |                                                                                                                                                   
tCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGTTACCGATACGTACCAGCgac                                                                                                                                              >  1:172276/1‑147 (MQ=255)
 cTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTg                                                                                                                                              <  6:213657/148‑1 (MQ=255)
  tACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGTTACCGATACGTACCAGCgacgg                                                                                                                                            >  7:66385/1‑145 (MQ=255)
   tctgtaggTTCTGCAACCCCTCTGGTCTTGAATTGGGGGCCTTGCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCTATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAAGAGCGTggg                                                                                                                                            <  2:94682/143‑1 (MQ=255)
    cAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGTTACCGATACGTACCAGCgacgg                                                                                                                                            <  8:66385/147‑5 (MQ=255)
    cAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGTTACCGATACGTACCAGCgacgg                                                                                                                                            <  2:172276/147‑5 (MQ=255)
    cAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTggg                                                                                                                                            <  5:11890/147‑1 (MQ=255)
    cAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGc                                                                                                                                           <  3:47124/148‑1 (MQ=255)
        aggTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTc                                                                                                                                      <  4:61209/149‑1 (MQ=255)
          gTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCGTTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGTTACGTATCGGTAACAGCGTGGGCTTGTcc                                                                                                                                     >  2:11796/1‑148 (MQ=255)
          gTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCg                                                                                                                                    >  4:6217/1‑149 (MQ=255)
          gTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATATCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTcc                                                                                                                                     >  6:82718/1‑148 (MQ=255)
             cTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCACTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGAc                                                                                                                                  <  8:74073/148‑1 (MQ=255)
              tGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTAGGTCGGAAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACACGTCGCTGGTAGGTATCGGTAACAGCGTGGGCTTGTCCGCCTg                                                                                                                                >  4:91891/1‑149 (MQ=255)
               gCAACCCCTCTGGTCTAGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTg                                                                                                                                >  5:77173/1‑148 (MQ=255)
                    cccTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGCCCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTACGACTGCACGtt                                                                                                                          >  8:100572/1‑149 (MQ=255)
                    cccTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGtt                                                                                                                          >  7:62189/1‑149 (MQ=255)
                     ccTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGtt                                                                                                                          >  8:137886/1‑148 (MQ=255)
                           ttCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATc                                                                                                                      <  2:185471/145‑1 (MQ=255)
                            tCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTc                                                                                                                  <  3:91891/149‑1 (MQ=255)
                              ttGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCAt                                                                                                                >  3:20594/1‑149 (MQ=255)
                               tGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGc                                                                                                                                           >  6:65558/1‑121 (MQ=255)
                               tGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGc                                                                                                                                           <  5:65558/121‑1 (MQ=255)
                                gAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGt                                                                                                                                               <  6:263104/116‑1 (MQ=255)
                                gAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGt                                                                                                                                               >  5:263104/1‑116 (MQ=255)
                                gAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATg                                                                                                               <  7:137886/148‑1 (MQ=255)
                                 aaTGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGAc                                                                                                             >  1:20020/1‑149 (MQ=255)
                                  aTGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACg                                                                                                            <  5:82718/149‑1 (MQ=255)
                                       gggCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGggtg                                                                                                       <  7:136957/149‑1 (MQ=255)
                                       gggCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGggtg                                                                                                       <  3:6217/149‑1 (MQ=255)
                                        gtccTATCTGGCCAGATGGGGCAACTTTCCATCCGTCGGCTGCAATTTCTACGTCACTACAGGCTGCATGGTTCTCGACCCGTCGCTGGTACGTAGCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGggtgg                                                                                                      <  8:93008/147‑1 (MQ=255)
                                                 ggCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGTCTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGa                                                                                             <  6:30242/149‑1 (MQ=255)
                                                  gCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGTCCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGa                                                                                             <  7:100572/148‑1 (MQ=255)
                                                   cccGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGAt                                                                                            <  8:62189/148‑1 (MQ=255)
                                                                              gcagcaATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGAt                                                                >  2:255781/1‑149 (MQ=255)
                                                                                atcaATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATtc                                                              <  6:77173/147‑1 (MQ=255)
                                                                                            gTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTg                                                                                                 >  5:268635/1‑102 (MQ=255)
                                                                                            gTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTg                                                                                                 <  6:268635/102‑1 (MQ=255)
                                                                                            gTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCAt                                                  >  2:61181/1‑149 (MQ=255)
                                                                                                      gCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAg                                        <  1:255781/149‑1 (MQ=255)
                                                                                                      gCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAg                                        <  1:61181/149‑1 (MQ=255)
                                                                                                           aaGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACaa                                   >  5:120016/1‑149 (MQ=255)
                                                                                                              tttCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTg                                >  7:103977/1‑149 (MQ=255)
                                                                                                                tctcGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGc                               <  2:20020/148‑1 (MQ=255)
                                                                                                                  tcGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATc                                                 <  8:69466/128‑1 (MQ=255)
                                                                                                                  tcGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATc                                                 >  7:69466/1‑128 (MQ=255)
                                                                                                                    gACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCGAGGACAATTGCTTCAt                          >  1:243195/1‑149 (MQ=255)
                                                                                                                      cccGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCGAGGACAATTGCTTCATCg                        <  2:243195/149‑1 (MQ=255)
                                                                                                                      cccGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCTTCATCg                        >  5:35549/1‑149 (MQ=255)
                                                                                                                         gTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCTTCATCGGc                      >  6:78749/1‑148 (MQ=255)
                                                                                                                                gTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCTTCATCGGCCATGGTg               <  6:35549/148‑1 (MQ=255)
                                                                                                                                  aCGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCt                              <  4:169596/131‑1 (MQ=255)
                                                                                                                                  aCGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCt                              >  3:169596/1‑131 (MQ=255)
                                                                                                                                          gTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCTTCATCGGCCATGGTGCGATCGTGATg    <  8:197152/149‑1 (MQ=255)
                                                                                                                                          gTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCTTCATCGGCCATGGTGCGATCGTGATg    >  7:197152/1‑149 (MQ=255)
                                                                                                                                           tAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCTTCATCGGCCATGGTGCGATCGTGATGc   <  4:20594/149‑1 (MQ=255)
                                                                                                                                            aaCAGCGTGGGCTTGTCCGACTGCACGTAGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCTTCATCGGCCATGGTGCGATCGTGATGcc  >  1:151780/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCTTCATCGGCCATGGTGCGATCGTGATGCC  >  NC_002947/3550303‑3550591

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