Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,550,435 G→C 7.5% R95P (CGT→CCT) ‡ PP_3134 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,550,4350GC 7.5% 164.1 / 4.1 53R95P (CGT→CCT) ‡PP_3134hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/27);  new base C (2/2);  total (24/29)
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.46e-01

TCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCTTCATCGGCCATGGTG  >  NC_002947/3550303‑3550578
                                                                                                                                    |                                                                                                                                               
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                                                                                                                                  >  5:263104/1‑116 (MQ=255)
                                gAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGt                                                                                                                                  <  6:263104/116‑1 (MQ=255)
                                gAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATg                                                                                                  <  7:137886/148‑1 (MQ=255)
                                 aaTGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGAc                                                                                                >  1:20020/1‑149 (MQ=255)
                                  aTGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACg                                                                                               <  5:82718/149‑1 (MQ=255)
                                       gggCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGggtg                                                                                          <  3:6217/149‑1 (MQ=255)
                                       gggCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGggtg                                                                                          <  7:136957/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                                                                                    <  6:268635/102‑1 (MQ=255)
                                                                                            gTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTg                                                                                    >  5:268635/1‑102 (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                 >  3:169596/1‑131 (MQ=255)
                                                                                                                                  aCGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCt                 <  4:169596/131‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                               
TCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAGCAATTTCTACGTCAATACAGGCTGCAAGTTTCTCGACCCGTCGCTGGTACGTATCGGTAACAGCGTGGGCTTGTCCGACTGCACGTTGATCGGTCATGACGGGGTGGTATTGCTGATCGAGCATCGCTTTGGCAAGCACCTGGATTCGGTCGGTTTCATCGACATCAAGGACAATTGCTTCATCGGCCATGGTG  >  NC_002947/3550303‑3550578

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