Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,147,973 G→C 8.0% S32R (AGC→AGG)  PP_1006 ← heme receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,147,9730GC 8.0% 153.4 / 5.3 50S32R (AGC→AGG) PP_1006heme receptor
Reads supporting (aligned to +/- strand):  ref base G (17/29);  new base C (2/2);  total (19/31)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

AAGCCGACCTGCCAGCCGGTGACCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGCTGTTCTCGT  >  NC_002947/1147839‑1148118
                                                                                                                                      |                                                                                                                                                 
aaGCCGACCTGCCAGCCGGTGACCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTg                                                                                                                                    <  2:902069/150‑1 (MQ=255)
              gCCGGTGACCTGGCTGAAGGCGTTGATGGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCACGCCCGGGCTGCTCGCTTGGCCTGGCGGGGTGGCCAGCGCCGTCAGGCTGCATTCGCCGCTGGCGAGCAGGGCAAg                                                                                                                      <  2:137474/150‑1 (MQ=255)
               ccGGTGACCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAgg                                                                                                                     <  1:940561/150‑1 (MQ=255)
               ccGGTGACCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAgg                                                                                                                     <  2:1119454/150‑1 (MQ=255)
                    gACCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCaa                                                                                                                <  1:931311/150‑1 (MQ=255)
                     aCCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCTGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAAc                                                                                                               <  2:466301/150‑1 (MQ=255)
                       cTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAg                                                                                                             >  1:927454/1‑150 (MQ=255)
                       cTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAg                                                                                                             <  1:810934/150‑1 (MQ=255)
                           cTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGAc                                                                                                         <  2:364018/150‑1 (MQ=255)
                              aaGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGa                                                                                                                                <  2:273448/124‑1 (MQ=255)
                              aaGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGa                                                                                                                                >  1:273448/1‑124 (MQ=255)
                                 gCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCTGCAGCCTGccg                                                                                                                                             <  2:632848/108‑3 (MQ=255)
                                 gCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCTGCAGCCTGccg                                                                                                                                             >  1:632848/1‑106 (MQ=255)
                                  cGTTGATCGACGACACCACCGGCTTCTGAACGATGGTGAAGAGGTAGTCGCCCATTCGCGGGCTACTGGCCTGAGGCGGCTCGTCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTc                                                                                                  <  2:73741/150‑1 (MQ=255)
                                    ttGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCt                                                                                                <  1:887514/150‑1 (MQ=255)
                                       aTCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGGCAACAGGGACAGTTGTCCTGTa                                                                                             <  2:418183/150‑1 (MQ=255)
                                                acCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGa                                                                                    >  1:330706/1‑150 (MQ=255)
                                                 ccAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCCGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGaa                                                                                   <  1:378791/150‑1 (MQ=255)
                                                  cAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAg                                                                                  <  1:614969/150‑1 (MQ=255)
                                                  cAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAg                                                                                  >  2:1110838/1‑150 (MQ=255)
                                                    gCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGGCAACAGGGGCAGTTGTCCTGTACGGCgggtcaggg                                                                                >  2:142470/1‑144 (MQ=255)
                                                      ggCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGATTGGCGCGGCCCGGCAGCCGCCTCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTg                                                                              <  1:142470/150‑1 (MQ=255)
                                                           ctgGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGCGCAAGGCAAGGGTCTAAAAGGCAAGTTGGCCGGGAGGGCGGGGAAGGTCGGAcca                                                                         >  2:197464/1‑148 (MQ=255)
                                                                cGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAg                                                                                                                            <  1:846942/94‑1 (MQ=255)
                                                                cGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAg                                                                                                                            >  2:846942/1‑94 (MQ=255)
                                                                cGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCggg                                                                    <  1:160818/150‑1 (MQ=255)
                                                                  aTGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTa                                                                  <  1:9347/150‑1 (MQ=255)
                                                                     cTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGa                                                               >  2:1077201/1‑150 (MQ=255)
                                                                      tGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAgg                                                                                 <  2:630436/131‑1 (MQ=255)
                                                                      tGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAgg                                                                                 >  1:630436/1‑131 (MQ=255)
                                                                       tcaGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAg                                                                 >  1:948029/3‑146 (MQ=255)
                                                                       gAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAg                                                                 <  2:948029/146‑1 (MQ=255)
                                                                             gTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGa                                                       >  2:1084562/1‑150 (MQ=255)
                                                                             gTAGTCGCACATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGTGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGa                                                       >  2:565754/1‑150 (MQ=255)
                                                                                    cccATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGa                                                >  2:1022113/1‑150 (MQ=255)
                                                                                     ccATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGaa                                               <  2:229774/150‑1 (MQ=255)
                                                                                      cATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTAcggcgg                                                                                       >  2:255387/1‑109 (MQ=255)
                                                                                      cATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTAcggcgg                                                                                       <  1:255387/109‑1 (MQ=255)
                                                                                         gcgcGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAgcgc                                           <  2:927454/150‑1 (MQ=255)
                                                                                            cgGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTcc                                        <  1:1186148/150‑1 (MQ=255)
                                                                                                 tgctgGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACTGCGGGTAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGAc                                   <  2:295816/150‑1 (MQ=255)
                                                                                                                    ggcggcCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAgg                  <  1:77311/148‑1 (MQ=255)
                                                                                                                    ggcggcCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAgg                  >  2:77311/1‑148 (MQ=255)
                                                                                                                        gcCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATgg            <  1:910477/150‑1 (MQ=255)
                                                                                                                         gcggcgcCTGCAGCCTGCCGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGc           >  1:272788/4‑150 (MQ=255)
                                                                                                                          cAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGCt          <  1:464928/150‑1 (MQ=255)
                                                                                                                           ggcgcCTGCAGCCTGCCGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGCTg         <  2:272788/149‑1 (MQ=255)
                                                                                                                            gcgcCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGCTgt        >  2:780552/1‑150 (MQ=255)
                                                                                                                                cGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGCTgttctc    <  1:1110838/150‑1 (MQ=255)
                                                                                                                                  gCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGCTgttctcgt  >  2:363836/1‑150 (MQ=255)
                                                                                                                                      gCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAg                                              >  2:462007/1‑102 (MQ=255)
                                                                                                                                      gCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAg                                              <  1:462007/102‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                 
AAGCCGACCTGCCAGCCGGTGACCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGCTGTTCTCGT  >  NC_002947/1147839‑1148118

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑